NGENIEUR-ARCHIV 


a CESELLSCHAFT FUR ANGEWANDTE MATHEMATIK UND MECHANIK 
; ZUSAMMEN MIT — 


_ A. BETZ- K. KLOTTER. E. METTLER - K. vy. SANDEN 
F. SCHLEICHER - E. SCHMIDT - E. SORENSEN 


. | HERAUSGEGEBEN VON. 
-R. GRAMMEL 
| XXIV.BAND | SECHSTES (SCHLUSS-) HEFT 1956 


a 


| SPRINGER-VERLAG - BERLIN / GOTTINGEN / HEIDELBERG 


|| Abgeschlossen am 20. November 1956 


! Postverlagsort Berlin Bees ear ae 


Be a i eet 


Das 


INCENIEUR-ARCHIV ee 


erscheint nach MaSgabe des eingehenden Materials zwanglos in einzeln bovetliattea Heften, die zu ‘Banden 
vereinigt werden. E 


Die fur das Ingenicur- Archiv. bestimmten Manuskripte sind unmittelbar an den Herausgeber 
Herrn Professor Dr. R: Grammel, Stuttgart N, Robert-Bosch-StraBe 101 ‘ 
oder an die Herren 
Professor Dr. A. Betz, Gottingen, Herzberger LandstraBe 39 A 
Professor Dr.-Ing. K. Klotter, Stanford (Calif.), 621 Alvarado Row. 5 
Professor Dr. E. Mettler, Karlsruhe-Durlach, Geigersberger Str. 12. 5 “2 ae 
. Professor Dr.-Ing. K. v. Sanden, Karlsruhe- West, HertzstraBe 16 af: H. West) ae | 


“Professor Dr.-Ing. F. Schleicher, Dortmund, Plauener Str. 44 a ae 

Professor Dr.-Ing. E. Schmidt, Technische Héchsebule: Miinchen, iRcibetE. Zhes Biaieg ic 

Professor Dr. -Ing. KE. Sérensen, Augsburg, Man Tee: : 
einzusenden. 


Die zum Druck angenommenen Arbeiten werden, soweit dies drucktechnisch moglich ist, nach der — 
‘ Reihenfolge ihres Eingangs beim Herausgeber verdéffentlicht. ‘ 


Die Mitarbeiter crhalten von ihrer Arbeit -zisammen 75 Sonderdrucke unentgeltlich. 


Fiir die Abfassung der Arbeiten wird auf das vom Deutschen Normenausschu8 herausgegebene 
Heft ,,€estaltung technisch-wissenschaftlicher’ Veroffentlichungen™ hingewiesen. Die Vorlagen fur Ab- 
bildungen sind auf besonderen Blattern erwiinscht und kénnen entweder in Reinzeichnungen oder in klar- —_ 
verstindlichen Handskizzen bestehen; die Beschriftung und nétigenfalls die eso erie a nimmt der — 
Verlag vor. ; 


Mit der Aunahme des Manuskriptes erwirbt der Verlag das ausschlieBliche Verlagsrecht fiir alle Sprachen : 
und Lander. : 


Im ,,ingenieur-Archiv erscheinende Arbeiten diirfen vorher an anderer Stelle atch veroffentlicht — 
sein und Such spater nicht anderweitig veréffentlicht werden. 


Photographische Vervielfaltigungen, Mikrofilme, Mikrophote von ganzen Heften, einzelnen Beitragen A 
oder Teilen dardus sind ohne ausdriickliche: Genehmigung des Verlages nicht gestattet. eS: 


$. 
Die Wiedergabe von Gebrauchsnamen, Handelsnamen, Warenbezeichnungen usw. in dieser Zeitschrift berethuigt = 
auch ohne besondere Kennzeichnung nicht zu der Annahme, daB solche Namenim Sinne der Warenzeichen und Marken- ~ 
schutz-Gesetzgebung als frei zu betrachten waren und daher von jedermann benutzt werden diirften. 


* 


SPRINGER-VERLAG - 
Heidelberg a : Berlin W 35 


Neuenheimer LandstraBe 28—30 / Fernsprecher 27901 ~ Reichpietschufer 20 / Fernsprecher Sammel-Nr. 249251 


: Inhalt: eae % 


Seite Ps 


Grammel, R., und H. Ziegler, Der kardanisch gelagerte schnelle symmetrische 
Kreisel mit Lagerreibung. Mit 8 Abbildungen . 


Trostel, R., Instationire Warmespannungen i in einer. Hohlkugel. Mis 4 Abbil- 
dungen 


Teer rman mr ae eee cen eg tat SS 373 
Volterra, E, The Equations of Motion for Curved ora Twisted Elastic Bars de- 
duced by the use of the .,Method of Internal Constraints“ SRE Rare Shes a 392. 
Zoller, K., Uber das Grammelsche Vertalira’ bei. Bigenschwingungsaulgaben Mit ; 
4, Abi oeauen ea Oe panes Cae trie ee ey ooo a RRS a et 401 i 


Paslay, P. R. and A. Slibar, The Motion of Automobiles in Unbanked Curves . 412 4 


Te 


INGENIEUR-ARCHIV 


XXIV. BAND SECHSTES (SCHLUSS-) HEFT 1956 
esse 


Der kardanisch gelagerte schnelle symmetrische Kreisel mit Lagerreibung 
Von R. Grammel und H. Ziegler 


1. Einleitung. Der EinfluB der Reibung in den Lagern eines kardanisch aufgehangten Kreisels 
auf seine Bewegung und also ihre Abweichung von der idealen Bewegung ohne Reibung ist bis jetzt 
nur unvollstandig abgeschatzt worden. Die bisherigen einfachen Ergebnisse! (Abnahme der Eigen- 
drehgeschwindigkeit, Zunahme der Prazessionsgeschwindigkeit und Sinken des Schwerpunkts 
beim schweren Kreisel) lassen sich in dieser Allgemeinheit nur zum Teil aufrechterhalten und 
miussen durch weniger einfache und in einigen Fallen sogar iiherraschende Ergebnisse ersetzt 
werden, wenn man bei einem Kardangehange, wie es in Abb. 1 fiir einen schweren Kreisel schema- 
tisch dargestellt ist, die bisher unbeachtete 
Koppelung zwischen der Reibung in den inne- | 


ren Kardanlagern AA und in den auf eren wl8 
Kardanlagern CC beriicksichtigt und die bis- s 3 = 
lang stets vernachlassigte Reibung in den SIRES 
mittleren Kardanlagern BB hinzunimmt. & 
Obwohl es der Rechnung schwerlich ge- Ali c 


lingen kann, den EinfluB der Reibung quanti- 
tativ genau zu erfassen, so wollen wir doch in 
der folgenden Untersuchung das Problem er- 
neut aufgreifen, mit dem Ziel, die Mannig- 
faltigkeit des Ablaufs der Bewegung wenig- 
stens qualitativ vollstandig wiederzugeben, 
und es somit im Rahmen bestimmter Voraus- 
setzungen theoretisch zu einem gewissen Ab- ess, 
schluB bringen. “Rrobnget. {8 

Wir beschranken uns hier auf den is 
schnellen symmetrischen Kreisel mit 
lotrechter auBerer Achse CC des als masse- 
los und starr angesehenen Kardangehanges, 
machen also die in der Theorie des schnellen 
Kreisels iiblichen Vereinfachungen (wenn auch 
spater ein wenig erweitert) und vernach- ee 
lassigen die Drehimpulsanteile der Kardan- ““% 
ringe. 

Die physikalische Art der Reibung in den 
Lagern, also die mathematische Form des 
Reibungsgesetzes, lassen wir so lange wie 
moglich offen, da ein wesentliches Ergebnis | 
unserer Untersuchung fiir alle Reibungsgesetze 
giiltig ist und das allgemeine Verhalten des ae) ar 

° . . Abb. 1. Rechtsdrehender Kreisel im Kardangehange mit den Achsen AA, 
Kreisels auch ohne Kenntnis des speziellen 33, Cc und den Reibungsmomenten M,, M,, und M, (bei positiven 
Reibungsgesetzes ziemlich weit verfolgt wer- Drehungen a, «,,, @p), dem Schweremoment My und dem Figendreh- 
den kann. Wo wir ein Reibungsgesetz be- tepals 
nétigen, um explizit rechnen zu kénnen, be- 
schrinken wir uns im wesentlichen auf zwei Typen von solchen Gesetzen, welche Idealisierungen 
der wirklichen Verhaltnisse sind, namlich eines fiir trockene Reibung und ein ziemlich allgemeines 
fiir fliissige Reibung. Die mit beiden Gesetzen erhaltenen Aussagen sind durch Versuche quali- 


tativ bestatigt worden (Ziff. 11). 
1 R.Grammel, Der Kreisel, Bd. J, S. 111, 2. Aufl. Berlin-Géttingen-Heidelberg: Springer 1950. 
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Ganz unberiicksichtigt lassen wir die kaum genau zu iiberblickende Abhangigkeit der Rei- 
bungsmomente von den Lagerkraften. Da beim schnellen Kreisel der Schwerpunkt nur schwach 
beschleunigt wird, so andern sich die Krafte in den Lagern BB und CC nicht merklich, wogegen 
das Reibungsmoment in den Lagern 4A durchaus yom Winkel 6 (Abb. 1) abhangen kénnte. Das 
Gegenmoment des Kreisels gegen das Reibungsmoment in den Lagern AA erzeugt in den Lagern 
CC ein Paar von seitlichen Kriaften, deren Starke vom Winkel 6 abhingt und also das Reibungs- 
moment in den Lagern CC in einer von 6 abhangigen Weise beeinflussen kann. Auch von diesem 


Einflu8 sehen wir ab. 


2. Bezeichnungen. Wir wihlen die Bezeichnungen wie in Abb. 1, benennen also die Figuren- 
achse e, die (waagerechte) Knotenachse k, die (lotrechte) Prizessionsachse p und die (in der 
lotrechten Querebene durch Figuren- und Prazessionsachse gelegene und zur Figuren- und 
Knotenachse senkrechte) Querachse q in iiblicher Weise. Abb. 1 gibt auch die positiven Rich- 
tungen dieser vier Achsen an, und wir setzen fest, daB der Schwerpunkt S des Kreiselkérpers, 
wenn er nicht mit dem Kardanmittelpunkt O zusammenfallt (wie beim schwerelosen, meistens 
, kraftefrei’ genannten Kreisel) auf der positiven Figurenachse liegt (sog. ,,schwerer‘‘ Kreisel). 

Die Zahlenwerte der Eigendrehgeschwindigkeit m, und der Prazessionsgeschwindigkeit w, 
(mit der sich also die Querebene und die Knotenachse um die Lotlinie drehen) zahlen wir beim 
rechtsdrehenden Kreisel (er ist in Abb. 1 dargestellt) positiv, wenn die zugehérigen Drehvektoren, 
jeweils mit der Rechtsschraubenregel gebildet, in der positiven Figurenachse und positiven Pra- 
zessionsachse liegen, beim linksdrehenden Kreisel, wenn sie in den negativen Achsen liegen; den 
Zahlenwert der Nickgeschwindigkeit w, = 6 zihlen wir bei beiden Kreiseln positivy, wenn der 
zugehérige Drehvektor dann in der positiven Knotenachse liegt. Fir die Zahlenwerte der zuge- 
hérigen Komponenten D,, D, und D, des Drehimpulsvektors sollen die gleichen Vorzeichen- 
Festsetzungen gelten, fiir die Zahlenwerte der zugehérigen Reibungsmomente M,, M, und M;, in 
den Achsen 4A, CC und BB jedoch gerade die umgekehrte Festsetzung. 

_ Der Vektor des Schweremoments M, liegt, wenn von Null verschieden, stets in der positiven 
Knotenachse und hat den Betrag 


M, = Qsino mit OE=cG. (2,1) 
wobei s = OS ist und G das Gewicht des Kreiselkérpers bedeutet. 


Wenn man die so definierten Vektoren des rechtsdrehenden Kreisels an einer Ebene parallel 
zur Querebene spiegelt und im Spiegelbild die Rechtsschraubenregel durch die Linksschrauben- 
regel ersetzt, so zeigt das Spiegelbild den gleichen Kreisel in Linksdrehung und mit umgekehrter 
Prazessionsdrehung w,, aber gleicher Nickdrehung @,. Weil die Eulersche Vektorgleichung der 
Bewegung invariant gegen Vertauschen jener beiden Regeln ist, so diirfen wir uns im folgenden 
durchweg auf die Behandlung des rechtsdrehenden Kreisels beschranken: alle Ergebnisse gelten 
dann ohne weiteres auch fiir den linksdrehenden, wenn wir die fiir beide Kreisel entgegengesetzte 
Festsetzung des positiven Zahlenwertes von wm, beachten und beim schweren Kreisel alle Aus- 
sagen tiber das Steigen und Sinken des Schwerpunkts unverdndert beibehalten. Wie werden 
dies spater auch noch an Hand der Bewegungsgleichungen bestatigen. 


3. Das Kriftespiel. Wir beabsichtigen, die Vektorgleichung der Bewegung des Kreisels 
D + [oD] =MN, (3,1) 


worin ® der Drehimpulsvektor, 0 der Vektor mit den Komponenten Wp und w, und Wt der Vektor 
des 4uBeren Momentes ist, in dem rechtwinkligen, die Drehungen ®, und wm, mitmachenden 
System aus Figurenachse e, Knotenachse k und Querachse qin Komponeneten zu zerspalten und 
miissen also vor allem die drei Komponenten M,, M, (auBer M,) und M, vont in diesem System 
finden. 
Die feste Umgebung ubt mittels der Lager CC ein Moment auf den auferen Kardanring aus, 
eee Se peter ig ES ee wee Querebene (Abb. 2), eine Komponente 
; nachse, und eine lotrechte Komponente M, 
zerlegt werden darf. Die Komponente M, ubertragt sich tiber die Lager BB und AA auf den 
Kreiselkérper und zerspaltet sich dort in eine Komponente Mj in der Querachse und in einen 
ersten Anteil M;’ zu dem Reibungsmoment M, in den Lagern AA. Die Komponente M, ist iden- 
tisch mit dem Reibungsmoment M, in den Lagern BB. Die Komponente M;j ist sdenniach mit 
dem Reibungsmoment M, in den Lagern CC und ubertragt und zerspaltet sich in eine Kom- 
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ponente M3 in der Querachse und in einen zweiten Anteil My zu dem Reibungsmoment M, in den 
Lagern 4A. Wenn wir zu den Vorzeichenregeln fiir M,, M, und M;, noch die Festsetzung hinzu- 


fiigen, da der Zahlenwert von M, positiv sein soll, wenn dieser Vektor in der negativen Quer- 
achse liegt, so gilt also 


M, = M, sin 6 + M, cos 6, | 


(3.2) 
M, = M,sind— M,cos 6. | 


Die drei Reibungsmomente M,, M, und M;, bedeuten, soweit Bewegungen @,, Wp, und w, 
iberhaupt vorhanden sind, die Gleitreibungsmomente und sind als solche durch die Lagerbe- 
schaffenheit bestimmt (und voraussetzungsgemaB nur durch sie, unabhaingig von der Raum- 
stellung des Kreisels). Man sieht sie bei ,,fliissiger“‘ Reibung als gegebene stetige Funktionen von 
We baw. wp baw. w, an, die gleichzeitig mit diesen Groen verschwinden, bei ,,trockener“ Reibung 
als vorgeschriebene Festwerte, bei ,,gemischter‘‘ Reibung, die wir aber im folgenden nicht be- 
handeln, als gegebene Funktionen von w, bzw. Wp baw. Wks 
die beim Verschwinden dieser Drehungen gegen von Null 
verschiedene Festwerte gehen. 

Die trockene Reibung (und auch die gemischte) kann so 
stark sein, da die eine oder beide Drehungen w, und w, 
gar nicht méglich sind oder im Laufe der Bewegung auf Null 
abgebremst werden. Wir sprechen dann von Sperrung der 
Figurenachse oder wenigstens einer ihrer Bewegungskompo- 
nenten; undin diesem Falle gehen die Gleitreibungsmomente 
M, und M, in die Haftreibungsmomente Mj und M;; iiber, 
von denen man nur wei, daB sie zwischen Null und einem 
Héchstwert liegen, den man mit den Gleitreibungsmomenten 
M, und M, anzusetzen pflegt. 

Die Reibungsmomente miissen gemaéfs unseren Vor- 
zeichenregeln den Forderungen 


VE = 0= sgn M, = sgno,. sgn M;, =sgna, (3,3) 


geniigen; und natirlich mu M, so klein sein, daB der Cha- 
rakter des Kreisels als eines schnellen Kreisels wahrend 
der untersuchten Bewegungsdauer nicht verloren geht, ge- 


schweige denn, da die Drehung ™, tiberhaupt gesperrt  Abb.2. Momentenzerlegung in der Querebene 
A 2 und senkrecht zu ihr. 
wurde. 


Me 


4, Die méglichen Typen der Bewegung. Wenn man in erster Naherung fiir den schnellen Kreisel 
den gesamten Drehimpuls D, wie iiblich, mit dem Eigendrehimpuls D, verwechselt (da dieser 
beim schnellen Kreisel die von wm, und a, herrithrenden Anteile yon Q weit iiberténen soll) und 
jetzt zu den Momentkomponenten M,, M;, und M, noch das Schweremoment M, (2,1) hinzu- 
nimmt, so zerspaltet sich die Vektorgleichung (3,1) in die drei Komponentengleichungen 


D, =—M., (4,1) 
D, w, sind = Q sind — M,, (4,2) 
DS yea 


q ? 


mit dem Wert von M,, der aus (3,2) durch Elimination von M, folgt; man findet 
M, sind = Mp, — M. cos 6. (4,4) 


Aus diesen Gleichungen kann man einige allgemeine Folgerungen ziehen, die noch von jedem 
speziellen Reibungsgesetz unabhangig sind. Wir stellen dabei vor allem fest, dafs diese Glei- 
chungen sich nicht andern, wenn man gleichzeitig D, und , gegen — D, und — w, vertauscht, 
wp, aber belaBt, da dann auch M, und M,, aber nicht M;, das Vorzeichen wechseln, und damit 
ist nochmals bestatigt, das man sich weiterhin auf den rechtsdrehenden Kreisel beschranken 
darf und damit auch das Verhalten des linksdrehenden kennt. 

Die Gleichung (4,1) driickt die erwartete Tatsache aus, da8 der Eigendrehimpuls D, und 
mit ihm die Eigendrehgeschwindigkeit «w, fortwahrend abnehmen, wenn ein widerstehendes 
Moment M, vorhanden ist. 
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Bei den Gleichungen (4,2) und (4,3) miissen wir offensichtlich weiterhin die beiden Falle 
sin) —0 ausscheiden, also erstens den sogenannten aufrechten oder stehenden Kreisel 6 = 0 
— sein Verhalten ist ein Stabilitatsproblem, das wir hier nicht behandeln wollen — und zweitens 
den hangenden Kreisel 6 = 180° — dieser Fall ist trival. Aber auch ein gewisser Bereich in der 
Umgebung von 6 = 0 und von 6 = 180° wird von unseren Naherungsgleichungen nicht erfaBt. 
Wie namlich schon beim reibungslosen Kreisel unsere Gleichungen die Nutationen der in Wirk- 
lichkeit i. a. nicht streng regularen, sondern nur pseudoregularen Prazession wp = Q/D, nicht 
enthalten, so verwischen sie auch bei Reibung die méglichen Nutationen der nunmehr (wenn 
Wp =F 0) als pseudoquasiregular zu bezeichnenden Prazession und kénnen also insbesondere auch 
die Nutationen, die in der Umgebung von 6 = 0 und 6 = 180° neben der Eigendrehung @, 
allein noch méglicherweise vorhanden sind, nicht wiedergeben. Wir wollen ein fiir alle Male 
verabreden, da wir mit der Aussage, die positive Figurenachse erreiche die Stellungen 6 = 0 
oder 5 — 180°, oder der Schwerpunkt steige bis zum héchsten oder sinke bis zum tiefsten Punkt, 
stets ausdriicken wollen, da®B der kleine Bereich um 6 = 0 oder 6 = 180°erreicht werde, der 
beim hinreichend schnellen Kreisel auch hinreichend eng und nicht oder kaum mehr wahr- 
nehmbar ist, in welchem die Bewegung der Figurenachse also nur noch aus Nutationen besteht, 
und dessen Einzelheiten unsere Theorie nicht mehr erfassen kann. 


Um nun weiter in méglichst allgemeiner Weise aus Gleichung (4,2) auf die Prazession und 
vermittels (4,3) sowie (4,4) auf die Nickbewegung zu schlieBen (wobei natizlich Koppelungen 
zwischen diesen Gleichungen zu beachten sind), sehen wir im Falle trockener Reibung zunachst 
von Sperrungen ab (mit denen wir uns erst in Ziff.7 und 8 beschaftigen werden). Im Falle fliissiger 
Reibung lassen wir alle Reibungsgesetze zu, bei denen M, und Mp monotone Funktionen von w, 
bzw. w, sind; ferner schlieBen wir reibungsfreie Lager (die ja in Wirklichkeit auch nicht vor- 
kommen) in dieser Ziffer noch aus. 


Um zuniachst den Fall des Kreisels ohne Schweremoment vorwegzunehmen, beachten wir, 
daB nach (4,2) mit Q = 0 die Vorzeichen von w, einerseits und M;,, mithin auch @, andererseits 
entgegengesetzt sind. Der Kreisel fiihrt also bei steigendem Schwerpunkt eine positive, bei 
sinkendem eine negative Prazession aus. Im ersten Falle (w, <0) ist mit w, auch M, positiv 
und nach (4,3) mit @, auch M, negativ, zufolge (4,4) also 6 < 90°. Im zweiten Fall (@, > 0) 
schlieBt man analog auf 6 > 90° und unter der Annahme, daB der Schwerpunkt weder steigt 
noch sinkt (@, = 0), auf 6d = 90°. Hiernach (und unter dem Vorbehalt, daf sich infolge von 
Sperrungen noch Modifikationen ergeben kénnen) steigt oder sinkt der Schwerpunkt unter dem 
KinfluB der Reibung bis in die héchste bzw. tiefste Lage, je nachdem er anfanglich héher oder 
tiefer liegt als der Drehpunkt, und nur bei anfanglich waagerechter Figurenachse bleibt diese 
stehen. 

Nachdem so der Kreisel ohne Schweremoment erledigt ist, setzen wir weiterhin Q > 0 voraus 
(um erst wieder in Ziff. 9 auf den Kreisel ohne Schweremoment zuriickzukommen). Aus Glei- 
chung (4,2) folgt dann, da®B fir M;, = 0 die Prazessionsgeschwindigkeit ®p mit derjenigen 


Ora 7 > 0 (4,5) 


des durchweg reibungsfreien Kreisels mit dem gleichen momentanen Eigendrehimpuls D, tiber- 
einstimmt. Diese ist nach unserer Verabredung positiv, und da D, infolge der Reibung in den 
inneren Kardanlagern AA allmahlich abnimmt, wachst @p in jedem Falle monoton an. Ferner 
folgt aus (4,2) und (4,5), daB fiir wo, <0, also M;, <0 die Prazessionsgeschwindigkeit Wp > Wp 
und fiir w, > 0, also M;, > 0 umgekehrt Wp <@p ist. Die (positive) Prazession erfolgt also 
unter dem Einflu§ der Reibung in den Lagern BB bei steigendem Schwerpunkt rascher, bei 
sinkendem langsamer als beim reibungsfreien Kreisel mit dem gleichen Eigendrehimpuls. Da- 


mit ist insbesondere Prazessionsumkehr oder iiberhaupt negative Prazession von Anfang an 
nur bei sinkendem Schwerpunkt méglich. . 


Wir fiihren jetzt mit M, > 0 das Reibungsmoment in den duBeren Kardanlagern CC ein, 
das zur Prazessionsgeschwindigkeit @, (4,5) gehért, und nehmen erstens an, da in einem be- 
stimmten Augenblick 

M,cos6 < M, (4,6) 


sei. Dies trifft stets zu, wenn der Schwerpunkt tiefer liegt als der Drehpunkt, bei héher liegendem 
Schwerpunkt dagegen nur fiir geniigend kleine Werte von M,. Nimmt man dann an, daB der 
Schwerpunkt nicht sinke, d.h. daB om, < 0 sei, dann ist auf Grund unserer Voraussetzungen 
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tiber die (trockene oder fliissige) Reibung wegen (W, = @, auch M, > M,, also nach (4,6) 


M,cosd < M,. (4,7) 


Aus (4,4) folgt daher M,> 0 und damit aus (4,3) w, > 0 im Widerspruch zu unserer Voraus- 
setzung. Der Schwerpunkt sinkt also momentan, und zwar sowohl beim rechts- als auch beim 
linksdrehenden Kreisel. Da dabei cos 6 abnimmt, wahrend M, nicht zu- und nach (4,5) Wp ZU-, 
mithin M, nicht abnimmt, so bleibt die Ungleichung (4,6) bestehen. Es ist daher klar, daB 
der Schwerpunkt (eventuelle Sperrungen stets vorbehalten) bis in die tiefste Lage sinkt (wobei 
die Prazession mindestens in der letzten Bewegungsphase miglicherweise negativ ist). 


Nimmt man zweitens an, daB in einem bestimmten Augenblick 
M,cos 6 = M, (4,8) 


sei, was nur bei hochliegendem Schwerpunkt und nur unter der Voraussetzung M, > M, moglich 
ist, dann fiihrt die gleiche SchluBweise unter der Voraussetzung w, < 0 auf das Ergebnis w, > 0 
und unter der Voraussetzung w, > 0 auf w, < 0. Die Annahme, da der Schwerpunkt steige 
oder sinke, ist damit als unzutreffend erwiesen, und es muB momentan w, — 0 sein. Bei trockener 
Reibung behalten M, und M, ihre Werte bei, und es entsteht also keine Nickbewegung. Dieser 
Zustand ist aber natiirlich labil — man denke an die méglicherweise vorhandenen Nutationen. 
Bei Fliissigkeitsreibung dagegen nimmt M, ab und M, zu; die Gleichung (4,8) geht dann in die 
Ungleichung (4,6) iiber, und der Schwerpunkt sinkt (falls nicht neve Sperrungen eintreten) wieder 
bis zur tiefsten Lage. 


Setzt man drittens voraus, daB in einem bestimmten Augenblick 
M, cos 6 > M, (4,9) 


sel, was wiederum nur bei hochliegendem Schwerpunkt und unter der Voraussetzung M, > M, 


zutreffen kann, dann folgt aus der Annahme «, > 0 wegen w, < w, baw. M, < M, nach (4,9) 
M, cos 0. M,., (4.10) 


und da hieraus nach (4,4) auf M, < 0 und daher mit (4,3) im Widerspruch zu unserer Annahme 
auf w, <0 geschlossen werden mu, kann der Schwerpunkt momentan nur steigen, und zwar 
wiederum sowohl beim rechts- als auch beim linksdrehenden Kreisel. Bei trockener Reibung 
behalten M, und M, ihre Werte bei; der Schwerpunkt steigt also bis in die héchste Lage. Bei 
fliissiger Reibung tut er dies, da M, ab- und M, zunimmt, nur unter der Voraussetzung, daf} 
die Ungleichung (4,9) nicht vorher in die Gleichung (4,8) iibergeht. Trifft dies aber zu, dann 
geht dieser Fall in den vorangehenden und damit das Steigen des Schwerpunkts vor dem Er- 
reichen der héchsten Lage in ein Sinken bis zur tiefsten iiber. 

Mit diesen Ergebnissen ist der Bewegungsverlauf im Falle fliissiger Reibung qualitativ vollig 
beschrieben, im Falle trockener Reibung bis auf etwaige Sperrungen. Im Hinblick auf die nun 
folgende Ausarbeitung der Einzelheiten sei nur noch vermerkt, daf der Winkel 6,, welcher die 
Bereiche momentanen Sinkens und Steigens des Schwerpunktes trennt, nach (4,8) durch 

M. 
cos 04 = M, (4,11) 
gegeben, also nur fir M, > M, reell, dann aber spitz ist. Da fiir diesen Grenzwinkel w, = 0 
und daher w, = ©, ist, kann man iibrigens statt (4,11) auch 
EMO) a (4,12) 
schreiben. 

5. Die Bewegung bei fliissiger Reibung. Der einfachste Ansatz fiir fliissige Reibung lautet mit 

drei nichtnegativen Festwerten €,, & und ¢, 


M, = € We » Mi, = Ek Wk » Mp = Ep Wp; (5,1) 


Er erfiillt die Forderungen (3,3). Mit diesen Werten und der Drehmasse (dem Tragheitsmoment) A 
des Kreiselkérpers um seine Figurenachse, also D, = Aw, gehen die Bewegungsgleichungen 


356 Grammel und Ziegler : Der kardanisch gelagerte Kreisel mit Lagerreibung Ingenieur-Archiy 


(4,1) bis (4,3) samt (4,4) tiber in 


: s 
ne | (5.2) 
D, wp sin 6 = Q sin 0 — Eh On» (533) 
D, @p 8in 6 = Ep Op — a D, cos 6. (5,4) 

Es ist zweckmafig, die letzten beiden Gleichungen umzuformen mit der Kote 
u = cos 0 (5,5) 
der sogenannten Kreiselspitze, als welche wir weiterhin einfach den Schwerpunkt ansehen wollen, 
da wir den Langenmafistab stets so wahlen kénnen, daB s = OS = 1 ist; wegen u = — 6 sin 6 

= —, sin 6 gehen (5,3) und (5,4) dann tiber in 

D, Wp (1 — v2) = Q (1— uv?) + &U, (5,6) 
D, is =—* Des — & Op » (5,7) 


und wir erinnern nochmals daran, daB diese Gleichungen fiir sind = 0, also u = + 1, keine 
exakten Aussagen mehr machen kénnen, aber immerhin bis in die allernachste Nahe von u = +1 
elten. 

; Die Integration von (5,2) ist zwar stets explizit moglich, diejenige der beiden anderen Glei- 
chungen (5,3), (5,4) baw. (5,6), (5,7) mit dem Wert D, aus (5,2) scheint es aber in geschlossener 
Form nicht zu sein, wenn é,, ¢ und ¢, alle drei von Null verschieden sind. Wir erledigen daher 
zuerst drei Sonderfalle. 

a) Der erste Sonderfall ¢, = 0 hat eine praktische Bedeutung fiir viele Kreiselapparate, 
z. B. fiir kardanisch aufgehangte Lotkreisel und kiinstliche Horizonte, falls diese mit Synchron- 
motoren angetrieben werden, deren Stander auf dem inneren Kardanring aufsitzt. In diesem 
Fall ist dauernd M, = 0, also ¢, = 0, da das Reibungsmoment in den Lagern AA durch das 
Gegenmoment des Motors aufgehoben wird, und dann behalt D, dauernd seinen anfanglichen 
(stationaren) Wert D, bei. Die Gleichungen (5,6) und (5,7) lauten jetzt 


Dy (1—w) =Q(1— wt) +i, (5,8) 
Dy u = — &) >. (5,9) 
Die Elimination von w, fiihrt auf 

[D5 (1 — uw?) + ex ep] u +epQ(1—u?) = 0. (5,10) 

Diese Gleichung hat die zu u = uy fiir t = 0 passende Lésung 
Ep Qt = De (uy — u) Sere, tn i eS : 
Weil gemaB (5,10) fiir e, > 0 und |u| <1 dauernd wu < 0 ist, so sinkt der Schwerpunkt in diesem 
Sonderfall von jeder Anfangskote uy (<1) aus dauernd [im Einklang mit Ziff. 4, da ja die 
Bedingung (4,7) stets erfiillt ist]; er erreicht, falls auch ¢, > 0 ist, seine tiefste Lage u = — 1 
erst mit tf = oo, also nur asymtotisch, dagegen, falls ¢, — 0 ist, nach der endlichen Zeit 

ep Var 40) De 
Ep Q : 

Im Falle ¢, = 0 bleibt der Schwerpunkt dauernd in seiner Anfangshéhe Uy [wiederum im Einklang 


mit Ziff. 4, da jetzt dauernd (4,8) erfiillt ist], auch dann, wenn ¢ > 0 ist, und zwar zweifellos 
jetzt stabil. 


Die Prazessionsgeschwindigkeit wird nach (5,9) mit dem Wert w aus (5,10) 


l-—-u? Q 
1 — wu? -+ (& &/D§) Do 


(5,11) 


r= (5,12) 


Opes (9513) 
es: & €p > 0 ist, nimmt Wp», wie man hieraus leicht abliest, von einer etwaigen positiven Anfangs- 
ote Uy aus zundchst zu, erreicht fiir die waagerechte Lage u = 0 der Figurenachse seinen 
Hochstwert 
DyQ 


DB + Ek oy (5,14) 


Op max = 
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und nimmt dann — oder tiberhaupt von Anfang an, wenn Uy < 0 war — monoton ab bis zu 
seinem asymptotischen Wert Null fiir die tiefste Schwerpunktslage u = —1. Falls jedoch & 
oder ep oder beide verschwinden, so behiilt Wp dauernd seinen Wert Q/D, bei, der ohne Reibung 
vorhanden ware. 

b) Der nachstwichtige Sonderfall ¢, =0 ist fiir schnelle Kreisel mit groBem Schwere- 
moment M,, also groBemQ, bedeutsam, weil bei solchen in der Regel die Prazessionsgeschwindig- 
keit w, zwar klein gegen w,, aber groB gegen jw, ist, so daB das Reibungsmoment M,, nicht wesent- 
lich sein mag. Wir diirfen hier ¢, > 0 voraussetzen (da der Fall ¢, = 0 schon erledigt ist) und 
wollen auBerdem ¢, > 0 voraussetzen (da wir den Fall ¢, > 0, e = 0 spater allgemeiner mit- 
behandeln werden). : 

Die Gleichung (5,2) kann man fiir sich integrieren. Man fiihrt dabei zweckmabig ein neues 
Zeitmafs ein, nimlich die ,,Exponentialzeit* 


&t/A 
Le 6), (5,15) 
welche mit ¢ monoton geht, und erhalt dann mit dem Anfangswert D, = D, fir t = 0, also 
T= T einfach 
re 
D.= = (5,16) 
und somit die erwartcte monotone Abnahiae des Eigendrehimpulses. 
Mit diesem Wert von D, und wegen e, = 0 folgt aus (5,3) 
Q 
Oo ——_——— 4h 5 Syslliee 
Daas (5,17) 
also eine monotone Zunahme der Prazessionsgeschwindigkeit. 
Diese Werte von D, und @, und der aus (5,15) folgende Wert 


oH a 


ee T (5,18) 
fiihren (5,7) tiber in 
du u 
—_——_ -—_- —— - —= 9 
it TF +t Ye = (5,19) 
mit der Abkiirzung 
epAQ é 
ug = = De (5,20) 


welche eine einfache Bedeutung hat. Sie geht namlich mit den Anfangswerten w,, = D)/A und 
pg = Q/Dy von w, und wy, iiber in UQ = Ep Wpo/Ee Wen, und das ist das Verhaltnis der beiden 
Anfangswerte M,, und M,, der Reibungsmomente M, und M_,, also 


M 
ug = M., . (5,21) 


Die GréBe ug ist somit fir M,) > Mp, die Kote u, = cos 0, des Schwerpunkts fiir den alsdann 
reellen Grenzwinkel 6, aus (4,12) im Augenblick t = 0, welcher die beiden Bereiche anfanglichen 
Sinkens und Steigens des Schwerpunkts trennt; fiir M.. < M,, hat ug nur die Bedeutung eines 
Parameters der Bewegung, da dann kein reeller Winkel 0, zu ug gehort. GemaB unserer Voraus- 
setzung €, > 0 ist nach (5,20) auch ug > 0. 

Das zum Anfangswert u = u, gehorige Integral von (5,19) lautet 


u.= (Uy + ug) T— ug T?, (5,22) 
woraus noch 

d ‘ , 

a7 = (uo + ug) —2u9 7 (5,23) 


folet. 
2 . « \e tf 
Wir miissen nun zwei Falle unterscheiden. Es sei erstens uy < ug. Dies bedeutet fiir ug < 1, 


da®B die Anfangskote u, unter dem wirklichen Wert ug liegen soll, fiir ug 21, da8B sie jeden 
beliebigen Wert (zwischen + 1 und — 1) haben darf. Jetzt ist schon fir T = 1, also bei Be- 
wegungsbeginn, zufolge (5,23) du/dT <0 und fiir T > 1 umsomehr du/dT < 0, und dies besagt, 


358 Grammel und Ziegler: Der kardanisch gelagerte Kreisel mit Lageireibung Ingenieur-Archiv 


daB der Schwerpunkt von seiner Anfangskote uy) an dauernd sinkt (im Einklang mit unseren 
Feststellungen in Ziff. 4). Die Exponentialzeit T_, bis zum Erreichen des tiefsten Punktes 
u = — 1 folgt aus (5,22) zu 


il; . 
a =r yl tug + Y(uo + wo)? 4 Fug]. (5,24) 


Sei zweitens uy > ug. Dies ist nur moglich, wenn 1> ug> 0 ist. Jetzt wird nach (5,23) 
du/dT > 0 fiir T = 1, und dies besagt, daB der Schwerpunkt zunachst steigt. Setzt man den 
Wert von T aus (5,23) in (5,22) ein, so kommt 


u = 7,,|(t - ug)? (in) | (5,25) 


und nun mu8 man zwei Unterfalle unterscheiden, je nachdem der aus (5,25) mit dujdT =P 
folgende Wert 


seeks VO) ea (5,26) 
4 ug 


V 


ist. Im Unterfalle u’ < 1 tritt eine Umkehr du/dT = 0 des anfanglich steigenden Schwerpunkts 
bei der Kote u’ wirklich ein, und es wird von da an (wiederum in Ubereinstimmung mit Ziff. 4) 
dauernd du/dT <0 [weil nach (5,23) von einem Augenblick mit du/dT =-0 an dieser Wert 
mit fortschreitender Zeit T nur abnehmen kann]. Die Exponentialzeit des Steigens ist nach (5,23) 
Ti ee (5,27) 
2 ug 

die Exponentialzeit des Sinkens folgt aus (5,24) mit u’ statt uy sowie demjenigen Wert von ug, 

der aus (5,20) mit dem zu T’ gehérenden Wert Dg des Eigendrehimpulses berechnet ist. 
Die Bedingung u’ < 1 fiir diesen Umkehrfall kann man nach (5,26) auch in der Form uy < u* 


schreiben, wobei 


u* = 2Vug — ug (5,28) 
cine nur von ug (5,20) abhangige charakteristische Kote ist, und zwar gilt stets 
Ug st acaeks (5,29) 


wie man rasch erkennt, wenn man darauf achtet, daB ug in (5,28) ein positiver echter Bruch ist. 
Die Bedingung dafiir, daB der Schwerpunkt von seiner Anfangskote uy aus zuerst steigt, bei 
einer Kote u’ < 1 umkehrt und bis zu seinem tiefsten Punkte u = — 1] sinkt, kann man somit 
in die Ungleichungskette 


Wor Ug =< wD) (5,30) 


zusammenfassen. Mit ug—> 1 geht iibrigens auch u* — 1, so daB dort der Bereich méglicher 
Anfangskoten uy mit spaterer Umkehr auf Null zusammenschrumpft. 
Im Unterfall u’ > 1 endlich steigt der Schwerpunkt ohne Umkehr bis zum héchsten Punkt 
uw = 1, wozu nach (5,22) die Exponentialzeit 
1 


a ug! 


i 


Uy + Ug V (uy ug)? 224 ug. (5,31) 
gehort. 

Hiermit ist eine in Ziff. 4 allgemein geschilderte Erscheinung an einem bestimmten Reibungs- 
gesetz auch rechnerisch bestatigt. 

c) Der letzte Sonderfall ¢, = 0 ist wohl der wenigst wichtige. AuGBer dem auch hier giiltigen 


Integral (5,16) fiir D, folgt jetzt aus (5,6) und (5,7) mit Ep = 0, dem Wert D, (5,16) und der 
Exponentialzeit (5,15) 


Wiese 1 E (5,32) 
und entweder 
O Ee€ Walas 
— a aig ‘ k 0 

ee Dy “AD, We uF eS) 

als Funktion der Zeit oder auch 

— ay fee, 

= p(t tag Toa) (to +9) (5,34) 


als Funktion der Kote u. 
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Die Deutung ist einfach. Je nachdem u, = 0, steigt der Schwerpunkt oder bleibt (labil) in 
. gleicher Hohe oder sinkt, und zum beschleunigten Hauptanteil QT/D, der Prazession wm, kommt 
noch ein ebenfalls beschleunigt positiver, verschwindender oder beschleunigt negativer Neben- 
anteil infolge der Reibung in den mittleren Kardanlagern BB. Fiir Uy > 0 behalt also die Pra- 
zession dauernd ihren Drehsinn bei; fiir uy < 0 dagegen wechselt sie ihre Drehrichtung nach der 
aus (5,33) folgenden Exponentialzeit 


Ee Ef 


1 Say f 
fie == ae (x — yl + 96 | mit OTe, 0° (5,35) 


_ welche wegen uy < 0 positiy ist, und zwar bei der aus (5,32) und (5,35) folgenden negativen Kote 
Up =x—Jl+x. (5,36) 


Falls e, > 0 ist, geht w,—> o fiir u > + 1, was nur wieder besagt, daB unsere Theorie bei 
-u = +1 nicht ohne weiteres anwendbar ist. Falls jedoch e, = 0 ist, bleibt @, wieder durchweg 
endlich. 

d) In dem nun noch iibrigbleibenden allgemeinen Falle ¢,> 0, ¢,> 90, €p > 0 kann man 
wenigstens einige qualitative Schliisse ziehen, ohne die Bewegungsgleichungen (5,6) und (5,7) 
- zu integrieren, was wohl nur numerisch méglich ware. 

Auch hier gilt das Integral (5,16) von D, mit (5,15). Wenn man sodann aus (5,6) und (5,7) 
das eine Mal wu eliminiert, das andere Mal @p und dabei noch (5,6), (5,18) sowie die Abkiirzung 
ug (5,20) beniitzt, so erhalt man die beiden Gleichungen 


D, 1 —w &, € \ 
pea 1S Sad eg ek ebeeoy 5,37 
hte BlsY BET), om 
(1 a i) TS — (u—ug T?) (1—v). (5,38) 
2 i 


Aus (5,37) folgt mit @, = 0, daB eine Umkehr der Prazession nur dann eintritt, wenn der 
Schwerpunkt sinkt, und zwar bei der negativen Kote u, (5,36). 

Aus (5,38) folgt, daB fir 
Uy Sug der Schwerpunkt 
dauernd sinkt, da mit wach- 
sendem T die rechte Seite von 
fs,30), die fir T—1 und 
Ug Sug Null oder negativ war, 
starker und starker negativ 
wird. Ug su* 

Was die Anfangskoten Ug=il" 
Uy > ug anlangt, so gibt Abb.3 
die Verhaltnisse schematisch Ug<Up<w* 
wieder. Dort sind jeweils zwei 


(u, T)-Diagramme der  Be- ig | eae : 
wegung vom gleichen An- 
fangswert uy fiir ¢, = 0 und Up<Ugy-———— _ 


€; > 0 miteinander verglichen. 
Fiir ¢, = 0 sind diese Diagram- 
me nach (5,22) Parabeln. Aus 
(5,38) schlieBt man, daf die 
Héchstwerte u’ von u sowohl 
fiir ¢, = 0, wie fiir ¢, > 0 auf 
der Parabel Abb. 3. Die (u, T)-Diagramme fiir 0 (genau) und fir &,> 0 (schematisch). 
bee Gee Me (5,39) 

liegen (in Abb. 3 strichpunktiert), und da die Kurve fiir ¢, > 0 bei gleicher Abszisse T links 
von der Parabel (5,39) jeweils eine geringere Steigung du/dT, rechts von ibr jeweils eine geringere 
Neigung — du/dT hat als die dortige Kurve fiir ¢, = 0. 

Fiir ¢, = 0 gehért, wenn die Umkehrkote wu’ > 1 ist, gemal (5,38) auch noch die Gerade 
u = 1 zur Lésung und stets die Gerade u = — I, und so kann fiir ¢, > 0 kein Zweifel iber den 
Verlauf der Kurven bestehen. 
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Fiir Anfangskoten u, mit der Bedingung (5,30) gibt es auch im Falle ¢, > 0 stets eine Umkehr 
des Schwerpunkts von anfanglichem Steigen zu endgiiltigem Sinken. Fiir Anfangswerte uy > u* 
ist die theoretische Kurve von e&, > 0 benachbart zu der Kurve von ¢, = 0, welche nun aus einem 
Parabelstiick und einer waagerechten Strecke u = 1 besteht, so daB also theoretisch tiber kurz 
oder lang eine Umkehr fiir ¢, > 0 zu erwarten ware. In Wirklichkeit existiert nun aber, wie 
schon erwahnt, bei u — 1 ein von Nutationen erfiillter Bereich, in welchem unsere Theorie nichts 
Exaktes mehr aussagen kann, und der in Abb. 3 durch die Schraffur angedeutet ist. Sobald 
eine Kurve in diesen Bereich eintritt, ist zu erwarten, daB die Figurenachse in ihm so lange 
festgehalten wird, bis der Kreisel seine Eigenschaft als schneller Kreisel verloren hat und die 
(fiir den Fall der Reibung noch unbekannte) Stabilitatsbedingung des nunmehr ,,stehenden“ 
Kreisels nicht mehr erfiillt ist. Deswegen gelten die Kurven ¢, = 0 fir u’ =1 nur, soweit sie 
in Abb. 3 voll ausgezogen sind, wogegen ihre dort gestrichelten Teile keine Exaktheit bean- 
spruchen kénnen. Fiir ¢, > 0 kann man, ehe eine strenge Theorie fiir die Umgebung von u = 1 
entwickelt ist, anhand von Abb. 3 qualitativ nur so viel sagen, dab fiir uy => u* und hinreichend 
kleine Werte von ¢, > 0 keine Umkehr des schnellen Kreisels eintritt, daB dagegen fiir hin- 
reichend groBe ¢, — so gro8 namlich, daB die Kurve nicht mehr in den schraffierten Bereich 
eindringt — Umkehr auch schon beim schnellen Kreisel méglich ist. 

In Abb. 3 sind auch die Exponentialzeiten T,, (5,31), T’ (5,27) und T_, (5,24) fur e, = 0 
eingetragen, und es ist durch den schraffierten Streifen bei u = — 1 angedeutet, dafs auch dort 
unsere Formeln keine exakten Werte mehr zu liefern vermégen. Die inneren Grenzen der beiden 
schraffierten Streifen sind natiirlich nicht scharf definiert, und das, was wir vorhin ,, Kindringen‘‘ 
in den Bereich bei u = 1 nannten, kann nur durch eine exakte Theorie wirklich festgestellt 
werden. Der Deutlichkeit halber sind die beiden Streifen in Abb. 3 viel zu breit gezeichnet. 

Man mufB hier schlieBlich noch bei einigen Formeln von Abschnitt b) eine kleine Einschrankung 
hinzufiigen, die um so unwesentlicher wird, je schneller der Kreisel ist, und zu welcher die zwar 
kleine, aber endliche Breite des schraffierten Bereiches bei u = 1 nétigt: Die strenge Grenze 
zwischen umkehrenden und nicht umkehrenden Bewegungen wu liegt nicht genau bei uy = u*, 
sondern ein wenig tiefer, und die diesbeziiglichen Ungleichungen wie (5,26) und (5,30) sind dann 
auch in dem Sinne aufzufassen, daB das Symbol ,, 2‘ soviel wie ,,merklich gréBer oder kleiner‘‘ 
bedeutet. 

Es erscheint zunachst befremdlich, da die Bewegung der Figurenachse beim Loslassen des 
Kreisels mit einer von Null verschiedenen Prazessionsgeschwindigkeit und i. a. auch mit einer 
von Null verschiedenen Geschwindigkeit wu zu beginnen scheint. Dies ist aber nur dadurch vor- 
getauscht, das wir die Nutationen beim schnellen Kreisel auBer acht gelassen haben. Tatsachlich 
beginnt natiirlich die Bewegung der Figurenachse ohne Ruck ganz stetig, spielt sich aber in den 
Nutationen sofort auf die Anfangswerte @, und u unserer Rechnung ein in dem Sinne, daB diese 
und alle Folgewerte die Mittelwerte der Bewegung sind, wie sie das Auge sicht, wenn es die 
Nutationen nicht beachtet, eine Erscheinung, wie sie ja auch sonst aus der Kreiseltheorie wohl- 
bekannt ist.+ 

6. Eine Verallgemeinerung. Die linearen Gesetze (5,1) fiir die fliissige Reibung stimmen mit 
Versuchen an Lagern nicht gut iiberein. Soweit die Reibung wirklich als fliissig und noch nicht 
als gemischt gelten kann, hat man bessere Erfahrungen mit allgemeineren Potenzgesetzen 


1 
IM Senne VA eee a Mp = €p 5 * (6,1) 


gemacht, wo é,, é, und €, nach wie vor nichtnegative Festwerte und I, m, n positive, nicht not- 
wendig ganze Zahlen sind. Merkwiirdigerweise lassen sich die Sonderfalle b) und c) von Ziff. 5 
mit diesen allgemeineren Reibungsgesetzen ebenso explizit durchrechnen; und obwohl wir uns 
dessen bewuft sind, daB unsere Theorie mehr qualitative als quantitative Richtigkeit anstrebt, 
so wollen wir doch die Ergebnisse auch fiir die allgemeineren Gesetze (6,1) hier zusammenstellen, 
— namentlich auch deswegen, weil explizit berechenbare Probleme mit nichtlinearen Reibungs- 
gesetzen ziemlich selten sind. 

In (6,1) und bei allen folgenden Potenzen mit nichtganzzahligen Exponenten soll stets der 
positive Hauptwert gemeint sein, und in jedem Falle sind natiirlich die Forderungen (3,3) nach- 
zuprifen, 

Die erste Bewegungsgleichung (4,1) lautet jetzt 

“ Das 


Ee 


ar D. (6,2) 


1) Siehe etwa R. Grammel, Der Kreisel, a. a. O. S. 76. 
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und gibt integriert zum Anfangswert D,, wenn wir weiterhin den schon in Ziff. 5 erledigten Fall 
1 =] ausschlieBen, 


De = T (6,3) 
mit der neuen ,,Exponentialzeit* 


1 
(l— 1) €e Di" , =i 


ee (6,4) 
wozu noch 
CIES RI) atl 
a eee i) 


gehért. Auch diese Exponentialzeit T geht von T = 1 fiir t = 0 aus und nimmt sowohl fiir 
l> 1, wie fiir] < 1 mit t monoton zu. Da sich weder der allgemeine Fall e, > 0, &, > 0, €, > 0, 
noch der Sonderfall a) von Ziff. 5 weiter explizit behandeln lassen, so wenden wir uns den beiden 
Sonderfallen b) und c) von Ziff. 5 zu. 

b) In dem Sonderfalle ¢, = 0 nehmen jetzt die Bewegungsgleichungen (4,2) und (4,3) als 
Verallgemeinerungen von (5,6) und (5,7) die Gestalt an 


D, w =0Q, (6,6) 
De ai Dia enor. (6,7) 

Aus (6,6) folgt mit (6,3) 
oa é T, (6,8) 


also wieder eine monotone Zunahme von w,, und damit ist verbiirgt, daB die Ansatze (6,1) in 
diesem Sonderfalle die Forderungen (3,3) erfiillen. 


Mit (6,5) und (6,8) geht (6,7) tiber in 


du u tat 
oe ee 1_ 9 (6,9) 
mit der Abkiirzung 
¢ Al Oz 
__ &p 
ug == “e, Diem ’ (6,10) 
welche sich ohne weiteres in 
Ep po __ Mp, Git 
ue e608, Moy (6, ) 


umformt, also die friihere Bedeutung hat. Da wir auch hier e, > 0 und ¢, > 0 voraussetzen 
diirfen, so ist ug > 0. Wir unterscheiden nun zwei Fille. 


1) 1 +n +1. Hier hat (6,9) das zum Anfangswert ug passende Integral 


1 " ee ¢ 
eee 1 — 1) Mo 7 uel P= ug Tt") oy) 
samt 
du 1 


i earl {Ll +n—1)uy + uo] —(l+n)ug T ae : s (6,13) 


Aus (6,13) folgt fiir den Bewegungsbeginn T = 1 


du 
=| = U)— Ug, (6,14) 
(a7) =“ 


also anfangliches Steigen oder Sinken des Schwerpunkts, je nachdem uy Z ug ist. Falls ug = 1 
ist, sinkt der Schwerpunkt von Anfang an und dann gemaB (6,13) dauernd bis zu seinem tiefsten 
Punkte. 
Falls ug < 1 ist und uy) < ug war, so sinkt der Schwerpunkt ebenfalls monoton bis zu seinem 
tiefsten Punkt; dagegen steigt er fiir uy > ug zunachst an, und zwar héchstens bis zu der aus 
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(6,13) mit du/dT = 0 folgenden Zeit 


1 

,__[+n—1]) uy + ue iat ae 615 

ie eae , (6,15) 

wobei in der Tat T”’ > 1 fiir uy > ug ist. Jetzt kommt es wieder darauf an, ob die zu T” ge- 

hérende Umkehrkote u’ 2 1 ist; denn wenn u’ > 1 ist, bleibt wu vor oder mit der Zeit T’ beim 

Wert u =1stehen. Wenn jedoch u’ < 1 bleibt, so sinkt fir T > T’ der Schwerpunkt dauernd 

bis zum tiefsten Punkt. Wie dabei die Zeiten T,, und T_, aus (6,12) zu berechnen sind, brauchen 
wir nicht weiter zu erlautern. 


Die Kote wu’, die zu T’ gehért, folgt aus (6,12) mit (6,15) zu 


tH 
(hep te De se 80 Pe (6,16) 
(I + n)ug 
Die Bedingung wu’ < 1 fiir eine Umkehr fiihrt nach (6,12) auf 


1 


U Up Te up 
L 


Ven 


(i-+n—l)uS(l+n)ug™"—ug (l+n2Z1) 


oder (wegen uy > ug fiir anfangliches Steigen) auf 


ug <Uy < ui fir J+n>1, (6,17) 
igh ug us mae) ean ea ed ; 


mit den Abkiirzungen 


mee ff Se Ss Ik. 


uy eee la n) ug. "— ug 
1 


(6,18) 
fur J+tnc<l, 


gl! ln 


und man stellt fest, daB fiir jeden positiven echten Bruch ug (wie es ja jetzt zutrifft) die eckigen 
Klammern in (6,18) positiv sind. Man priift auBerdem leicht nach, daB mit wachsendem ug 
auch uj und u3 monoton wachsen, und daf hier mit ug —> 1 ebenfalls uf —> 1 und ux — 1 gehen. 


2) 14+-n =1. In diesem Falle lauten die Ergebnisse (mit denselben Gleichungsnummern) 
“ = T (u—ug ln T), (6,12*) 
du 
IT=%— uw(l+mnT), (6,13*) 

du 
8) asa an 

ty Ug 
Teen Oe (6,15*) 
Wy ug 

Wi Up 1 ao € “e 4 (6,16*) 
Wola, <a (6,17*) 
u3 = ug (1 — In ug) (6,18*) 


mit gleichen Folgerungen wie im Falle ] + n = 1. : 
Zusammenfassend ist zu sagen, daf} alle wesentlichen Ergebnisse von Ziff. 5 b) auch fiir das 
allgemeinere Reibungsgesetz qualitativ erhalten geblieben sind — wie ja nach Ziff. 4 zu erwarten 
war. 
c) Im Sonderfalle e, = 0 gelten nach wie vor die Gleichungen (6,2) bis (6,5). Die beiden 
anderen Bewegungsgleichungen nehmen nun statt (5,3) und (5,7) die allgemeinere Gestalt an 


D, wp sin d = Q sin 6 $F e& |6|™ (520), (6,19) 


. Ee 2 
Dia ne, as (6,20) 
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womit wiederum die Forderung (3,3) erfiillt ist. Denn aus (6,20) folgt mit (6,3) und (6,5) 


ian pad Me (6,21) 


also ein monotones Steigen oder Sinken des Schwerpunkts, je nachdem die Anfangskote u, = 0 


war, und eine (labil) waagerechte Prazession der Figurenachse, falls diese anfanglich waagerecht 
war. 


Wegen u = cos 6, also § = — u/sin 6 hat man statt (6,19) mit (6,3) 
Dy ‘ : jum 4 
Pp Mpsin 6 =Q sin db + en Gam 9 (uz 0)" 


- oder wegen der aus (6,21) und (6,5) folgenden Beziehung 


é oD: ; 
REA acer yaa he 


und mit dem aus (6,21) folgenden Wert sin? 6 = 1 — ui T? 


Q 5 TA-Dm+1 
Got De Ta ie (Uy = 9) (6,22) 
(l==a2-72)2 ; 


mit der Abkiirzung 


ee NEUE dae 


a= Alm (6,23) 
oder auch 
_ Q |cosd |cos 6|@-Dm+1 
Op Ds | Uo | Se b sin™+1 A 4 (up Zz 0) (6,24) 


mit der Abkiirzung 


Se 


b Alm 


(6,25) 


In (6,22) ist w, als Funktion der Zeit, in (6,24) als Funktion des Winkels 6 dargestellt. 


Fir uy, > 0, also steigenden Schwerpunkt, behalt die Prazession ihren Drehsinn dauernd bei; 
fiir uy <0, also sinkenden Schwerpunkt, kann sie ihr Vorzeichen bei einem Winkel 6, > 90° 
wechseln, wenn der aus (6,24) mit w, = 0 zu ermittelnde Wert 6, reell ist. 

Auch fiir diesen Sonderfall sind also alle wesentlichen Ergebnisse von Ziff. 5 c) qualitativ er- 
halten geblieben. 


Wir wollen schlieBlich nur noch bemerken, da alle Formeln von (6,6) an auch dann gelten, 
wenn man wieder speziell ] = 1 wahlen will; man mu dann nur unter T die Exponentialzeit 


(5,15) statt (6,4) verstehen. 


7. Bewegungen und Sperrungen bei trockener Reibung. Wenn die Reibung in allen Lagern trocken 


ist, so sind in den Bewegungsgleichungen (4,1) bis (4,4) die Reibungsmomente M,, M; und M, 


feste Werte. Nimmt man diese Festwerte weiterhin als positiv an, so muf man wegen (3,3) die 
Bewegungsgleichungen in der Form schreiben 


Da Me (7,1) 
D, wm, sind = Q sin 6 FM, (Mp = 0), (72) 
D, w,sin 6 = + M, — M, cos 6 (@, z 0). (G3) 


Diese Gleichungen gelten nur fiir ungesperrte Bewegungen, und nur fiir ungesperrte Bewegungen 
sind M, und M,, die vorgegebenen Gleitreibungsmomente. 

Da beim schnellen Kreisel keine Sperrung der Eigendrehung «, in Betracht zu ziehen ist, so 
ist M, stets als vorgegeben anzusehen und kann (7,1) wieder fiir sich allein integriert werden; man 
erhalt, wenn der Fall M, = 0 zunachst ausgeschlossen wird, mit der neuen ,»,Exponentialzeit“ 


= (7.4) 
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wozu noch 


oer (7,5) 
0 
gehért, das zum Anfangswert D, passende Integral 
dD. = = ? (7,6) 


also, weil ja auch T monoton mit t geht, die erwartete monotone Abnahme von D,. 
Die Gleichung (7,3) nimmt, in die Kote u (5,5) und in die Zeit T umgeschrieben, wegen (7,6) 
die Gestalt 


wa ON) ie 
Dae =u-- uy (Mp = 9) (To 
an, wobei wir die Abkiirzung 
ig a (7,8) 


mit ihrer bei (4,12) erklirten Bedeutung beniitzen. Fiir manche Zwecke ist es besser, noch ein 
zweites Zeitmah 


os Do 
T= In Es == In D, Da AP (7,9) 
einzufiihren, ebenfalls monoton mit t, soda man (7,7) auch in der Gestalt 
d 
z= =u uy (Wp = 0) (7,10) 
schreiben kann. : 
Die Gleichung (7,2) schlieBlich geht mit 
‘= a (7,11) 
und (7,6) uber in 
ee Sp s > Pere . 
op = ET (1 sae! 4 (2 0;u = 0). (7,12) 
Dabei definiert fir 0 < M, <Q die Abkiirzung s, einen reellen Winkel 6, mit 
sin bp = ae (7,13) 


wogegen s, fir M, >Q nur ein Parameter ist. Wir miissen nun vier Falle unterscheiden. 

a) M,>M., M, >Q. Hier ist nach (7,8) und (7,11) u;>1 und s,>1, also 6, nicht reell, und 
nun fiihren die Gleichungen (7,10) und (7,12) auf einen Widerspruch. Denn falls u>0 ware, so 
mite nach (7,12) w,>0 sein, was aber nach (7,10) auf du/dr < 0, also auch wu <0 fiihren wiirde. 
Aus u< 0 wiirde in gleicher Weise u>0 folgen. Dieser Widerspruch ist nur zu beseitigen, wenn 
man! ees also u->—0 und wm,-—>+-0 gehen laBt und in (7,10) und (7,12) statt u, (7,8) und 
sp (7,11) die Parameter 


q M. ’ 
einfiihrt, was mit dem jeweils vorhandenen Wert u = cos 6 nach (7,10) und (7,12) die Haft- 
reibungsmomente 


My = Q sino, Via Ve, (7,14) 


ergibt und also vollige Sperrrung der Schwerpunktsbewegung wu und der Prazession w, bedeutet. 
Die Haftreibungsmomente Mf und M; sind nach Betrag und Vorzeichen einleuchtend: Mj 
mu das Schweremoment kompensieren, M; die Riickwirkung des Reibungsmomentes M,, als 
Vektor Mj? fiir u> 0 nach unten, fiir u<0 nach oben weisend. 
In einem Phasendiagramm mit den Koordinaten u und du/dr wird diese véllige Sperrung der 
Figurenachse durch das stark ausgezogene Geradenstiick in Abb. 4 dargestellt. Fir diesen Fall 
ist das Phasendiagramm allerdings noch trivial und nur der Vollstandigkeit halber angegeben. 


4 Dabei soll x > + 0 baw. x + —0 heifben: x geht von positiven bzw. von negativen Werten her gegen 
Null. 
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b) M,p<M., Mi>Q. du 
Jetzt ist uy <1 und at 
BL.’ Dies. fibrt fir 
Ug > w > —u, auf den 


gleichen Widerspruch Ope 
wie vorhin, fiir u> u, " im 
dagegen nur dann auf Abb. 4, Phasendiagramm fir M, > M,, M,,> Q. 

keinen Widerspruch, 

wenn w>Q0 und also ie 

Wp > 0 ist, fir u< tly. 


wenn u<0 und also 
Wp <0 ist. Im Phasen- 
diagramm gibt das die 
stark ausgezogenen Ge- 
radenstucke von Abb.5. {_——_— 
Dort stellt die obere Ge- ust 
rade die ,,obere Glei- op 
chung (7,10), die untere 
Gerade die  ,,untere‘‘ 
Gleichung (7,10) dar. ea 
Man hat volligeSperrung ae 
von u und Op im Bereich Abb. 5, Phasendiagramm fiir M,, < M,:; M,, Oe 
zwischen u, und — u,, 
und zwar mit den Haft- 

reibungsmomenten iWeb 
(7,14), dagegen Steigen at 
des Schwerpunkts mit 
positiver Prazession w, 


fir u> u, und Sinken 

des eee punks mit € ete 

negativer Prazession @, ae 

fir u < — uy. Se Uy woe: 
i MoM, <0. cement nS ae 

Nunist u,>1unds, <1, 


also 6, reell und die zu- 


| 
| 
| 
gehorige Kote u, =cos0, 
| 
| 


Ye 


<1 (und positiv). Man 
 pruft leicht nach, daB 
kein Widerspruch nur 
fir u<0 und u, =u ae 


=>—u, entsteht, dal per 


also der Schwerpunkt 
: . s Abb. 6. Phasendiagramm fir M, > M,, M), ale 
im Bereich zwischen u, P 


und —u, sinkt, und 

zwar mit positiver Pra- tau 1s] 
zession @,,wogegen ober- 
halb u, und unterhalb 
—u, Sperrung eintritt. 
Das Phasendiagramm 
zeigt Abb. 6. 

d) M,< M., M.<0. 
Hier ist uy <1 und s, <1, 
also 6, reell mit u,< 1, 
und nun hat man zwei 
Unterfalle. 

1) u, > up. Man er- 
halt jetzt ohne Wider- 
spruch u > 0 und Op > 0 Abb, 7, Phasendiagramm fur My < M,, M, < Qund Ug > Up: 


Ve 
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fiir u > u,, ebensou <0 und, > 0 fir up >u>— up, und endlich u <0 und w, <0 fur 
u <—vuy,, in den anderen Bereichen dagegen Sperrung. Das Phasendiagramm Abb. i zeigt 
dies iibersichtlich: Steigen des Schwerpunkts mit positiver Prazession. fiir u e uy, Sinken mit 
positiver Prazession im Bereich zwischen u, und —u, und mit negativer Prazession fiir u << — u,, 
sonst Sperrung. 

2) uy <u, Man stellt leicht fest, daB nun beide Sperrbereiche des vorigen Unterfalls ver- 
schwunden sind und das Phasendiagramm Abb. 8 entsteht. 

Zusammenfassend kann man aus Abb. 4 bis 8 entnehmen, daf fiir eine Anfangskote uy = uy 
nur (labiler) Stillstand der Figurenachse mdglich ist, fiir uy < uy auBer (gesperrtem) Stillstand nur 
Sinken des Schwerpunkts, fiir uy > u, nur Steigen des Schwerpunkts in Frage kommt, und da 


bei steigendem Schwerpunkt nur eine positive Prazession, bei sinkendem Schwerpunkt dagegen 
auch wieder ein Umschlag 


zu negatiyver Prazession még- 
ds lich ist. 
Die zur u-Achse normalen 
Teilstticke in den Phasen- 
diagrammen (in Abb. 6 bis 8 
gestrichelt) mit den Abszissen 
Up, —Up, stellen Spriinge in 
den Geschwindigkeiten u, w, 
dar. Diejenigen mit der Ab- 
szisse + u, entsprechen dem 
bei labiler Sperrung infolge 
einer kleinen Stérung  ein- 
setzenden Sinken des Schwer- 
Abb. 8. Phasendiagramm fiir Mp < M,, M,, < Qund Mg <p punktes. Diejenigen mit der 
Abszisse — u, treten bei sin- 
kendem Schwerpunkt von selbst auf (da u nicht konstant bleiben kann, solange u von Null ver- 
schieden ist) und haben (Abb. 6 und 7) eine plétzliche Sperrung oder aber (Abb. 8) den Umschlag 
von der positiven zur negativen Prazession zur Folge. 

Die Grenzfalle, in denen statt der Ungleichheits- die Gleichheitszeichen gelten, sind aus den 
Phasendiagrammen ohne weiteres abzulesen. 

Die Haftreibungsmomente Mj und M; fiir die Sperrungen in den Fallen c) und d 1) kann 
man analog zu (7,14) finden. 

Natirlich ist auch in den Diagrammen Abb. 4 bis 8 wieder ein schmaler Streifen je bei u = + 1 
von unseren Gleichungen nicht mehr streng erfaBbar, und fiir die scheinbar plétzlich einsetzenden 
Bewegungen, Sperrungen und Spriinge von u und q, ist auf die SchluBbemerkung von Ziff. 5 zu 
verweisen. 

Wahrend die Prazessionsgeschwindigkeit w, (soweit ungesperrt) schon durch (7,12) darge- 
stellt ist [mit Vorwegnahme des erst aus der folgenden Gleichung (7,15) zu entnehmenden Wertes 
0 als Funktion von T), so liefert (7,7) den zeitlichen Verlauf von u (soweit ungesperrt) zu 


u=+u,+ (ug uu) T (@, 20). (7,15) 
und hieraus folgt in den Fallen von Abb. 5, 7 und 8 die Exponentialzeit des Steigens bis zum 
héchsten Punkte u = 1 


= 


Ye 


TT. =——— (uo > u, <1) (7.16) 


und etwa im Fall von Abb. 8 die Exponentialzeit des Sinkens des Schwerpunktes bis zur Kote 


u=—U, der Prazessionsumkehr 
=i 


Py (uy < uy) (7,17) 


P Ug — Ug 
und dann vollends noch die Exponentialzeit des Sinkens des Schwerpunkts bis zum tiefsten 
Punkte u = — 1 


(ype pe) ee (1 — ug) (up + ug) 
P (up — tq) (ug — Ug) 


(Uy < Ug, pes Uy): (7,18) 


Wahrend die friiheren Sonderfalle M, = 0 bzw. M, = 0 schon in den bisherigen Diagrammen 
und Formeln mit s, = 0, u, = 1 bzw. u, = 0 enthalten sind, so bedarf der Sonderfall Vis) 
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des synchron angetriebenen Kreisels einer besonderen Behandlung. Jetzt sind die Bewegungs- 
gleichungen statt (7,1) bis (7,3) bzw. (7,12) 


D, = Dy, (7,19) 
= Q 1 Sp ss a 
Op = De fb sine (0 ee 0. u = 0) 5 (7,20) 
. My 
w= + oa (wo, 20). (721) 


Man stellt leicht fest, da8 fiir M, > Q, also Sp > 1, stets ein Widerspruch vorhanden ist, also 
vollige Sperrung eintritt, und da8 fiir M;, < Q, also Sp <1 nur bei up > u > —u, und u < 0 
mit Wp > 0 kein Widerspruch vorliegt, was Sinken des Schwerpunkts mit positiver Prazessions- 
drehung im Bereich zwischen u, und —u, bedeutet und Sperrung fiir u> up undu <—u ,. In der 
Phasenebene wird dies dargestellt durch Abb. 6, worin jetzt t statt 7 geschrieben werden darf. 
Fiir das Sinken des Schwerpunkts gilt nach (7,21) 


M 
U = Up a t (u, > u> —u,). (7,22) 


0 


8. Weitere Sperrungen. Aufer den in Ziffer 7 gefundenen Sperrungen bei trockener Reibung 
gibt es noch weitere Sperrungen, die aber aus unseren Bewegungsgleichungen (7,2) und (7,3) nicht 
entnommen werden kénnen. Dies riihrt davon her, daB diese Gleichungen nur eine erste Nahe- 
rung der exakten Vektorgleichung (3,1) fiir den schnellen Kreisel sind. Wir miissen sie mit dem 
nachst héheren Grad von Genauigkeit ansetzen und nun in dem Glied ay) auch noch die beiden 
Komponeten D, = Bod und D, = B @, sin 0 von ® in der positiven Knoten- und Querachse hin- 
zunehmen, welche wir bisher gegen das i. a. groBe Fiihrungsglied [p D] vernachlassigen durften. 
Dabei ist B die aquatoriale Drehmasse des Kreiselkérpers um jede Achse durch den Kardan- 
mittelpunkt O senkrecht zur Figurenachse. In dem Glied [p D] dagegen darf nach wie vor D, allein 
die Rolle vonD tibernehmen. Diese Erweiterung ist deswegen nétig, weil es, wie wir sehen werden, 
Sperrungen gibt, bei denen die Komponenten von [0 ®] in der Knoten- oder Querachse ver- 
schwinden, so daf dann die dortigen Komponenten von ® nicht mehr gegen sie vernachlassigt 
werden diirfen. Man erhalt so statt (7,2) und (7,3), wenn man auch hier wieder die Komponenten- 
gleichung in der Querachse nachtraglich mit sin 6 erweitert, 


Bd + D.wp)sind =Qsind = M, (620), (8.1) 
— Bsin 6 g (w, sin 6) + D, om sind = + M, — M, cos 6 (wp z 9), (8,2) 


wogegen die Gleichung (7,1) fiir D, unverandert gilt, so daB wir in (8,1) und (8,2) D, als bekannte 
Zeitfunktion (7,6) eingesetzt denken wollen. 


Eine erste Sperrung betrifft die freien Pendelungen des schnellen Kreisels um die Knotenachse 
bei festgehaJtener Prazessionsachse, also bei Sperrung in den auBeren Kardanlagern CC. Wir 
suchen vor allem Antwort auf die Frage nach dem erforderlichen Festklemmungsmoment in 
diesen Lagern. Die Gleichungen (8,1) und (8,2) lauten fiir diesen Fall mit m, = 0 


Bd =Qsind + M, (6 = 0), (8,3) 

+ M3 = M., cos 6 + D, dsind. (8,4) 

Hier ist (8,3) die bekannte Pendelgleichung (mit trockener Reibung M,). Uhre Integration liefert 
O(t) und O(t), und dann folgt aus (8,4) als mindest erforderliche Lagerreibung der Héchstwert von 


Mj nach beiden Drehrichtungen. 
Das Gegenstiick zu dieser Sperrung stellt die Prazession mit gesperrter Drehung 6 dar. Hier 


lauten die Gleichungen (8,1) und (8,2) mit a, = 6=0 


ae Mp + Me cos 6 
Pee B sin? 6 


Mk = (Q — D, w,) sin 0. (8.6) 


Wir miissen nun wieder vier Falle unterscheiden. 


(@, z 9), (8,5) 


(49) 


26 
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a) M, > M., also 6, nach (4,12) nicht reell. Hier fiihrt (8,5) mit beiden Vorzeichen von M, 
auf einen Widerspruch. Denn mit w, > 0 wiirde jetzt die rechte Seite negativ; dies gabe negative 
Beschleunigung @, und also (da q, jetzt natiirlich mit dem Wert m, = 0 beginnt) wm, <0; mit 
py <0 aber wiirde die rechte Seite von (8,5) positiv und somit ®, > 0 und w, > 0. Es tritt 
daher bei jedem Winkel 6 Sperrung auch der Priazession w, ein, und zwar mit den Sperrmomenten 
(7,14). Dieser Fall ist schon in Ziff. 7 a) auf anderem Wege gefunden worden und gilt auch noch 
fir M, ‘= M;. 

b) M, < M., 0 < 0, reell aus (4,12). Jetat stellt (8,5) nur mit w, > 0 keinen Widerspruch 
dar. Die Prazession , mit festem Winkel 6 <6, lauft also positiv beschleunigt an und weiter. 
Fiihrt man den dann aus (8,5) folgenden Wert 
— Mp ah Me cos 0 , (8,7) 

B sin? 0 


Op = 


in (8,6) ein, so hat man das Haftreibungsmoment Mj gefunden. Natiirlich gilt dies nur bis zu 
Werten von ~,, die noch klein gegen die Eigendrehgeschwindigkeit ~, sind, weil sonst der Krei- 
sel nicht mehr als schnell angesehen werden diirfte. 

c) M, < M,, 6, <6 < 180° —6,. Man findet fiir beide Vorzeichen von wp in (8,5) einen 
Widerspruch, und dies bedeutet wieder Sperrung auch fiir die Prazession mit den Sperrmomenten 
(7,14), und zwar auch noch in den Grenzfallen 6 = 6, und 6 = 180° — 6,. 

d) M, < M., 6 > 180° —6,. Jetzt stellt (8,5) nur mit w, < 0 keinen Widerspruch dar. 
Die Prazession jw, mit festem 6 > 180° — 6, lauft also negativ beschleunigt an und weiter. Das 
Sperrmoment folgt wie im Falle b), und auch diese Uberlegung gilt nur fiir nicht zu groke Werte 
von @,. 

Wiirde man (, tiber der Abszisse u = cos 6 auftragen, so entstiinden Diagramme von der Art 
der friiheren Phasendiagramme fiir u, und zwar gemafB Abb. 4 fiir den Fall a) und gemaB Abb. 5 
fir die Falle b), c) und d) zusammengenommen. Aus ihnen kénnte man dann fiir jeden festen 
Wert u die Prazessionsbeschleunigung @, ablesen. 


9. Der Kreisel ohne Schweremoment. Wenn der Schwerpunkt S mit dem Kardanmittelpunkt O 
zusammen fallt und also nach (2,1) Q = 0 wird, so kann man fiir diesen Kreisel, der ohne Reibung 
kraftefrei genannt wird, den EinfluB der Reibung aus unseren friiheren Ergebnissen ohne weiteres 
entnehmen, wenn man dort Q = 0 setzt. Man darf dabei den Kreisel stets als rechtsdrehend an- 
sehen. 

So findet man fiir fliissige Reibung im Falle a) von Ziff. 5, also fiir ¢, = 0, d. h. den synchron 
angetriebenen Kreisel, aus (5,11) und (5,13) 


w= Nes Opa 0 (e50== 0), (9,1) 


“und das besagt, daf die Figurenachse stillsteht und die Reibung ohne Wirkung ist, ein durchaus 


einleuchtendes Ergebnis. 
Ebenso findet man fiir e, = 0 aus (5,17), (5,20), (5,22) und (5,28) 


Op = 0s ug = 0; wut; i = (e,—O)R (9,2) 


und das besagt, daB das auBere Ende der positiven Figurenachse sich dem auBeren Kardanring 

auf kiirzestem Wege ohne Umkehr und ohne Prazession nihert, es sei denn, da® sie senkrecht zu 

ihm stand, in welchem Falle sie dann stehen bleibt (natirlich labil). 
Im Falle ¢, = 0 findet man aus (5,32) und (5,33) 

OG, Uy die 


AD, 1—w3 T? 


ipod 1, O, = (e,.= 0) (9,3) 
also dic gleiche Bewegung der Figurenachse, nun aber mit einer Prazession, deren Drehvektor Op 
mit der positiven Figurenachse einen spitzen Winkel bildet, oder aber (labilen) Stillstand der 
Figurenachse senkrecht zum duBeren Kardanring. 

Im allgemeinen Falle ¢, > 0, & > 0, €) > 0 gilt Abb. 3 mit den Werten Ug = u* = 07s und 
alle Kurven ¢, = 0 sind Gerade durch den Nullpunkt des Koordinatensystems, so daB die benach- 
barten Kurven ¢, > 0 stets auch in den schraffierten Streifen bei u = 1 eindringen werden, zumal 
da sie wegen (5,38), worin ug = 0 zusetzen ist, offenbar die Geraden wu = + 1 zu Asymptoten 
haben. 

Die Ergebnisse von Ziff. 6 lassen sich ebenso leicht auf Q = 0 spezialisieren. 


ele 
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Fiir trockene Reibung ist jetzt nach (7,11) s, = co, so daB kein reeller Winkel 6p und keine 


zugehérige Kote u, existiert. Von den Phasendiagrammen kommmen also nur Abb. 4 fiir M, > M, 
und Abb. 5 fiir: M, < M, in Frage. Das Ergebnis lautet: Ist M, = M_,, dann tritt cole 
Sperrung der Figurenachse ein; ist M, < M., dann tritt véllige Sperrung der Figurenachse fiir 
alle spitzen Winkel zwischen der positiven Figurenachse und dem nachstliegenden Halbstrahl der 
Prazessionsachse ein, welche gréfer oder gleich 6, sind: fiir die kleineren Winkel nahert sich die 
Figurenachse dem duferen Kardanring auf kiirzestem Wege so, wie vorhin geschildert. Im Falle 
der Sperrung gilt die zweite Gleichung (7,14) fiir das Sperrmoment M*. 

Von den beiden Sperrungen in Ziff. 8 kommt jetzt natirlich nur die zweite in Betracht, und 
auch von ihr nur der Fall b) und der Fall d) (da die anderen Falle schon soeben aufgefiihrt worden 
sind), also eine Sperrung 6 mit dem Sperrmoment M# (8,6) (mitQ = 0) und dem ungesperrten 
Wert w, aus (8,5) (mit dem oberen Vorzeichen von M,) fir jeden spitzen Winkel zwischen der 


positiven Figurenachse und dem nachstliegenden Halbstrahl der Prazessionsachse, welcher 
kleiner als 6, ist. 


10. Der Einflu8 der Nutationen. Die bisher diskutierten Lésungen stellen Naherungen fiir den 
schnellen Kreisel dar, die einerseits auf anschaulich begriindeten Vereinfachungen beruhen und 
andererseits die beobachtharen Bewegungen insofern idealisieren, als sie die Nutationen der 
Figurenachse nicht wiedergeben. Es wird daher niitzlich sein, wenn wir zum Schlu8 nun auch 
die exakten Kreiselgleichungen aufstellen, um erstens ein genaues Bild von den Vernachliassi- 
gungen zu gewinnen, welche unsere bisher beniitzten Gleichungssysteme enthalten, und um 
zweitens — wenigstens fiir den Fall fliissiger Reibung — mit Hilfe einer Stérungsrechnung zu 
zeigen, da die bisher vernachlassigten Nutationen die friiheren Ergebnisse nicht modifizieren. 

Die exakten Gleichungen kénnen entweder mit dem Lagrangeschen Verfahren oder aber mit 
dem Drehimpulssatz in der Form (3,1) gewonnen werden. Wahlt man den letzten Weg, dann 
empfiehlt es sich, (3,1) im Koordinatensystem (e, k, q) zu zerlegen, das (Abb. 1) aus der Figuren- 
achse e, der Knotenachse k und der Querachse q gebildet wird. Beachtet man dabei, daB die 
Winkelgeschwindigkeit p durch vektorielle Addition von @, und @, erhalten wird, dann gewinnt 
man mit Hilfe von Abb. 1 die Komponentendarstellung 


D =(w, cos 0, W,, Wp sind). (10,1) 


fir 9. Da das gewahlte Koordinatensystem ein Hauptachsensystem des Kreiselkérpers ist, und 
da sich die Winkelgeschwindigkeit des Kreisels aus 0 und w, zusammensetzt, ergibt sich fiir den 
Drehimpulsvektor analog 


ee (A (mw. + @,cos 0), Bu,, Bow,sind), (10,2) 
und schlieBlich hat man nach Abb. 1 sowie (4,4) 
mM =( M0 60h, Mo 4 (10,3) 
In Komponenten aufgelést und geordnet lautet (3,1) daher 
A (We — Wx Wp Sin 0 + @®, Co8 0) =— M., | 
B oy, — Bw, cos 6 sin 6 + A (@, + Wp cos 6) mW, sind = Q sind — M,, (10,4) 
B @, sin? 6 + 2 B ax wp cos 6 sin 6 — A wy (We. + @p ©08 O) sin 6 = M, cos 6 — M,. | 


Ist der Kreisel schnell, so kann man, wie dies seine Definition nahelegt, @,/@,. und Wp/W- 
als von erster Ordnung kleine Groen behandeln und (10,4) durch Vernachlassigung von Gréfen 
héherer Ordnung auf 


A(@_ + @, cos 0) =— M., 
B oy, + A @, @p sin d = Q sin 6 — M;. (10,5) 
B wp sin? 6 — A @, Wy, sin 6 = M, cos 6 — M, | 


reduzieren. Beachtet man, dai w, = 6 und nach (10,2) D, = A (We + Wp cos 6) ist, dann stellt 
man leicht fest, daB® in der beobachteten Naherung und bei Beschrankung auf trockene Reibung 
die drei Gleichungen (10,5) der Reihe nach mit (7,1), (8,1) und (8,2) identisch sind, womit nun 
nachtraglich die Ergebnisse von Ziff. 8 von der analytischen Seite her als Naherungen erster 
Ordnung gerechtfertigt werden. Das Gleichungssystem (4,1) bis (4,4) stellt dagegen eine grobere 
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Naherung dar und folgt aus (10,5) mit der weiteren Annahme, daf die GréBen Cx) e und @p/M¢ 
mindestens von zweiter Ordnung klein seien, womit dann auch (wie zu Beginn von Ziff. 5) 
D, = Aw, gesetzt werden kann. Diese Zusatzannahme mag etwa mit der Feststellung gerecht- 
fertigt werden, da8 bei annahernd konstantem @, (d. h. etwa im Falle M, = 0) a/@¢ und Wp/W-¢ 
nur dann iiber gréBere (d. h. von erster Ordnung grofe) Zeitintervalle klein von erster Ordnung 
bleiben kénnen, wenn @/w? und @,/w: klein von héherer Ordnung sind. Mit Riicksicht auf den 
Vorbehalt wm, ~ konst. sowie auf die Méglichkeit, daB c,/a@~ und @p/@, oszillierenden Charakter 
haben, kann aber dieses Argument keine strenge Giiltigkeit beanspruchen, und in der Tat werden 
ja mit dem System (4,1) bis (4,4) gerade die raschen Oszillationen in Form der Nutationen nicht 
erhalten. 

Wir gehen deshalb von der strengeren Naherung (10,5) aus, wobei wir uns aber vorbehalten, 
die eben diskutierte Zusatzannahme allenfalls spiter wieder einzufiihren. Das System (10,5) 
besitat, falls keine Lagerreibung vorhanden, mithin M, = M, = M, = 0 ist, ein partikulares 
Integral in Form der regularen Prazession mit den konstanten Bestimmungsstiicken 


We = Deg, ) = Oo 5 Wp = oe = Wp og e (10,6) 
Die zugehérige Bewegung tritt freilich nur unter speziellen Anfangsbedingungen und nach der 
dritten Beziehung (10,6) nur unter der Voraussetzung auf, daB Q/A wéo von erster Ordnung klein 


sei. 

La®Bt man jetzt kleine Abweichungen von den erwadhnten speziellen Anfangsbedingungen 
und gleichzeitig Reibungsmomente von der GréSenordnung des Stiitzmomentes Q zu, so dah 
auch M,/A wec, M;/A wo und M,/A wo von erster Ordnung kleine GréBen sind, dann sind aus- 
schlieBlich Lésungen zu erwarten, die sich in einem hinreichend beschrankten Zeitintervall nur 
wenig von (10,6) unterscheiden. Man kann sie in der Form 


We = Wey + Wey» 6=6,+6,; Op = Mp9 + Opt (10,7) 
ansetzen, wobei die GroBen @,1/@eo, 6, und @p;/@.9 als klein von erster Ordnung angenommen 
werden diirfen. In dieser Naherung gilt 

sin 0 = sin dy + 0, cos by, cos 0 = cos 6) — 0; Sin by; (10,8) 


und wenn man jetzt (10,7) in (10,5) einsetzt, (10,6) sowie (10,8) beriicksichtigt und alle von 
hoherer als erster Ordnung kleinen GréBen unterdriickt, erhalt man fiir @,,, @,, = by und Wp, 
die linearisierten Differentialgleichungen 


A (Wey + Wp1 cos 0o) ae ae M. ; 
B @p, sin? dy — A Weg Wp Sin 6g = M, cos 6) — M,. 
Wir gehen jetzt von den Winkelsgeschwindigkeiten w, = Wey, W, = 0, (Wp = 0 aus, die dem mit 
Q = 0 permanent rotierenden Kreisel entsprechen, und setzen die Reibungsmomente als Funktionen 


der Winkelgeschwindigkeiten voraus, die sich in der Umgebung dieser Werte in so rasch kon- 
vergierende Potenzreihen 


M.(a,) = Mey rye rs [Meo = M.(@e) > 0] ’ 
M,,(@x«) = Mio OP AR 25 [ Mi, 0 — M;,(0) 205 (10,10) 
M,(0p) = Mio (po + px) °° [Mig = Mi(0) > 0] 


entwickeln lassen (wobei Striche Ableitungen bedeuten), da8 im Rahmen unserer ersten Naherung 
nur die (angeschriebenen) ersten Glieder zu beriicksichtigen sind. Mit dieser Annahme scheidet 
die trockene Reibung aus, dafiir bleibt aber die geschwindigkeitsproportionale fliissige Reibung 
gemaB Ziff. 5 samt Verallgemeinerungen zugelassen, welche diejenigen von Ziff. 6 in gewissem 
Sinn erganzen. 
Setzt man (10,10) unter Beschrankung auf Gréf®en, die klein von erster Ordnung sind, in 
(10,9) ein, dann kommt 
A (@p; C08 6g + @) + My =0, 
B On, + Mig Oxy + A eq Wp; sin dy = 0, (10,11) 
B cp, Sin? 6g — A Wey Wp Sin by + Mig (@po + @p;) — M.o cos 6) = 0. 
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Aus der zweiten Beziehung (10,11) folgt jetzt 


. / 
B Oy, + Myo Ox, 
A @,9 sin 09 


Op, =— : (10,12) 


und wenn man damit in die dritte eingeht, erhalt man fiir w,, die Differentialgleichung 


Oey +2Y Wey +22 we, +A =0 (10,13) 
mit den Beiwerten 
ih , Le nie ese 
Y =a Bantd, Meo + Meo sin’do) , | 
5 1 . o / / 
= B? sin? 6, (42 iro sin? Oo ++ Mio M, 0) ’ (10,14) 
A We , 
es Bsn O (Me 608 0g — Moo po) » | 
von denen die ersten beiden auf Grund der Nebenbedingungen in (10,10) positiv sind. 
Die allgemeinste Lésung von (10,13) lautet 
Op, =— J aera, (10,15) 


wobei die Punkte eine gedampfte Schwingung andeuten, die offensichtlich den Beitrag der 
Nutation zur Nickbewegung verkérpert, wahrend der ausgeschriebene Term je nach dem Vor- 
zeichen von A ein monotones Steigen bzw. Sinken des Schwerpunktes darstellt. Setzt man (10,15) 
in (10,12) ein, so kommt ferner 

M;o i 


Soa ER EE 
A Weq sin 69 % 


on , (10,16) 


wobei auch hier die Punkte fiir den Nutationsanteil in Form einer gedimpften Schwingung 
stehen. SchlieBlich liefert die erste Beziehung (10,11) analog 


Meo eo 
abit (10,17) 


Ve. =— 


wobei die Punkte nun den Einflu8 der Nutationen auf die Eigendrehgeschwindigkeit andeuten 
sollen. 

Nachdem sich damit herausgestellt hat, daB die durch (10,15) bis (10,17) dargestellte Losung 
erster Ordnung die Nutationen enthalt, beschranken wir uns auf die Diskussion der (ausgeschrie- 
benen) sakularen Bestandteile, die ja offensichtlich von den Nutationen unabhangig sind. Nach 
(10,10) ist der zweite Klammerterm in der dritten Beziehung (10,14) mit dem in Ziff. 4 defi- 


nierten Moment M, identisch, und da man nach (10,10) und (10,8) in unserer Naherung fiir den 
ersten Klammerterm auch M, cos 6 schreiben kann, ist 


AZ=0 fir M,cosd—M,=0. (10,18) 
Damit folgt aber aus (10,15) und daraus, daB nach (10,6) wz, = 0, also wm, = wp, ist, 
eat) pte pe te (10,19) 
und analog aus (10,16) 
p50 fir  M,cosdS M,. (10,20) 


In (10,19) und (10,20) stecken alle Aussagen von Ziff. 4 iiber die Nickbewegung und die Pra- 
zession. Diese allgemeinen Resultate bleiben demnach in vollem Umfang giltig, auch wenn 
man auf die Zusatzannahme, da @,/w29 und @,/@e¢o von héherer Ordnung klein seien, ver- 
zichtet und auf der strengen ersten Naherung (10,5) basiert. 

Im Gegensatz hierzu ist festzustellen, da® die sikularen Terme in (10,15) bis (10,17) keine 
Lésung der Gleichungen (4,1) bis (4,4) darstellen, sondern lediglich gegen eine solche konver- 
gieren, wahrend die Reibungsmomente gegen Null gehen. Es folgt hieraus, daB die quantitativen 
Ergebnisse von Ziff. 5 bis 7 zwar durch die Nutationen nicht modifiziert werden, aber trotzdem 
auch innerhalb einer Theorie erster Ordnung nur Naherungen darstellen, die freilich bei kleinen 
Reibungsmomenten hinreichend genau sein diirften. 
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11. Versuche. Um die Ergebnisse dieser Arbeit experimentell zu iiberpriifen, wurde ein von der 
Firma Schuler & Séhne in Ziirich gebauter, durchweg mit Kugellagern ausgeriisteter Kardan- 
kreisel ohne Eigenantrieb mit Bremsen versehen, welche die Reibung in den Lagern AA, BB 
und CC zu vergréBern und damit gegeneinander abzustimmen gestatten. Zur Zeit liegen erst 
die Ergebnisse einiger qualitativer Vorversuche vor, die indessen die allgemeinen theoretischen 
Resultate von Ziff. 4 ibereinstimmend bestatigen. 

So zeigt z. B. der Versuch, da der Schwerpunkt des Kreisels ohne Schweremoment (Q = 0) 
unter einer positiven Prazession bis in die héchste Lage steigt und unter einer negativen Pra- 
zession bis in die tiefste Lage sinkt, je nachdem er urspriinglich gehoben oder gesenkt ist. Das 
Winkelintervall um 6 = 90°, in dem allenfalls Sperrung eintritt, ist ohne kiinstliche Bremsung 
klein. Wird die Reibung M, in den inneren Kardanlagern AA verstarkt, dann verlaufen beide 
Bewegungen, wie auf Grund von (4,2) bis (4,4) zu erwarten ist, rascher. 

Wird cin Schweremoment Q angebracht, dann verschiebt sich die Grenze zwischen Steigen 
und Fallen des Schwerpunktes in Ubereinstimmung mit (4,12) zu spitzen Winkeln 6,; fiir einen 
bestimmten Wert von Q liegt sie beim ungebremsten Kreisel beispielsweise etwa bei 70°. Auch 
jetzt werden sowohl die Prazession wie die Nickbewegung bei Verstarkung der Reibung in den 
Lagern AA rascher, wahrend gleichzeitig 6, im Einklang mit (4,12) gegen 90° geht. 


(Eingegangen am 2. Februar 19506.) 


Anschriften der Verfasser: Professor Dr. Richard Grammel, Stuttgart-N, Robert-Bosch-Str. 101, 
und Professor Dr. Hans Ziegler, Riischlikon bei Zurich, Weiherweg 6. 
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Instationire Warmespannungen in einer Hohlkugel 
Von R. Trostel 


1. Einleitung. In der vorliegenden Arbeit wird der Spannungszustand infolge instationdrer 
rotationssymmetrischer Temperaturfelder in einer Hohlkugel untersucht, deren innere und 
auBere Oberflachen unter dem EinfluB zeitlich und bis auf die Achsensymmetrie auch 6rtlich 
beliebiger thermischer Randwirkungen stehen. 

Die Ermittlung eines Warmespannungszustandes bedingt zunachst die Bestimmung des 
zugehorigen Temperaturfeldes, und daher wurde in Abschnitt 2 zundchst einiges tiber das 
achsensymmetrische Temperaturfeld vorangestellt. Hieran schlieBt sich in Abschnitt 3 die Er- 
mittlung der Warmespannungen in einer dicken Hohlkugel an. Das Ziel war die Ermittlung eines 
oberflachenspannungsfreien Warmespannungszustandes, der schlieBlich durch Uberlagerung der 
aus dem sog. thermisch-elastischen Verschiebungspotential! hervorgehenden Warmespannungen 
mit gewissen, durch Lovesche Funktionen charakterisierten Randspannungszustanden gefunden 
werden konnte. In Abschnitt 4 werden die Lésung fiir den kugelsymmetrischen Fall und einige 
Naherungsformeln angegeben, wahrend Abschnitt 5 einige kurze Hinweise auf das achsensymme- 
trische Problem an diinnwandigen Kugeln (Kugelschale) enthalt. 


2. Das achsensymmetrische Temperaturfeld infolge zeitlich und értlich beliebig vorgegebener 
thermischer Randwirkungen an den Oberflichen. a) Grundlagen der Wirmeleitung; die 
Warmeleitungsgleichung. Fir den rotationssymmetrischen Fall (00/dp = 0) lautet die 
Differentialgleichung der Warmeleitung in Kugelkoordinaten 
(Abb. 1) 

eo, 200 1 ( ao ad 


Avg. = - DT CV200 an 00 
AD = aa 4 aot 


eto a) — | ef aa 
ore Y og Gye) Ree or (1) 


fe Tr 


wobei c [cal/gr °C] die spezifische Warmemenge je Gewichtsein- 
heit, y [gr/em*] das spezifische Gewicht und A [cal/cmsec °C] die 
Warmeleitzah] hedeuten. Gleichung (1) driickt das Gesetz der 
Erhaltung der Warmeenergie aus und zwar mit Verwendung des 
Fourierschen Warmeleitungstheorems. MHiernach ist die pro 
Zeiteinheit durch die Flachencinheit hindurchtretende Warme- 
menge q (cal/cm? sec) proportional dem Temperaturgradienten 
a0. 1 a8 
q = {43 dy; 0} = —A grad 0 = 14 Aoi ay? Of (2) 


Fir die Formulierung der Randbedingungen hinsichtlich eines Warmetberganges in ein anderes 
Medium benétigt man noch eine zusatzliche Beziehung, die uns das Newtonsche Abkihlungs- 
theorem liefert. Entsprechend dem Fourierschen Theorem wird hier die aus dem erwarmten 
Kérper in ein umgebendes Medium der Temperatur # y tibergehende Warmemenge je Flachen- und 
Zeiteinheit ; 

Rand = QRand Vl = — & (Om —- OR) . (3) 
also proportional zur Differenz der Temperaturen des Mediums (%y) und der Kérperoberflache 
(Or) gesetzt. Hierbei ist n der Vektor der aufSferen Normalen der Kérperoberflache und a [cal/ 
em?sec °C] die Warmeiibergangszahl. 


b) Die Rand- und Anfangsbedingungen. «) Randbedingungen. Die in der aufzufinden- 
den Lésung enthaltenen Randwirkungen an den Kugeloberflachen r =r; bzw. r =r, sollen 


mathematisch formuliert sein konnen: 
1. Durch zeitlich und 6rtlich beliebig veranderliche vorgegebene Randtemperaturen @;(y, t) 


bzw. d.(y, t) der Kugeloberflachen, 
2. durch zeitlich und 6rtlich beliebig veranderliche vorgegebene Temperaturen Oui, t) baw. 


Oma(yp, t) der im Innern bzw. am Au®enrande befindlichen Medien, 


1 FE. Melan und H. Parkus, Warmespannungen stationirer Temperaturfelder. Wien 1953. 
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3. durch zeitlich und Grtlich beliebig veranderliche vorgegebene, durch die Oberflachen der 
Hohlkugel hindurchtretende Warmemengen q;(y,t) baw. qa(p, t) sowie 
4. durch Kombinationen zwischen ]., 2. und 3. 


Im Falle 1. lauten die Randbedingungen 


Or; Ws t) = OY, t) ’ O(Tas Ws t) = Val, t) (4 a) 

und im Falle 3. mit (2) 
ce ees Bu a al res Ab 
a ‘Or a wy, t ie gil, 4 ; A or Ti, Ye TAY ‘ ( ) 


Zur Darstellung der Randbedingungen im Falle 2. kann man zunachst nach (3) fiir die Oberflache 
r=7;, indem man Qrand = Q(ris Pt) = {ar(Tir Y 13 Q(T Y> 3 Of} und n= {ips Rye aly, 
= {—1; 0; 0} setzt, Grong 1 = —9(r:. y, t) =— % [Imidyp, t) — O(r;, y, t)] und entsprechend 
fiir die Oberflache r = r, schlieBlich q,(ra, yp, t) = — « [Oma(y, t) — O(a, y, t)] folgern, und 
damit erhalt man mit Beachtung von (2) 


(ri, p, t) = & [Omi(y, t) — Ori Y, )] 


OP | ae Od | 
Se) Oe - be QATas Ys t) =— 4 [Omalp, t) — Ara, yp, t)] = — A ze ae 5 
lh A ad| 
A 08 | done te ho emaneeey oe 
Ora, W, t) iS O or ir = O mals t) ’ (Tis W, t) Gl Gr bebe a Omi, t) . (4c) 


Wir fassen nun sdmtliche durch (4a) bis (4c) beschriebenen Randwertprobleme und damit 
auch den Fall 4. durch folgende allgemeinere Randwertaufgabe, in der jene als Spezialfalle ent- 
halten sind: 


a8 
= alps t)s by Arn yst) + by 


Ta> Yt 


oo 
a, Hra, py, t) + dy a 


=e). 5) 
THY, t 
Hierbei sind a,, a,, 6; und 6, zunachst willkirliche Konstanten, und indem man beispielsweise 
a, = b, =1, ag =b, = 0, aly, t) = Daly, t) und —i(y, t) = d,(y, t) setzt, folgt sofort der Fall 1. 
usw. 

Mit den Randbedingungen (5), von denen im folgenden ausgegangen wird, sind demnach 
beliebige thermische Randwirkungen an den Kugeloberflachen erfafBt. 


B) Die Anfangsbedingung. Die Lésung soll einen zur Zeit t = 0 gegeniiber seiner Umgebung 
nicht erwarmten Kérper zum Gegenstand haben. Bezeichnet man die anfangliche, iiberall 
als gleich vorausgesetzte Temperatur der Umgebung mit Null, so lautet also die Anfangs- 
bedingung 

Hr, y, 0) = 0. (6) 

c) Allgemeine Bemerkungen zur Ermittlung des Temperaturfeldes. Eine un- 
mittelbare Entwicklung der den Rand- und Anfangsbedingungen (5) und (6) angepaBten Lésung 
fiir das Temperaturfeld ist nicht méglich. Hierfiir mu vorbereitend die Aufgabe gelést werden, 
das Temperaturfeld fiir den Fail zeitlich von t = 0 an konstanter Randwirkungen an den Kugel- 
oberflachen zu bestimmen. Die Anfangs- und Randbedingungen lauten in diesem Falle in 
Spezialisierung von (5) 

oo 


ob % 
a, W(Tar Ys t) + ay 3 =faly), b; O(ris p, t) + be ae =fdy) firt>0, (7) 


Tar yt Ts, Ws t 

Hr, yp, 0) = 0. (8) 
Diese Aufgabe wird dadurch bewaltigt, daB man zunachst das Randwertproblem (7) durch die 
sog. stationére Lésung %,(r, y) befriedigt, und zur stationdren Lésung dann noch eine Zusatz- 
lésung 3,(r, y, t) instationdren Charakters hinzufiigt, die die bereits durch die stationdre Lésung 
erfiillten Randbedingungen nicht mehr verandert, jedoch nun auch noch zusammen mit (Tr, y) 
der Anfangsbedingung (8) Geniige leistet. In der nachfolgenden Ziff. d) erfolgt die Ermittlung 
der stationdaren Lésung #(r, y), wahrend die Ziff. e der instationaren Zusatzlosung @, (r, y, t) 
gewidmet ist. In Ziff. f wird dann durch Uberlagerung von #,(r, y) mit 0,(r, y, t) die Lésung 
fiir die Rand- und Anfangsbedingungen (7) und (8) gewonnen, und schlieBlich gelingt in Ziff. g) die 
allgemeine Lésung fiir zeitlich und értlich beliebige Randwirkungen an den Kugeloberflichen nach 
Mafgabe der Rand- und Anfangshedingungen (5) und (6) 


. 
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d) Das stationire Temperaturfeld. «) Allgemeines: Produktlésung der Wéarmeleitungs- 
gleichung. Im stationaren Falle folgt aus (1) wegen 0%,/dt = 0 die Differentialgleichung 
Pee a ee 20, , a,\ 
Ad, Or? e Or | re (re ? Oy + ay) aby 


Fur ihre in Produktform darstellbaren Lésungen folgen nach Einfihrung des Produktansatzes 
Y(t, p) = R,(r) Yo(p) gemaB der Daniel-Bernoullischen SchluBweise mit der Konstanten n (n +1) 


tei) ln aw = die beiden gewohnlichen Differentialgleichungen 
d?Ron dRon 
Pp a ee Fe n(n + 1) Ro»(r) = 0, (9a) 
GEM d¥, 
“Ty + ctg yp A + n(n + 1)¥%, =0. (9b) 


Gleichung (9a) ist eine Kulersche Differentialgleichung, fiir deren allgemeine Lésung man schlieb- 
lich mit dem Ansatz R,,(r) = A, r*. 


Ro, (7) = Ain rt + Bm (10) 
erhalt, wahrend Gleichung (9b) 
1 d 1 dP, oe 
any ay (sin y a +n(n + 1) ¥%,(p) =0 


mittels der Transformation 


: os dé : 
Pony) = PonrlE(p)] mit é(y) = cosy, d-h. 1—& = sin? y baw. aes siny (lla) 
in die Legendresche Differentialgleichung 
d ave 
mall _s) | Boy (nel en (eye (11b) 
mit der allgemeinen Lésung 
Y% mes) =— on (cos yp) = Ay n PAé) ali Bon Q,.(¢) (11c) 
uberfihrt werden kann. P,(&) bzw. Q,(&) sind hierbei Kugelfunktionen erster bzw. zweiter Art. 
Nun sind wegen lé| = |cos y| <1 die Kugelfunktionen zweiter Art im vorliegenden Falle nicht 
definiert, so daB anstelle von (llc) als Lésung von (11b) 
Po n(6) = Won (cos y) a ots PE) (11d) 


verbleibt. Mit A;, 42, — A, und B,, A2, = B, hat man dann die Produktlésung 


‘BA 
Ron(t) Lon(y) = (4. mi ; Pee 
die wir nun noch in der Form 
. — B, 
Or. cos y) - I(r, é) == ody Ro,(r) Yon(y) — » (4. Te ln ” | PAG) (12) 
n=0 —) 


verallgemeinern kénnen. Die Konstanten A, und B, (n= 1, 2,3...) sind nunmehr den Rand- 
bedingungen des Problems anzupassen. 


B) Lésung der stationéren Randwertaufgabe. Von der unter Ziff. da) gefundenen Produkt- 
lésung (12) soll nachfolgend entsprechend den Ausfiihrungen der Ziff. c) die Erfiillung folgender 
Randbedingungen gefordert werden 


i 


OF, : q dv, | Q 
ay Oo(Tas Y) + Tear , z =faly), 5 Dolris y) + be ar | = fily) » (13) 


TE 


und da die Randfunktionen f, und f; ohne Schwierigkeit in der Form fa = fu(€) = fa (cos py) baw. 

fi =f© =f. (cos y) darstellbar sind, so erhalten wir, indem wir (12) in die Randbedingungen 

(13) einftihren, mit | 
Aa pir. ) By (ary OF (n+ 1) a, ry Aged) ie angen | 


(14) 
A,, (by 1? + nb, re —1) + By (b, 7 +N _ (n+ 1) bret?) = oye f 
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folgende Beziehungen: 


So Pe) =f), So Prlé) = filé)- (15) 
0 0 


eS n == t 


Mit Beachtung der Orthogonalitatsrelationen! 


foal Q) tile Wise De 
( £ pat 16 
J Palé) Pu) dé |s eave (abi Tp) == 1) (16) 
gewinnt man nunmehr aus (15) folgende Beziehungen fiir die ie baw. c?: 
— J 
p A Pip e | 
c= SRE | fle) Pale) dB oP = SS a fale) PalE) dé (17) 
Mi 


womit diese und damit auch nach (14) die 4, und B, als bekannt anzusehen sind. Man erhalt 
schlieBlich mit 


AO) = (ar + may 72 —}) (by r+) — (nm 4 1) by rp +) 


| 
=(o —@+1)—_(n + I) agr-@+ *)) (b, ry + nb, re 
ow = ~ (b, ra (n +1) __ (n + 1) by rae + 2)) 
i - (18) 
xP — eC r— +1) (n + 1) a, eee *)) : 
: I 
Os saath rm + nb, i ae m= Gow (47 TS etd, ra) i, 


fur die Konstanten A; und 155, 
An Ne U S ) a. C a B B 2) Bye Bee (a | 9 
(i) c i) © “ ) 2G Bee ( ) 


Die Lésung der stationéaren Randwertaufgabe ist damit vollstandig bestimmt 


Ello 6b te oe) +” (BP oP + fe oy r+) AC) 

s (20) 
=— > yr(r) P(E) r) | ) 
wenn man noch zur Abkiirzung (4? ch? + of? cM?) r™ + (BY? cP + BS ci) r—~@ +) = ¥,(r) setat. 


e) Die instationare Zusatzliésung 0,(r, yp, t). 0) Allgemeines; Crane der Warme- 
leitungsgleichung. Die maBgebende Gleichung fiir das instationére Temperaturfeld @, (r, y, t) 
ist die Differentialgleichung (1), und wir gewinnen hieraus fiir die in Produktform darstellbaren 
Lésungen mit dem, die zunachst willkiirliche Konstante @? enthaltenden Ansatz 


SAAD) na (") “R(r) Py) (21) 
1 /d2R 2 dR eal dY d2 Vv P 
R G “4 r ar) ! if (ete y dy a | ee@es 
woraus mit der Bernoullischen SchluBweise 


a7Ry 2 bel i a dee 
R,, a APE = = ae Ry) See (cts a dy. ae dip? =n(n-+ 1) 


r 


die Beziehung 


d. h. die beiden gewohnlichen Differentialgleichungen 
2 
d?R,, 2..dK,, ! (* __ n(in+)) 


r= 


JR = 0 ’ (22a) 


d?¥ aha 
= pete a +n(n+1)¥, =0 (22b) 


' Siehe z. B.: R. Rothe und I. Szabé6, Hohere Mathematik, Bd. VI, § 8.4. Teubner 1953. 
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hervorgehen. Wahrend Gleichung (22b) mit (9b) identisch ist, also eine Legendresche Differen- 
tialgleichung mit der Lésung 


Vicon)? (e) —=—A4,,P(é) (23) 


darstellt, laBt sich die Differentialgleichung (22a) mittels der Transformation von Lommel? 
auf eine Besselsche Differentialgleichung reduzieren. Setzt man nimlich 


Rul) = fl) (24) 


so erhalt man hiermit aus (22a) die Besselsche Differentialgleichung 


el 
eg pe eer, r) = 0, (25) 


or (w r)? 


wobei die Striche Ableitungen nach dem Gesamtargument wr bedeuten, mit der allgemeinen 
Lésung 


7 
< “ 


filo r) == (iy J ast (w r) + Can ve is 3) (@ r) : (26) 


Mit (26) steht damit nach (24) die Lésung von (22a) fest, und man erhalt das Fundamental- 
system 


Lo = = ie " J ne (or) + Cond _ (> a (w ” ; (27) 


so dah nunmehr gemaf (21) fiir die in Produktform darstellbaren Lésungen von (1) mit A, Ci, 
= A, und A, C2, = B, verallgemeinernd 


Gere 0 ees [AJ o)+ BI ylon|Pie ee 


iz - 
EP ome ee 


hervorgeht. 


fp) Das Randwertproblem des Temperaturfeldes J), (r, y, t). GemaB den Ausfiihrungen der Ziff. c) 
wird von der durch (28) gegebenen instationaéren Zusatzlésung verlangt, daB sie bei der nach- 
foigenden Uberlagerung mit der stationdren Lésung @,(r,y) die bereits durch die stationare 
Lésung erfillten Randbedingungen (13) nicht mehr verdndert. In diesem Sinne hat man also 
von #,(r, y, t) die Erfiillung folgender Randbedingungen zu fordern: 


ea 


2 | 
or ian wy, ¢ 


a, O,(ra, yp, t) +- ay ee) =0, A(T yp, t) + 6 oe (29) 


Diesen kann mit der Lésung (28) offenbar dann geniigt werden, wenn man fiir jedes cinzelne 


Glied (n) der Reihe 
a R,(Ta) + dy a alias 0, b, R,(ri) ar by ae =! (30) 


| a 


dr 


oder, ausgeschrieben, das in A, und B,, lineare homogene Gleichungssystem 


A, le J. il{@r,) + o (2 Ci si (QT 5) = J | 1 (@ ra) 


Ein alse Yat ry 


a) (wm ra) Se ke / i (w r))| = 0, 


ati By Gg J (« i =) (@ Ta) ia 52 (2 W Tad (" af — (n + 
> a st ray ee 


9 
oy, 


(31) 


AY 2, all i (@r;) + 2 (2 ord a 


2 


wie Ba |b oe c 1) (@ ri) + Pe (2 Orn J. y (Or)imas (x if ne rd) =i 
eae an. on 


iF Siehe z. B.: R. Rothe und I. Szabé6, Héhere Mathematik, Bd. VI, § 15.10. Teubner 1953. 
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folgern kann. Damit (31) fiir von Null verschiedene A, und B, (n > 0,4; 2-2) erfulltise 
ist es notwendig, das die Nennerdeterminante des Systems (31) verschwindet 


A,(w) = la ¥f a (@ ra) + af (2 Orad 1(@ra)—JS na (@ rs) 
LD id = 


27, Tati a 


x pce) ee o rl ae a 1) (7) 


2 


Nears w+ — ; 


fet a Lb ( : ss (@ ra) + ie (? Ta ye S (w ra) — it (n . =) (w ra) 


27; n+ 


x G ip Auten (: Bi Je (ee eee ") re (32) 
n+ n+ a 


2 
2 


womit eine transzendente Gleichung fiir die Eigenwerte w gewonnen ist, deren Lésungen kn 
(k = 0, 1,2...) fiir das folgende als bekannt vorausgesetzt werden sollen. Aus dem Gleichungs- 
system (31) folgt iiberdies eine Abhaingigkeit der Konstanten A, und B, in der Form 


ie i na (ort 72 (20 nS, Pe (or—J_, 1%) 
UL erSery a n PAE 2 


B, =— Ana 
a J ( i a ( ra) + zr,(? ory J” (1 : 4) (@ ra) —J _ E 5 1) (@ 2) 


2 / 


so da man mit der neuen Konstanten 


A, 
(pp nee aml — : 


a, ee ) a é aye =e OTe a he eet. ye) (w rd) 


Pek 
SD a 


somit 


R,(r) == poe AGT) a B. ae ae | 


w Dy 


= a an UP) aye? Cred (ae ogee etic ft 8 rs) 


v T 2 5) 


as ha yal : Or) gs, (2 Oe 1 (@ ra) 
fae ee : pb 


Orr. 
2 2 


ee @ rd)| 4 . 5 (@ ”) = €,O,(o7T) (33) 


n + — 
2 


und damit fiir die Lésung (28) 


Ne 

O(r, as 1) sane (") Ce Pr (on) oP AE) (34) 
n=0 

erhalt, und da (34) fiir alle nach (32) berechenbaren Eigenwerte «,, (k =0,1,2.. -) die ge- 

stellten Randbedingungen (29) erfiillt, ist die Lésung (34) in der folgenden Form verallgemei- 

nerungsfahig : ‘ 


Okn\* 
oD oO play ees Fe 
O,(r5p, ) = 3) DY) Cun @p (Onn t) Palé)-e ( a : (35) 
k=0 n=0 
Die noch zunachst unbestimmten Konstanten C,, sind der bisher nicht beriicksichtigten An- 
fangsbedingung des Problems vorbehalten. 


y) Das Anfangswertproblem des Temperaturfeldes §,(r, y, t). Wie bereits in c) erwahnt, 
erfiillt die stationdére Lésung (20) zwar die Randbedingungen (7), jedoch nicht die Anfangs- 
bedingung (8), und es muf daher von der instationaren Zusatzlosung ?,(r, yp, t) verlangt werden, 
daB sie der Anfangsbedingung 


I(r, p, 0) = — Br, y) (36) 
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geniigt, d.h. zur Zeit t = 0 gleich der negativen stationdéren Lésung wird, so daB sich die Wir- 
kungen von %,(r, y) und 9,(r, y, 1) zur Zeit t = 0 in der Tat aufheben. Indem man in (36) die 
Beziehungen (20) und (35) einsetzt, erhalt man zunichst 


foe) 


DES Gata hee eS extn) B12). 


n==0 [k=0 n=0 


woraus fiir jedes Reihenglied (n) 
>» Cen ®,, (Mkn r) = Yn(T) (37) 
=0 


gefolgert werden kann, und beachtet man nun die Orthogonalitatsrelationen der EKigenfunk- 
tionen ®, (@zn 1) 


r 
a 


‘if PD, (Mxnr) Dy (@in r) dr 


| O4ir eth, 
in ron, (ein -———=11) (nu is 1)) DP, + Went ®,, D,, = Okn iP ®,}| 


CG ie eee | (38)! 
so erhalt man schlieBlich fiir die in (37) verbliebenen Konstanten C;,, indem man beide Seiten 
von (37) zunachst mit r? ®,(w,, 1r) multipliziert und hernach in den Grenzen r; <r <r, inte- 
griert 


Cin Sa ae | Ta Hr(7) ®,, (Win r) dr (39) 


und damit letztlich fiir die instationaére Zusatzlésung aus (35) mit (39) 


& 


Bi) tes eS ee 
v, (r, wv, t) = — > > ie ®,, (@knt) P,(&) e ( é | r? y,(r) ®,(@kn 1) dr. (40) 


k=0 n=0U S 
U 


f) Das Temperaturfeld fiir die Anfangs- und Randbedingungen (7) und (8). 
Dieses wird gemaB der vorangegangenen Uberlegungen durch Uberlagerung der stationdren 
Lésung @,(r, y) mit der instationaren Zusatzlésung 2,(r,y,t) gewonnen. Dementsprechend erhiilt 
man mit (20) und (40) 


Or, y,t) = I(r, y) + A(r, y, t) 


FS - et hn . vs = 
= > hl — > Be e | i ) D,, (Wk nT) | rn,(r) ®,(@xzn 1) dr| P,(é) . (41) 
A ae oar os : 


g) Die allgemeine Lésung fiir das achsensymmetrische Temperaturfeld nach 
MaBgabe der Rand- und Anfangsbedingungen (5) und (6). Fir die Bestimmung dieses 
Temperaturfeldes wird von der Lésung (41) ausgegangen, die man mit 


W(r) Shere py poet) und a) Sgr? ak BO ra er) 
in der Form 
e 1 —(*E*) ° 
O(r, Y, t) — >; Ww (r) —_— om é a PD, (Mkn r) | r? nn (r) D,, (Mkn r) dr i 1 (é) 
n=) L k=0 ie a 


L 


tt 
a 


7 ene 
e > aa = a € ( 3 ) P,, (Mkn 7) | 72 yr) Dy (Win T) ir 6o4 PAE) (42) 


n=0 L k=0 = ; 
L 


1 Die Striche deuten wieder Ableitungen nach dem Gesamtargument ,,7r an. 
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schreiben kann. Sie charakterisiert ein Temperaturfeld zur Zeit t infolge der zur Zeit t=0 
einsetzenden und zeitlich konstanten mittels der Funktionen f,(€) baw. f(é) gekennzeichneten 
Randwirkungen. Im Falle zeitlich verinderlicher Randwirkungen @a(&,t) baw. gi(€,t) sind nun 
die Fourierkoeffizienten c® bzw. c? nicht mehr konstant, sondern Funktionen der Zeit, was 
sofort aus den Entwicklungen von —,(é,t) baw. y;(&,t) nach den Eigenfunktionen P,, (é) hervor- 


geht — man erhalt aus 


Pa (é, t) ee SS cn Ae) > gil, t) ee OS eu PAs) (43a) 
n—O n—0 
wegen (16) 
4-1 4-1 
cp =P OES | galEst) Palé) dE = 000), o6? == 2 I | gE, 2) PylE) dE = ex)» (4b) 
ail rey 


also zeitabhangige K oeffizienten — und es entstehen nun wahrend eines belie bigen Zeitelementes dt 
die Zusatzwirkungen 


of) Op Ovi 
0: 9(6, 0) = ae dt bzw. 0:9:(& 7) = a dt , 
(a)_ don? 
den = Gg At d.h. die Zuwachse 


dc dc 
de® =—" dr bzw. de? = dr. 
: dt Uh dt 


Betrachten wir nun vereinfachend zunachst allein 
Anderungen am auBeren Rande, und nehmen wir 
an, da eine zur Zeit t entstandene Wirkung dc{(r) 
von nun an zeitlich konstant weiterwirken wiirde 
(s. die obere Figur der Abb. 2), so wiirde hierdurch 
das zur Zeit t (nachdem also t—t Sekunden nach 
Einsetzen der Randwirkung dc\” vergangen sind) 


vorhandene Temperaturfeld mittels der Lésung 
(42) dadurch erfaBt werden kénnen, da man in 


Fey ihr c durch de® sowie t durch t—T ersetzt und 
c= 0 setzt: 
a) 2 r 
wo 00 = at) (t--1) : ] 
) ) eh ( a dc 
DOr, 1) > c (r) — > ta, e D, (kn 1) | rine (t) D,(OneT) at | Pale) re 
n=0 kz=0, 7" : ; 


Tr. 
v 


Die Gesamtwirkung zur Zeit t kann man sich nun entstanden denken durch Summation samt- 
licher, wahrend der einzelnen Zeitintervalle dt entstandenen Elementarwirkungen dc zuziiglich 
der Anfangswirkung c“(0) (Abb. 2), womit man 


co | co es "a 
O° (TW, b) = > x 


1 : | 
4@(r) — yi ee P,,(Wkn 1) | rn (r) On (Win t) dr 
n=0 L k= i 


Sa) t AG: kn 2 " r 
3 " , : : e ‘we (a) 
He > [ cE (ses De € PD, (Okn r) | r2 © P,, (Mkn 1) dr P,(é) = ce 
iE n | 


P(E) en?(0) 


n=0 « P=) dt 
u ii 
Oke lide 
nae ; > i yi [ eel ( 
G(r, Ws t) fa <a 7 ay r ce t) fear! a ®,, Onn T Cc Oye a J = ia =F e® 
2, |) EO) — D5, Palen 1) {o"0) t fe ae at 
0 


r 
wt 


< | PaO) Ba (Onn) a 


r. 
L 


P,(¢) 
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erhalt. Beachtet man nun noch die Relation 


kn \* t bn \? peer 
(a) ey aye t Ral == i (t—r) dc a) 2 ‘ al (- | (¢—T) 
de \*) + [e ee Oe) ae 
es (Ge) ME 


4 a 
0 


die sich mittels der partiellen Integrationsregel verifizieren laBt, so folgt 
r 
a 


a. a . 1 
Or, y,1) = > a = Dy Tis D,,(Oxn P) | 7? (r) Dy (kn 7) i| c (t) P,(E) 


roe NE 


50 


i 


oO es oo? t 4 ee 7S) 
kn ce ol 
= OH DS ad, a (Oinw tr) P(E) | | Tr ne'(T) On (unr) Co (r) € ( si dr dt. 


n=0 k=v 


Urs 
i 


eas ON an jedoch das erste Glied dieser Formel, da der in diesem auftretende Summen- 
ausdruc 


r 


io) 1 re 
Dy apg Palersn) [ PHi2(C) Palrsn yd 


T. 
t 


gerade die Entwicklung der Funktion 7“(r) [siehe (37) und (39)] nach den Eigenfunktionen 
®,,(@xn7) darstellt, und damit verbleibt 


2 
a2 6 
pO ek= kn 


) = 2 wr - Es (ee Fes 7) 
I(r, y.t) = >) S'S (unr) Pal) | | 72 (r) Dy (en 7) C(t) € | ‘ dr dv. 
k=v 


Onrrs 
i 


Ersetzen wir in dieser Formel c“(t) durch c{?(t) sowie 7?(r) durch 7{?(r), so haben wir die 
Lésung fiir das Temperaturfeld bei alleiniger thermischer Innenrandwirkung, und schlieBlich 
gewinnen wir fiir das Temperaturfeld infolge des Einflusses zeitlich und bis auf die Achsen- 
symmetrie auch 6rtlich beliebig veranderlicher thermischer Randwirkungen an den beiden 
Hohlkugel-Oberflachen 


t ie 


eo oO w Soe 
esp t) = >) > Gapy, Palanr) PalB) | [ * [0i%(2) uP) + 
n=0 k=0 me 
On ax 
ae; _( 0 
+ ep(t) MnP(7)] On(@enr) e dr dt 


mit 


baw. mit den Abkiirzungen 


2 = ) (G7) 


drdt (44) 


— a Ok n 


Fe deeid 
2 
Din = On eM ipa (3) | | 72 [e@(t) 4@(r) + en? (2) 7 (r)] Pa(@nn 1) & 
i) Pe 


einfach 


(To ot) = = 2 Pin D Opa ticle) ft aatt) . & = cos ye (45) 

3. Der achsensymmetrische Warmespannungszustand der dicken Hohlkugel. a) Allgemeines: 
das thermisch-elastische Verschiebungspotential. Das im Abschnitt 2 behandelte 
achsensymmetrische Temperaturfeld verursacht einen achsensymmmetrischen Warmespannungs- 
zustand. Samtliche Spannungs- und Deformationsgréfen sind somit unabhangig vom Polar- 
winkel p, so daB sich die allgemeinen Spannungs- und Deformationsgleichungen bereits erheblich 
vereinfachen, zumal die Rotationssymmetrie das Verschwinden der Tangentialverschiebungen u 
(im Sinne des Winkels y) sowie der Schubspannungen 1,, und t,,,, bedingt (Abb.3). Um zu den 
fiir die Behandlung des Problems maBgebenden Beziehungen zu gelangen, werden die Hookeschen 
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Gesetze sowie die Gleichgewichtsbedingungen am Element herangezogen. Die Massenbeschleuni- 
gungen infolge des zeitlich veranderlichen Temperatur- und damit Spannungsfeldes bleiben 
unberiicksichtigt, d.h., es werden zu jedem Zeitaugenblick statische Verhaltnisse vorausgesetzt. 

Es yerbleiben im rotationssymmetrischen Falle die Gleichgewichtsbedingungen, wenn man 
noch, um unnétige Indizes zu ersparen, t,,, = T schreibt, 


2 (sin por) +1 (sin yt) —rsiny (6, + oy) = 0, | | 
‘ (46) | 
= i" sin y 6») ser ple siny Tt) +rsinyt—rcosy oa, = 0 i 
sowie die thermisch-elastischen Gleichungen 
Se 9 fied Dats ring * pecan. 
&r = a = Fe lor—7(o, + 9,)] +o sopra an eae gy =_lo—? (Gg; + 0,) Ot | a 
47 
v 1 dw 1 Ow 1 dv w 2(1+ 7) | 
Nee tees\tetea Bp lO ae See Og) at O22 Te See ay rt, 


wobei E den Elastizititsmodul, y die Querkontraktion und «, den Warmeausdehnungskoeffi- 
zienten bedeuten, und mit v und w die Verschiebungen in radialer bzw. im Sinne von py tangen- 
tialer Richturg bezeichnet werden. Aus diesen sechs Gleichungen kénnen die sechs Unbekannten 
Or, Oy, Oy, T, v und w hberechnet werden, wobei hier entsprechend der geringeren Anzahl unbe- 
kannter Verschiebungsgréfen der Weg der Elimination der Spannungen beschritten werden soll. 
Demgema&8 werden mittels der Deformationsgleichungen (47) zunachst die Spannungen in Ab- 
hangigkeit von den Verschiebungen dargestellt 


Or 


2 


E E v (5 v 1 dw 


== w l+y 
~~ ile ETE A Ie Se Hh ar tt 2 r r peal rote y)— paged 


7) 


v Ov v 1 dw u 
( cae r oy r 


z E peat oe } x l+yv 0 
0 Tay poe cr ee Toe rar ee A re gear 0 


48) 
Feng’ v 1 dw v Ov v 1 Ow w 1+» ( 
eee eet r ow ede ee r ur r op | F etg yp) — 79,048 
a5 E Ow 1 dv w 
~ 20 +r)\ er ' + dp "| 


und nunmehr in die beiden Gleichgewichtsbedingungen eingesetzt, womit nach einiger Rechnung 
folgende zwei gekoppelte partielle Differentialgleichungen fiir die Radial- und Tangentialver- 
schiebungen v und w entstehen: | 


paee aaah oT 1? 1—2y 1 a iltapeeee 
males ra r Oy Se ete) 1 58 Ley [ete y +5 lap ae) Oe 


= —y 2r? 
1 oOo [/dv 19 ta 1 dw ow . iL skip ala Loe ipeloea ta ay (49) 
(3 rio) Pee Og Leyes | ' l—y 2r or le) — Fr Sai 


Kin partikulares Integral dieser Gleichungen gewinnt man durch Einfiihrung des sog. thermisch- 
elastischen Verschiebungspotentials F = F(r,y,t) in der Form 


9 


T 


, mor erelter 
a Oh Wen ROO a) 
Ol lity 
Ov v 1 dw 
or r 2 r ie oy ctgy =AF, 


womit die Gleichungen (49) in 


Q Ley i @) l+y 
a [AF — 75019] = 0, ene mer — () 


iibergehen, so daB wir hieraus die Poissonsche Gleichung 


Ise 
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folgern kénnen. Ist zunichst ohne die Beriicksichti 

ee p cksichtigung der Spannungs-Randbedingungen ein 
eve i, Integral der Gleichung (51) gefunden, so ergeben sich nach (48) mit Bricleinececs 
von (51) folgende, im allgemeinen den Randbedingungen nicht geniigende Warmespannungen: 


a) EGE » roi, ge her) aa 
fs = ae t 9, (Gb) 4 or Eo 
Hee 0F Tay * On] +9 ie Ge ee on) nese ‘| 
Ae) ey EES) Si ct Eo 9 o E 0/1 oF oS 
: ie \er or ay) — ety i ely. e aoe | 
/ 


tu ae Randspannungszustand; die Lovesche Funktion. Da der nach (52) be- 
as en are Warmespannungszustand 1 im allgemeinen die Randbedingungen nicht erfullt, ist 
1esem ein zusatzlicher Spannungszustand zu tiberlagern, der die infolge des Zustandes 1 hervor- 


) 
u 
4 (3 
d 
Sire Crt 
\ Les AN dt g 
t+ 2dr 
KOEN 
O \ \ 
Lea 
Dey, G \ ope 2eg \ 
/ Lie \ Ag ¥ \ 

\ 
Soe 
tale ae 

eee | / 
eee 
oy 
| aye 
me 
Nai 
ee: 
eae 


gerufenen Randspannungen an die tatsdchlichen Gegebenheiten anpabt. Dieser Zusatzspannungs- 
zustand rihrt also allein yon Oberflachenspannungen her, und seine weitere Verfolgung ist somit 
ein reines elastizititstheoretisches Problem. Die maBgebenden Gleichungen sind demnach durch 
die beiden Verschiebungsgleichungen (49) gegeben, sofern man dort 7 = 0 setzt. In diesem Falle 
fiihrt man zur Lésung der Verschiebungsgleichungen die Lovesche Verschiebungsfunktion L(r,y,t) 
ein. Man definiert sie mit 


<eote eet ates 
OD = cosy > — sin yp: — a (53a) 
in folgender Form durch ihre particllen Ableitungen 
2(1—v) Her 1 @ tor ae ee te i 
a a eT op ee OAh) 
und erhalt zunachst 
gone 6 =r pen oe or ee ae 
Bp Oe (rw) nem ey [eos > By + sin p = AL, | 
; F (53e) 
Ov v 1 dw w ley) lee | 
or sae r om r Ow a r oS = 1—2» OA ramp 4OL, 
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und damit sowie unter Beachtung von! 
Ab =A me (53 d) 
aus den beiden Verschiebungsgleichungen (49) die Beziehungen 
cosy AAL=0, 
—siny AAL=0, 
woraus wegen der Beliebigkeit von y die Bipotentialgleichung 
AAL == 0 (54) 


gefolgert werden kann. Ist eine geeignete, den Randbedingungen des Problems geniigende 
Lésung der Gleichung (54) gefunden, so folgt fiir die zugehérigen Spannungen aus (48), worin 
wieder 3 =0 zu setzen ist, mit (53a) bis (53d) 


E 
Or =~ +) 1 — 29) 


Pa ») cosp AL— 2 OL +r401], 


E i fg ctg y 0 
% = Tt) 29) pagt—— (a+ : JOE]: 
(55) 
E é IG loa il @e 
Oy = (Il +») — 2») p 4a L—2(1—m) siny- 5.4L — (+5 + api) OH], 


E i : 0 On/aleec 
1 == (i +») (1—2») la— (cos wap — sia w 5 )4E — a, (yay) OE]: 


c) Das partikulare Integral der Poissonschen Gleichung (51). Die Auffindung 
einer partikularen Lésung F(r,y,t) wird durch die Struktur der Lésung (45) fiir das Temperatur- 
feld auSerordentlich vereinfacht. Zerlegt man die Funktion F in Teilfunktionen F),, in der Form 


oOo wo 


F= DD Fans (56) 


n=0 k=0 


so kann (51) mit (45) auch in der Gestalt 


CO (ee) l oe we a 
A Dy Dy Fen aa Tap Oy 2 Din D, (kn r) ENE) T;. »(t) (57) 


geschrieben werden, und dieses ist immer erfiillt, wenn fiir jedes Glied (k, n) der Doppelsumme 


1l+y 
l—y 


Ae = Ot Di ®,, (Mkn r) igAleS)) Tx n(t) (58) 


besteht. Zur Lésung von (58) wird mit der zunadchst unbestimmten Konstanten y;, 
F, n — Vkn ®,, (Mkn r) PAS) T;1(t) (59) 


gesetzt, und da fiir die Funktionen ©,(q@;, r) und die Kugelfunktionen P,(é) nach (22a, b) die 
Differentialgleichungen 


2D, 2do, eer Geo) ieicEEDy pe 
aye me a. Ge (win — = )Pnornn Tr), r dy? + 2 ctg yp dy = 


bestehen, ergibt sich, sofern der Ansatz (59) in (58) eingefiihrt wird, mit Beachtung von (60) 
schlieBlich fiir die Konstante 


: oe gl ie, Din 
Vo = =p t 


Opn?” 
und damit fiir das partikulare Integral der Poissonschen Gleichung (51) 


ea Ite. SD 
F = >) D>) Fin =—p =u > Dies ®,, (nt) Pr(E) Tent) - (61) 


n=0 k=0 n—=0sk=0 


* Was durch entsprechendes Ausdifferenzieren bestatigt werden kann. 
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Fir die Spannungen dieses durch das th 
ermisch-elast 
Gate oe Re ote en an -elastische Verschiebungspotential gegebenen 


(1) E« 
0; the oS! BS D,, n [P,.(Mkn r) oF DP, (Min r)| P,(&) Ty n(t) 9 


n=U K=0 


Ea < 
oP — 1, BID Pan (loneed ta 


n=U k= On 


5 D (Onn ”) P,(é) 


1 dP, 
+ pare EY Gp Pann t)] Pall 


62) 
1) Ea, 1 ( 
ae ESS, ~ 2 Di, My IG n (Or n r) a Onkt P,, (Wx n ») BAe) 
n= k= 0 nk 
1 P, 

wate Ose] Tal 

Ea 
ies ers, t _ il Ae 

1—y 2 So, lao ~ Di (Wkn 7) are D1 (Wkn n| “Te Tt). s 


und speziell folgt daraus fiir die Radial- und Schubspannungen an den Oberflachen r =r; bzw. 
r=Tq der Kugel mit den Abkiirzungen 


2 a 
@) t Ea, 
di, — i, Din [®,, (Min? r;) + OF (kn r:)]5 dy a Tay Den [Pn (kn Fa) Bie Di (Wknta)] 
ef? Se E Oe By. a D),(kn ri) an ®, (Ox n ri) e® pS - Din PD (Onn ta) D, (Onn Ta) (oy) 
l—y Onn; (OnnT;)? |? es pe Ope ln (Opa) e 
einfach 
NERS oo 
oP (ryt) = >> > OF kn (Tin Ys t) DVIS eG 
n=0 k= n=0 k=0 
foe) e.8) co 
6(ry.1) = > Dy ohnlr rast) =— >) >'d® PAE) Thal), 
n=0 k= n=0 k=0 
ONY, (64) 
(D/», seth Lae dP, 
TE (rin pt i >) Tia iz, 5 Dee ve 7 Chin 7 {hae 
o8) ioe) 


) psy ees 
tn (72,4), t) = De (Tas > t) aa ie Haye dy De n(e)2 


wT kat i 


d) Der Randspannungszustand zur Erreichung randspannungsfreier Ober- 
flachen. Im folgenden wird es sich darum handeln, mit Hilfe eines bestimmten, zu tiberlagern- 
den Randspannungszustandes 2 resultierend einen oberflachenspannungsfreien Warmespannungs- 
zustand zu erzeugen, d.h., die durch (64) gegebenen und an der Kugel ein Gleichgewichtssystem 
bildenden Oberflachenspannungen zum Verschwinden zu bringen. Dieses kann man erreichen, 
wenn man der Oberflachenspannungsfreiheit bei jedem einzelnen Glied (k,n) der in (64) vor- 
liegenden Reikenentwicklungen Rechnung trigt. Die nachfolgenden Ausfiithrungen beschranken 
sich auf ein Glied (k,n). Die Verallgemeinerung auf siimtliche Reihenglieder bietet dann keine 
Schwierigkeit. 

Die maBgebende Beziehung fiir die Untersuchung eines solchen era ee ee ist 
die Differentialgleichung (54) fiir die Lovesche Funkiion LE(r, y, t). Zur Lésung von (54) wird 
zunachst 


AL®. = L2(r,y, 1) (65) 
gesetzt, womit (54) die Form 
AAL® = A (AL@) = ALE = (66) 


1 Die Striche deuten Ableitungen nach dem Gesamtargument (Pr an. 
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annimmt. Fir L®, (r,y,t) wird nun ein Produktansatz in der Gestalt 


ERC, y.t) = Teal) S [phir et + pr OF] PE) (67) 


u=0 


benutzt, wobei die p,, und p,, willktirliche Konstanten sind. Dieser Ansatz befriedigt wegen 


(a nie 2 ae (Pp n) pit 2 pS" n) p— (w+ V\ = ESS we ple (ke, ny pl Seay) n) p— (ut ”) 
jas ie r 


und 


1 /d? Pu : dPi\ u(u-+ 1) 
al dy? + ctgy z) am r2 Pa 


die Differentialgleichung (66) identisch, und man hat demnach fiir die Lovesche Funktion nach 
(65) die Differentialgleichung 


AL® = Tun(t) 3 [phe r# + pe") r— +0] P(E). (68) 
u=0 
Zu ihrer Lésung wird mit den noch unbestimmten Funktionen RY” (r) der Ansatz 
Li2(r, pst) = Tent) SY Ry’ (r) Pe) (69) 
u=0 


versucht, und indem man diesen in (68) einfiihrt, kann man fiir jedes einzelne Reihenglied yu 
A [Re %e) Pe = (ols + pM) Py 


folgern. Hieraus erhalt man schlieBlich unter Beachtung von 


il (GPZy dP, u(u+)) 
be ( Peet gy iy) eee te 


folgende gewéhnliche Differentialgleichung fiir R{")(r): 
(k,n n 
CNT wut) p 


dr? hi i: r 


(k, eG \= Pie 7 ee Dae n) pe }-1) (70) 


Thre Lésung setzt sich aus dem Fundamentalsystem der homogenen Gleichung, woftr mit den 
Wicwarli k,n kyr 
willkiirlichen Konstanten ne und Pye ) 
k,n sn) »—(ut+1 
pi dre a Pe ) Took ) 


geschrieben werden kann, und einem partikuléren Integral der inhomogenen Gleichung (70) 
zusammen, das mittels der Methode der Variation der Konstanten gefunden wird. Man erhalt 
schlieBlich fiir die gesamte Lésung der Gleichung (70) 


RY n) See ) poe po —(u—1) (k, Jee lt (k,n) —(#+1) 
pS eta) On, RE ae ee 
so dafs fiir die Lovesche Funktion L;,,,(r,y,t) nach (69) 
0 pe” e pt n) z 
EEG a 2; sae Re eG ee 


(71) 
hervorgeht. Nach (55) kénnen nun mit (71) die Spannungen angegeben werden. Man erhilt 


beispielsweise fiir die Radial- und Schubspannungen nach einigen Rechnungen und unter Be- 
achtung der Rekursionsformeln! 


(@ +1) Pld) —2 (e+ 1)EP,8) +a Pea) =O, =) Pa ee ey 
§ Pi) — Pi®) =" PE), Phinlé)— Pia =(2u+ 1) Pe), 
Prsslé) —E PAE) = (u +1) Pe 


' Siehe z. B.: J. Lense, Kugelfunktionen. Leipzig 1950. Die Striche deuten Ableitungen nach dem Argu- 
ment €=cosy an. P_,(&) ist Null zu setzen. 
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schlieBlich mit den Abkiirzungen 
(k, n) 


E(u +1) Py”) 3u+5 9 u—l P3y, 
(l+ 7) (1— 2) onli 2y)r + 3enpy 2” — (“4 +1) (ut 4)) ro) 


+ py”) + 2) (wu + 3) r- = = feri(r) , 


Eu Bie } 
i one a LW ye (k,n w+l a 
rics ee 7 (29 NE 8) TP Peril eon asie i 
=H — 9 W— 3) 0-4] Foy, 
E noe +1) ay 2 
‘ ath bees (He Ne Pra (4 + 1) ra (wet? 
(l+%) 0 — 2») eal =) yc ie — eee ae 


ft 


pt) (u el )(u =e 3) (u4 4) ae — g(r) : 


- k,n) 
E Pie. Th i pt] 
Py ! us k,n be ( ie —(u 
des) d= 29) Bee ¢ ”) os 0 loge te sf al ae 
Pee (we — 2) | = g(r) 


| 
ae 


folgende Relationen: 


2 wally, Fk, n} 
Orkn = Tralt) >) (FOG) Parr +fO%) P,-) 


u=0 


= Tal) [79 Pe + DTI +A) PA 


| 


oO 


P _ CH pts = avin CN pas 
T, ca Ty n(t) aa - a) ae iv ye (7) oF 


u =i 


a 


dP,, 


Se Tanlt) 2 gin) sie OU) eae 


Indem man nun von diesen cee hen des Zustandes 2 fordert, dab sie an den Kugelober- 
flachen r=—r; bzw. r=r, die entsprechenden Reihenglieder (k,n) der Spannungen des Zustandes 1 
zum Verschwinden bringen sollen, erhalt man mit (64) die Gleichungen 


rig "\r Dx Y k, n)} (k,n 
fe i) eae fe a ( 4 FEM; EF = di) P| 
ot} en ee of} rie d.h. oo (74) 
n +> 1/ k,n 
f® ra) Py (fs tee eG) i> fen ‘(r,)) le) = = a n 1a y) 
i= 
a k, n) ob, dP x (i) dP, 
P be is + gt, bl (ri) )) “dy re, n dy g 
aad) anys a eae (3 
n k,n [Pets (a Gin 
Srna sinivayte antl 
b= 


und hieraus kann man nach Vergleich der Koeffizienten der Funktionen P,, baw. deren Ablei- 
tungen zunichst aus (74) 
1. fir den Fall n = 0 


i fe 
fOr) = ath, eda 
+f 


lk, n) 
furl) ta -4 Uti ai hus OO 
lk, n 

eee fe 0 


(76) 


a 
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2. fur denFall n+ 0 
fo) =0, PAG )+fS20) = af, 
Fer) =H0, FEM) Pad) = dee, 
FN) + Sic % ri) =0 
FENG) + FEN) = 0 


(77) 
CUPS eile 


und aus (75) 
o (k n = n t n =(k,1 = a 
ce) g(r. )ee ete gi i Mri )= ey) zi ge (ra) 35 Bn4 | (Ta) rr ey 2 | 
wri) + Be Mn) = (78) 


BO Nre) +.B Mrs) = 


folgern. Aus (76) und (78) bzw. (77) und (78) sind mit (72) die Konstanten pes) (it 
uu = 0...) eines jeden Teil-Zustandes (k,n), je nachdem man die Reihenentwicklungen (yw) 
abbricht, mit beliebiger Genauigkeit ermittelbar. Die Lovesche Funktion Lj? ist damit nach 
(71) bestimmt, und nach Uberlagerung sdmtlicher Teilzustaénde auch der resultierende Rand- 
spannungszustand 2; 


0 
fir un 
0 


[OS SOS es (79) 


n=0 k=0 


4. Der kugelsymmetrische Fall. a) Das Temperaturfeld. Die thermischen Randwirkungen 
gi=Gi(t) und a= Galt) sind 1 in diesem Falle von wy bzw. € = cosy unabhangig, so daB man 
aus (43b) wegen 


Gives WNP sis (NA a8 


sofort 
ce? =e = 0 itp pesky Me gos Ce Wilt), 8c. (t) ==, 8) (80) 


erhalt. Somit wird nach (44) T;,,(t)=0 fiir n=1, 2, ..., wodurch in (45) die Reihenentwicklung 
hinsichtlich n fortfallt. Es bleiben allein die Glieder (k,n=0) tibrig, so daB man aus (45) schlieB- 
lich ‘mit D,o = Dz, O10 = Ox, Tee = Tz sowie mit P, — 1 


Or, t) = p> D;, ®y(ox 7) T;,(t) (81) 
=0 
gewinnt. 


Hierbei sind die Eigenwerte wzo = @, im Sinne von (32) durch die Gleichung 


AGG) = a, Ji (@ ra) + 5c (2 WTaJ 1 (@ Ta) — J 1 (w ra) 
2 ‘ a) 22) 


3 eee 


2 


Sé [6 J i(@ ri) + gt (2 ornJ s1(@n)—J i (@ r)| 
é tay, 


= la J (or) (2 Gita (or) Ts kes rs) 
2 ¢ epee 2 


2 2 


x G Ji (mr) + y (2 WT; J (w ri) — Ji (@: )) = 0) 


bzw. mit den Relationen! 


aoe 
Ji (or) =V. sinor, LC Nes 
2 


* Siehe z. B.: R. Rothe und I. Szabé, Héhere Mathematik VI. Stuttgart 1953. 


cos wr 


XXIV. Band 1956 Trostel: Instationire Warmespannungen in einer Hohlkugel 389 
und damit 
Pa 
@D = =a sin or 
o(@ r) Se [@; Ta COS & Ta — Ay (H Tq SiN W Ty + COS OT i)| — a 
— [a,r,sinor, + Ay (@ rq COS ® Ta — Sin ra) | ots “) (82) 
Or 
durch 
tee ri Ta (ay by — a, by) — - Ay by (ra —1;) 
‘ i © (a,b, b, + a, by b wo) Fi Ta 1 Gy.6, — a, b, r, —a, 6, 7; BD 
festgelegt. Mit (38) erhalt man aus (44) 
2 w? 
D,o = D, = ek ; : (84) 
a I" (2 + 7 + —@, 2%)| sd 
OO), Tr uF 
und auch die Zeitfunktion vereinfacht sich. Zunichst sind 
s ‘ r M (0) (a) 
nO) = af + BP r- 0) =o + AE, ype) = af + 
bzw. mit (18) 
ee “s B 
ial ihoslarn 
n(r) = yr) = — ee 
a, by ra 7; (Ta — 71) — G, bg 72 + ab, 1?” i @, by ra 1; (Ta — 1) — & bg 1G + Gy OF; 


und damit folgt aus (44) unter Beachtung von (80) 


4 bye o(t) = T;(t) 


SC i eo oy 


i) (2) 
drdr. (85) 


| | Ge =k 

oe = r WET) e 
G0, 7, 14%, —7,) 3, 0,72 a,b, 12 o( ) 

0 rj 


b) Die Warmespannungen. Daim kugelsymmetrischen Falle nur Radialverschiebungen v 


auftreten kénnen, erhalt man aus (48) mit w = 0 und >= == (( 
a E CU a 2) ee Eo, o 
Peay ee r| F300 | & 
E v Ov E 4, 29 | 
Ce OS (eer). 29) ei a a ee 
wobei fiir die Radialverschiebungen v selbst aus der ersten Gleichung von (49) — die zweite 
ist im kugelsymmetrischen Falle von selbst erfiillt — die Differentialgleichung 
ad{l1 oa l+v_ ad 
oe ieee @ = Sey a ey) 


hervorgeht. Fortlaufende Integration liefert mit den willkiirlichen Zeitfunktionen C,(¢) und C,(t) 


} Gs 
v= Eee a r? (r,t) dr + C,(t) ee =f y ; (88) 


1—»v r? 3 


und man erhalt schlieBlich damit nach (86) fiir die Spannungen 


2Eo 1 ert 
= — 5 : = rod + Eley es ey | | 


Ea ee G Cre 
Cn Coe alee: [0 ar 7 Elsqtay + eee 


(89) 
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Die Zeitfunktionen C,(t) und C,(t) werden aus den Randbedingungen o, (r;, t) =0 und o,(rq, t) =0, 
bestimmt, und man gewinnt letztlich nach einiger Rechnung 


Page es AS] 
ea aaa Poa} | eoeoa—f—(4F] [og uf, 


r 


con Eafe —atalleye feo ep | feo | 
—p i. I al 


(90) 


ri r 


c) Naherungsformeln fiir den kugelsymmetrischen Fall. Entsprechend den Aus- 
fiihrungen des Verfassers im ersten Heft! des Ing.-Archivs 1956, S.1ff. kann man nun auch hier 
fiir den Fall nicht zu schneller Anheizvorgange an diinnen Hohlkugeln Naherungsformeln fiir die 
Spannungen herleiten. Insbesondere erhalt man aus (90), wenn 0;(t) bzw. #,(t) die Innen- bzw. 
AufSentemperatur der Hohlkugel, R, ihr mittlerer Radius und h ihre Wandstarke bedeuten 
mit A(t) = I(t) — 3,(t) fir die extremalen Radialspannungen 


min o,(t) = ae — : i (91) 
und fiir die extremalen Ringspannungen 
1) am Innenrand 
min ¢,,; = min o,, = — a ae (1 | = oo ; (92) 
2) am AuBenrand 
Macon, = Max, a 3a 2 (1 ; ae . (93) 


Aus (91), (92) und (93) ersieht man, daB die Radialspannungen im Vergleich zu den Ringspan- 
nungen klein sind. Bei Anheizung des Innenrandes z. B. treten die gréBten Spannungen in Form 
von Druckspannungen am inneren Kugelrande auf. Sie sind, da max 430 >A#, (A 0, = stationare 
Endtemperaturdifferenz), stets gréBer als die entsprechenden Spannungswerte am Ende des 
Anheizvorganges. Uber die Ermittlung von max A? orientiere sich der Leser in dem bereits 
zitierten Aufsatz ,,Instationare Warmespannungen in Hohlzylindern“. 


Pr (x,t) 


Bele?) 


ly 
EEN 
Noy Yt) =Oy -r 


Abb. 4, 


5. Achsensymmetrische Probleme an diinnwandigen Hohikugeln. Um die durch (64) gegebenen 
Oberflachenspannungen zum Verschwinden zu bringen, genugt in diesem Falle ein zum Grund- 
zustand I zu tiberlagernder Membranzustand 2. Hierunter versteht man einen jeweils langs der 
Schalendicke h nicht veranderlichen Spannungszustand, dessen Schnittkrafte, worunter die 
Spannungsresultierenden je Schnittlangencinheit verstanden werden, jeweils in die an die Schale 
gelegte Tangentialebene fallen. Bezeichnet man mit PY, t) bzw. Pi (y,t) die auf die Schale 


* RK. Trostel, Instationére Warmespannungen in Hohlzylindern mit Kreisringquerschnitt. 
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wirkenden auBeren Radial- bzw. Schubkrifte je Flacheneinheit (Abb. 4), so erhalt man aus den 
Gleichgewichtsbedingungen fiir die Membranspannungen op baw. of schliefBlich 


y 
Gy (Yt) = = a | tee ) sin ¥ cos ¥ — pr(z, t) sin? x] dy, | 
tom (94) 


” 
oP (p,t) = aie [Pr'(x-8) sin x cos 4 —~ p(y, ) sin’ | ee | 
4=0 


und indem man hierin gemaf (64) 


POY, t) = — [oP (ra yt) — OP (ry, 0] = ey: (di, — din) Pal) Tal 


— 


; (95) 
; dP, 
Pi (y, t) =>— lee (res W, t) anal aa (r; YW, t) | = s > (en — exin) dy Ti, n(t) | 


n=1 
einsetzt, sind die Spannungen des Zusatzzustandes 2 gefunden. Hinsichtlich der Verteilung der 
Radial- bzw. Schubspannungen of und t+ langs der Schalendicke wird man wohl mit geniigen- 
der Genauigkeit jeweils ein lineares Gesetz in Ansatz bringen diirfen. 


(Eingegangen am 28. Februar 1956.) 


Anschrift des Verfassers: Dr.-Ing. Rudolf Trostel, Berlin-Charlottenburg 2, 
Technische Universitat, Hardenbergstr. 34, Lehrstuhl fiir Mechanik. 
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The Equations of Motion for Curved and Twisted Elastic Bars deduced 
by the use of the ,,Method of Internal Constraints“ 


By E. Volterra 


1. Introduction. The object of the present paper is to show how the ,,Method of Internal 
Constraints‘, which was used in a previous paper! for deriving the equations of motion of 
curved elastic bars, can also be used for deriving the equations of motion for curved and twisted 
bars taking into account the effects of shear and of rotatory inertia. 

The problem of vibrations of curved and twisted bars finds many important practical appli- 
cations among which are the vibrations of helicoidal springs and the vibrations of turbine buckets. 

The isotropic homogeneous elastic bar is supposed to be an elongated solid with a curvilinear 
baricentric axis. It is referred to a system of curvilinear right-handed coordinates (0, X, Y, Z) 
in which the origin 0 coincides with a generic point on the central axis of the bar, X with the 
central axis (and has as positive direction the positive direction of the central axis), Y coincides 
with the principal normal and is directed positively towards the center of curvature, and Z 
coincides with the binormal at 0. 

By assuming that during motion the sections of the bar originally normal to its axis remain 
plane, the following condition of constraint: 


V(x, y, 23 t) =9(x, t) + yA(x,t) + 2 p(x, 2) (1) 


must be satisfied. 


In equation (1) V represents the elastic displacement vector at a generic point of the bar, 
while ¥, A and yw are three vectors, a priori unknown, functions of the variables x and t. Calling 


ui, A‘, u', v' (it = 1, 2, 3) respectively the controvariant components of the vectors Vey andy 
on the X, Y and Z axes, from equation (1) equations (2) follow 


ul(x, y, 23 t) =91(x, t) + vy At(*,4) + 2 (x, 1) ; 
u?(%, ¥, 23 t) = v(x, t) | yA?(«, t) + 2 y7(x, 2), (2) 
u(x, y, 23 t) = v(x, t) + y A%(x, t) + 2 w3(x, t) . 


The geometric significance of the internal constraints expressed by equations (1) and (2) 
has been explained in the previous paper. In the same paper, also, the square of the line element 
ds* has been derived in the more general case in which the central axis of the bar has both first 
and second curvatures, which will be denoted respectively by the letters c and t. It has been 
supposed in deriving the expression of ds* that, in a first approximation, the second powers 
of the quantities (pure numbers) yc, zc, yt and zt can be neglected in comparison with unity. 
The equations of motion which shall be derived apply then to the general case in which curvature c 
and twist t are given functions of the variable x, provided that 


dR(x) ; dT(x) 
dx and ax 


[R(x) and T(x) being respectively the radii of curvature and twist of the central axis of the bar] 
are finite quantities. In fact, in deriving the equations of motion through the variational principle 
in the given order of approximation the terms 


do(x) ad dR(x) 
dx [R(x dx * 

dt(x) 5 1 : dF(x) 
dx KGCIP Gk 


can be neglected. 


1 See E. Volterra, Ing.-Arch., 23, (1955), p. 402. 
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From the fundamental (Riemannian) metric tensor it is possible to derive the covariant and 
controvariant components of strain and then it is possible to derive the expressions for the 
Kinetic and Potential Energies of the bar in terms of the nine unknown controvariant components 


i yl, v2, v3; Al, 22, 23; wh, we, 
of the three vectors », A and w. 

The application of Hamiltons Principle furnishes the nine differential equations of motion 
and the boundary conditions which univocally determine the nine unknown functions y!, »?, y3; 
AM, 27,13; ww, u?, wu? in the case of free vibrations. 

In deriving these equations of motion it will be assumed, for the sake of simplicity, that the 
axes Y and Z coincide with the principal axes of inertia of the cross-section of the bar and that 
the inertia characteristics of the cross-section of the bar are constant, although the method of ana- 
lysis presented in the paper may be applied also to the general case in which the inertia charac- 
teristics of the cross-section of the bar are known functions of the variable x. 

From these general equations of motion the equations of motion for the particular cases of 
twisted bars (c = 0), of curved bars (t = 0), of straight bars (c = t = 0) areimmediately derived. 


2. Expression of ds*; covariant and controvariant components of strain. The expression of the 
square of the line element ds? in the orthogonal curvilinear system O X Y Z is 1 


ds* = dx* (1 — 2 ye) + dy® + dz? —2 yr dx dz 4+ 2z2tdxdy. (3) 
The fundamental tensor will in this case be espressed by the metrics ?: 
[ay =1—2 ye Ch == AU ag = —2 yt | 
(Ay, = 227 Ag, = 1 ay3 = | (4) 
ae A dz. = 0 O33 oe | 


From which it follows with the same degree of approximation that 


a=I1—2yc, ja=l1—ye, ioe (5) 
la 
jat—1l4+2ye a® =—22t al3 = 2 yr | 
= = 2 27 Gp) = a3 — 0 | (6) 
give 2 yr Ch gail) Ge 


The covariant components of strain are expressed by * 
1 
yh rs (wij + Uxji) (7) 


where, due to the definition of covariant differentiation of tensors * 


StH " 


( ‘| representing Christoffel’s symbols of the second kind. It follows that 
J 


Ail im 


a Se r ay; SUH Bee 
i ‘ H a asta ax" ar Ox" 
and that : 
; ns arin oA _ Ot yt 
Yu =—e (A? ¥ [ez a) (1 —2 ye) (y ax r 3 dx 5x) 
ae ou? Ov Pein ay 7 ee 10 
+ 2ar(y ax * Bx ol ead Ox ree ee Gra She 


Yo, = 227A + ,?, 
eg = — 2 yt + Ls 


1 E. Volterra, L. c. pete 
3 the Bacio used in the text is the same as the one used by T. Levi-Civita in his book: The Absolute 


i ial Calculus, London 1926. ; ; 
Eh aay Viton Rendiconti Accademia dei Lincei, 20 (1934) p. 424, 463; 21 (1935) p. 14. 
4 T. Levi-Civita, L. c. 
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yin = TAP y ap pee 8) ~(1—2 ye) A + 202? yt A 
pee Ot Out =) aml r(y ae Ou? 2), 


Y Ox res dx | Ox) | Ox |” Ox Ox)’ 

P 1 / 

Vis =—tAPy + wre +) 4 oy (iL = Daye) ae Sra eee (10°) 
on. , oe ov! 1 {aA ae avs 
—ye(y ax r ” 0x ! ‘a ! 2 Ce fay tp a 


| ay LE 
Yes Bi i Ae 


The controvariant components of strain can be derived either from the formula 


Sue Pha . 
Mise == = (uil® oe uk/*) s (11) 
or from the formula 
yr= Si aiathy,  G,h=1,2,3). Re 
3. Expressions of the Kinetic and Potential Energies. The Kinetic Energy per unit volume 
of the bar will be expressed by 
eats om ete s,r =1,2,3 
P= 5 ae: SEO Ob. 
o being the density of the material of the bar, or by 
Ae 7 (eu oe Ou? _,ourdue , (dur _ /due\? 
Ty = $290) (5 ales pear rea akan Gt e + (5 ce!) 
The Total Kinetic Energy of the bar will be expressed by 
T = [Ya T, dx dy dz (14) 
V 


V being the volume of the bar. 
By substituting into equations (13) and (14) for u!, v2 and u® the expressions given by 
equations (2), by performing the integrations with respect to the variables y and z, calling 


if estes dy dz, es =a 2 dy dz 


the principal moments of inertia of the cross-section of the bar and A the area of the cross-section 
of the bar, one obtains the following expression for T: 


aon OAR\2) Our\2 GN OA Oyt 7 oe ov? 
pute. + 1,(F 7] 3 4(5, Cire ae ees, 
i 
ue Ovi 0A dv 03 Ovi 0A? \2 
~ A fun EO Se Wey ee Ses [ae 
ue “OL sO pol oane Ot ot *\ Ot (15) 


Ou O\2 OA? Ov? 0A3\2 Ou3\2 
[ | a 
1,28) A(3) , ae (57) ! 1,(3) 


+ 4(@) 2 rel 2 a 


at at 
I being the length of the bar. The Strain Energy per unit volume of the bar will be expressed by 
dh rap 
OO + Gy (16) 
where @, divergence of the elastic displacement vector V, is given by ! 


5 i : 
ee he _ 1 ova “) Cee (17) 


i\dx" Va @x' 


while wy, tensor of sehen is expressed as a function of the covariant components of 
strain by 
S: 45 
Ve na) AOpen (18) 


1 T. Levi-Civita, L. c. 
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S, ay . . . . 
Ais Lamé’s constant and G the modulus of rigidity of the material of the bar related to Young’s 
modulus E and to Poisson’s ratio ¢ by the relationships 


ee oH 7 =e E 
SS CSOs aa aera 


The Total Potential Energy of the bar will be espressed by 


Wice= | Va ee | Cy) dx dy dz. (19) 


V 


By substituting into equation (19) the expressions given by equations (17), (18), (2), (5), (6) and 
(10), by performing the integrations over the cross-section of the bar, one obtains the following 
expression for U: 


= A OA1\2 i Oyt\2 jh avl\2 
Se Ce ee 
1 
A : A opt Ovi 
+(z + aca + (5 +6) Ary -AABE + AWS + AAI 
G 0A2\2 GG Ou? Ov?\2 
+S A(a)? + 1 =| : 1,( a) ! A 5) 


G G 
2 2 

ow  G G 03\2 iG, C/u\2 
G 


Ox Ox 


2 3 
: (2) + GA 4% 4 


é 


ua + 40s) 
c 1,7? es 


2 Ou 
—(A+ 26)cl,u Dx 


1 
Bn (jo 26)cAe% Je AP AeA pi — (A+ 26) elas 


: an 41 O22 jo ot 
gee (a 1 oe Oe 
grt iw Celso Gels 


l 
( 
oA 
= ee 
Get las | 
rary ed pees ba 

= CL a a re 2 oa A } y Ae | 
on Om? 3 OP 

+ 267t LBW +26 Ty pea + 267A vo 


Qt Ov? 
Ox 


1 Ayyz2 
Be rt oere 


a GJ 
—26rI,#7#— +26rl, Ae ag 


Ox : 
2 Oe? 


Ox 
j1 Ch 
+ 26t 1,3 — 261 LP a —2G6tIyu Ag 


: 3 O28 dat av’ 
== OG Al ye MENTE 1,3 — 26% Pa er 


oe Cle ieee oF 2671, 0 | dx. 
267 1,550 + 264 tsp ae ahs y ag J 


4. The Equations of Motion for Curved and Twisted Bars. On the assumption that no external 
forces (body or surface forces) are applied to the bar, by applying Hamilton’s Principle under 


the form 
ty 
6f(U+ T)d& =0 
ty 
to the whole bar supposing the displacements to be zero at t, and t,, one obtains the following 
) 


i i ion i i 7 tiens: 
nine equations of motion in the nine unknown func 


4 2 3. Tomato. 
PAs Aee [Eases Poo) 
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of the variables ~ and t: 


2, r 
(1+ 26) A822 _ (a4 26)er4 2678 
» OA? Ov ; 1 aoe ous 
lan ae A EA ae a 
b Zaye 
—CP 2619 aor — 300 —2gen oe 
p 3 
6nZ® + Ce ee (2a + Gert 4 4Ger i 
x 
2 Om , Ovi ‘ ; A 077? 
+261 (ry 4 ee es (A + 26)? —Ap® + he =or ae 
ae ee 
z or ” 
2 0713 2 Ov x 4 a0 Vs 2 OW 
G sag te Vs Hea 2G6 lea 4G Ts a. 
: ane a 243 
Sk CU(ry Ts) ae OF raae OH CN OF am: 
Aipe 2y1 
Sir ae ae 
an 
(A+ 26) 4 26cr rat 26) 267 (re 4 
— ov? ¢ ee 1 2 Out 2 Oy? 
+Gecr oe 26% CU 0 ty, i aot yee Ae 
3 
Gy oS SOME ee es 4) a ea t (21) 
a AP yp. Q2yl 
ery re pe — GI =o SE + 2orjt—, 
Gr} 4 9G r (rt 4+) —46cr Ne den 
yt Ae 
Ax (A+ 26) uW2=—APR ier 08 none 
ol > OAL we 5 OFA 
2 ay 241 
reo Je = 426) agent 
OP? 2 O7A3 
+ 2o0try 52 — 20TH aS > 
2 2. 1 
62+ Geri 426% ioe + 467% — +6o 
x 
2yy2 272 
2068 2a bdo ter ak griees 
72 Ow 
20T Ty ee 
ay? e aryl Opie Out 
are pear e Olsens Opener | 
ary? 2 0743 SOzAL 
eT AL =o = Orem 20rT ae 


with the following boundary conditions at x = 0, x —I: 


an one 200 
499 (7 4G) er 2 2 Car ra eee 
Pee boa 42443 ed 67 =" 


an ae" | 
met eget 2epi— 272 — opt] di = 0, 
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ja+26)% +267™ (A+ 26) ep? 267] ut 0, | 
ea (22) 
= 2% pi —27 a3 One a 


a+ 26) — A+ 26)¢ re ee 


EZ (b2 Grey Aa & ie byt 0, 


2. 

B one blr i + 2xy8 drt = 0, 
3 3 

E ones Qo eet wt —2 rr] bye = 0. 


In these equations, together with the nine unknown functions 
ESE Sa VUES TRAD a erg iki a 


of the variables x and t, appear the radii of inertia r, and r, of the cross- -section of the bar which 
have been supposed to be constant, and the first and second curvatures c and t of the central 
axis of the bar which have been supposed known functions of the variable x. 


5. The Equations of Motion for Twisted Bars. In the case of twisted bars c = 0, and the above 
equations ae and (22) become 


AL267 o— 264 Bel 4 26r(9 + 1) en — 261 =o, 
(A+26)r ps or 26x (rh + 1) 6 4 2618 6 
a Bor , | (23) 
Bee ee rot Gr = OC CH 4 26 ryt =e + 2¢ poe 
eee o 6% 267i =9 Qo rte 
and 
ene 4 4crn™ 4 2678 +H AB 4 26) Apt =o Ze 
on 2h 260i —46r nS 4 26x (1 + ye 
— 6p —C =o ri 2orit oe, 
Gri HH 4 26x 26 c(h + & 
$4673 Gy ee a os 
ert Fer (rt +) B46 cr BAS A + 26) pA = 07 Ge 
P4208 p 261g Se 26rn Se + aR ie 


yl Oru sy OPI 
= Ose t 20T ty Ge 20 ts ae 
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with the following boundary conditions at x = 0 « =I: 


(A 4 26) 5" 2G ri —26rryl| ol = 0, 


| 


[ae +2622 — 264] jut = 0, 


sas om (25) 

Se Poa pies ora an 2:0 yoo i On, 

es as vend ae p27 dv? = 0 

| Ox Ox 

and 
ge t Beat + 2129) 5 = 0, 
Ov} 

ace tarps — 2422] 818 = 0, 

2=0, (26) 


ee 


a 
Ox 
ee ot 27 u2— 2A?! 6 
ay? 
E 
G +26) 42673 _26rn@ 4 12 +4 p5| oo =0. 
As equations (23), Ce. eo) and (26) show, during vibrations of a twisted bar the functions 
A}, w}, v? and v3; A?, 2%, u?, uw? and y1 are coupled. 


6. The Equations of Motion for Curved Bars. In the case of curved bars t = 0 and the above 


equations (21) and a become 


+2672 _@t 26)c ge 
Ov? 72 Ou! e241 O2y3 
C= oe Gi=or 25 —3ecri-, 20 CT ae ; 
2 2? 2 Oy? oe Ovi 
Cri + Ger; err te (eA G) ie elas 
6 : 02,2 : O2y2 
yy 2 Pi esa pie ee eee 
(A ae se a +ic1 OT; One OT:€ aes 
1 24,3 
ov » 2u (27) 


Gr oe 44 og Ea 8 Cs ee oe 
a ue SP Be ONY On? 


Cyl o2A1 72 3 
(A+ ne 7 — A+ 26)e SS 41S 4 ee 


Ox? 
Ov? 0?y1 a2 Al 
Sa (ee. CPi g ats 2 
2) yg OG mae aieeaee, 
1p? ea aM 
Cres 2 Gie ie pee 2G6)e + Ae 22 
5 Oy? 5. GE 
oa eS | a8 Spi ae 
Sp Ae 0a OTC, J 
and 
a2A3 . Bee : 2ayd 
Gri ae + Gen Cen Gu? CM ore or ~, 
Fa 
(A+ 26) 7S +26) eM + Gene = G0 Gee eee 
x t 
2 OP? 2 Out . Oru (28) 
G rym + (A+ 2G) ery ae Cae GAP =o eae aR ; 
0’y3 0713 Ou o*y3 2 OAs 
G Si 2 oe eS iy é 
ak Ole erat Oe OG Ona 


Yt eee 
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with the following boundary conditions at x — Vit sel i 
| A+ 26) F 4 A+ 26)0 mM + erty 4200 + o)ertie] on, 
Bs s oF ea é2 = 0, 
is Ou Oe (29) 
v b 26) _ 4426 ee) ae Ans dt =0, 
ee afar 
and: 
lee Fe eH] 82° = 0, 
E — c8 Our = 0r, 
(30) 


As equations (27), (28), (29) and (30) show, during vibrations of a curved bar the functions 
2,2, 8,9? and v?; 23, ut, u? and y3 are coupled. 


7. The Equations of Motion for Straight Bars. In the case of straight bars c = 0, tr — 0 and 
the above equations (21) and (22) become 


oral Ov? o2A1 
+2 a Ree Pe 
tee Ox? 2 Ox Ch = 0's ae | 
C ey? Ch ee Ove | 
‘ 0x2 | Ox. = on ” 
2 a2 Ovi = > O22 
ae 
Oru Ovi a 23 
G y Ox2 A Aye (A 26) u3 —))? =O ye ; 
ta aw ah 
(A + 2, G) On t A ae | A ve oe? 
, 2 OAs 243 
Gr; A? Gw— GA = Le | 
2/2 
op oe a | 
"OL 2 
2 Out ov : 07 
aoe SC cera: 
OP OEE OY 
dt | 7 ax © Oe 


with the following boundary conditions at x = 0, x =I: 


| 


i = 
| =U, 
| 5,8 — 0 
Bene 
[¢ P26) 4 aR 4 
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(31) 


(35) 


(36) 
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(37) 


wn | 
| =f | Oy? 0: | 


As equations (31), (32), (33), (34), (35), (36), (37) and (38) show, during vibrations of a 
straight bar the functions: J! and y?; J?, w® and v1; A? and w?; uw) and y? are coupled. 

These equations are identical with the equations which were derived in the previous paper 1. 
In fact, in this particular case the controvariant components 41, J?, 13; wi, w?, uw; vt, v? and vy? 
coincide with the cartesian components /,, Ay, Az3 Ws ys Uzi Vxo Vys ¥z Of the elastic displacement 


vector V. Equations (31) and (35) define the flexural eigenvibrations in the XY plane of the 
bar, equations (32) and (36) the longitudinal eigenvibrations along the axis of the bar, equations 
(33) and (37) the torsional eigenvibrations, and finally equations (34) and (38) the flexural 
eigenvibrations in the XZ plane of the bar.? 


(38) 


(Eingegangen am 15. Marz 1956) 
Anschrift des Verfassers: Prof. E. Volterra, Troy N. Y. (USA) 1567 Tibbits Ave. 


1 EF. Volterra, L. c. 


The material in the paper 1s based on an investi ti y Cc itute 

ga on, conducted at Rensselaer Pol techni Inst tut 

Troy, New York under the sponsorship of the Office of Or nance R C 0- 1 -O D 709 
ry d n esearch, Contract No. 3 il 5 R ’ 
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Uber das Grammelsche Verfahren bei Kigenschwingungsaufgaben 


Von K. Zoller 


1. Einleitung. Fir Eigenwertaufgaben, wie sie bei der Ermittlung der Eigenfrequenz von 
Schwingern mit kontinuierlich verteilter Masse vorliegen, verfiigt man iiber eine gut ausgebaute 
mathematische Theorie!, die sich vorwiegend auf die Differentialgleichung des Problems stiitzt 
und die u. a. auch zu den Naherungsverfahren fiihrt, die iiber das Hamiltonsche Prinzip und 
den daraus herzuleitenden Rayleighschen Quotienten etwa in den Galerkinschen Gleichungen 
baw. uber die quellenmafige Darstellung in den Grammelschen Gleichungen gegeben sind. Vom 
Standpunkte der Mechanik aus befriedigt dieses Vorgehen aber, wie schon eingehend untersucht 
und begriindet worden ist?, nicht voll, da diese Theorie, wenigstens in ihrem gegenwartigen 
Stand, schon einfache technische Schwingeranordnungen nicht umfaBt, wahrend die Mechanik 
immer noch ein wesentlich allgemeineres Variationsprinzip, das Hamiltonsche Prinzip anwenden 
kann. Auch beim Grammelschen, Verfahren laBt sich fragen, wieweit es nur ein besonderes mathe- 
matisches Verfahren zur Lésung bestimmter Eigenwertprobleme ist, und wieweit es sich dartiber 
hinaus auf ein allgemeines Grundgesetz der Mechanik griindet. Diese Frage beantwortet sich 
nicht von selbst, weder aus der Herleitung dieses Naherungsverfahrens iiber eine quellenmaBige 
Darstellung, noch aus seinem Zusammenhang mit den Schwarzschen Quotienten*. Damit ist 
natirlich nicht gesagt, daB man die beim Grammelschen Verfahren vorkommenden Ausdriicke 
bei bestimmten Problemen nicht etwa mechanisch deuten kénne. 

Bei Schwingungsaufgaben ist es nun méglich, das mechanische Grundgesetz, das auch hinter 
dem Grammelschen Verfahren steht, aufzuzeigen. Es ist dies wieder das Hamiltonsche Prinzip. 
Wahrend dieses aber, auf Eigenschwingungsaufgaben unmittelbar angewendet, bei eingliedrigem 
Ritz-Ansatz immer auf den Rayleighschen Quotienten fiihrt, muf man den mechanischen Schwin- 
ger erst einmal umdeuten, ehe es jenen andern Quotienten des Grammelschen Verfahrens liefert, 
den wir im folgenden der Kiirze halber den Grammelschen Quotienten nennen wollen. Diese 
Umdeutung besteht darin, daB man dem Schwinger ein widerstandsreziprokes, ebenfalls schwin- 
gungsfahiges Gebilde mit denselben Eigenfrequenzen zuordnet®. Soweit man diese Methode der 
Schwingertransformation beherrscht, kann men also auch die in vielen technischen Anwendungen 
besonders vorteilhaften Grammelschen Gleichungen aufstellen, baw. lat sich damit in Zweifels- 
fallen nachpriifen, ob die auf iiblichem Wege gewonnenen Ansatze in sich folgerichtig sind. 

Im folgenden beschranken wir uns, was unwesentlich ist, auf die Herleitung des Grammelschen 
Quotienten bei eindimensionalen, einlaéufigen® Schwingern. Diese Voraussetzungen sind 
wesentlich, da sie durch die beabsichtigte Schwingertransformation bedingt sind. Mehrlaufige 
Schwinger unterscheiden sich von einlaufigen dadurch, daf sich ihre potentielle Energie nicht 
mehr mit Hilfe des Produkts aus einer Schnittkraft samt der ihr zugehérigen Verformung 
darstellen 1aBt. In diesem Sinne ist der Biegestab im Fliehkraftfeld, aber auch die radial schwin- 
gende Scheibe nicht mehr einlaufig. Beide sind kinematisch noch eindimensionale Schwinger. 
Soweit der in solchen Fallen zusatzlich auftretende Ausdruck in der potentiellen Energie aber nur 
eine Korrektur darstellt — und darum handelt es sich z. B. beim Einflu8 der Fliehkraft auf die 
Eigenfrequenzen eines Biegestabes —, kann man den mehrlaufigen Schwinger naherungsweise 
immer als Uberlagerung mehrerer einlaufiger Schwinger ansehen. Auch wenn man damit auf 
eine Fehlerabschatzung verzichtet, so hat sich dieses Vorgelen, das auf dem Satz von Southwell 
beruht, bei technischen Problemen durchaus bewahrt. 

Wir werden es vorziehen, den skizzierten Gedankengang nicht in aller Allgemeinheit an einem 
abstrakten eindimensionalen einlaufigen Schwinger, sondern anschaulicher an den Beispielen der 


1 Vgl. z. B. L. Collatz, Eigenwertaufgaben mit technischen Anwendungen, Leipzig 1949. 

2 E. Stiefel und H. Ziegler, Z. angew. Math. Physik 1 (1950) S. 111. 

3 R. Grammel, Ing.-Arch. 10 (1939) S. 35. — Vgl. auch C. B. Biezeno und R. Grammel, Technische Dynamik, 
1. Bd. S. 177, 2. Aufl. Berlin/Gottingen/Heidelberg 1953. 

4 L. Collatz, a.a.O. FuBn.1, S$. 203. 

> K. Klotter, Ing.-Arch. 18 (1950) S. 291, insbes. 5S. 299. 

6 Zu diesem Begriff vgl. R.Grammel, Uber die Lésung technischer Eigenwertprobleme, Forsch.-Hefte 
Stahlbau, Heft 6, S. 36, Berlin 1943. 
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Torsions- und Biegeschwingungen von Staben naher auszufiihren. Auf diese Weise lassen sich die 
Besonderheiten, die durch die verschiedenen Randbedingungen auftreten kénnen, in ihrer Viel- 
falt besser darstellen. Dabei werden sich auch einige Fingerzeige dafiir ergeben, unter welchen 
Bedingungen das Grammelsche Verfahren besonders zu empfehlen ist, und wann der Rayleighsche 
Quotient den Vorzug verdient. 


2. Torsionsschwingungen von Staben. a) Stetig verteilte Drehmasse (Abb. 1). Mit der 
Dichte 9, dem polaren Flachentragheitsmoment I(x) und der Torsionssteifigkeit C(x) geniigen 
die freien Torsionsausschlige p(x, t) und das dabei geweckte Torsionsmoment M(x, t) auBer 

gewissen Randbedingungen den Gleichungen 


ey  oM 
a = Ol(x)ae = 3, (1) 
0 J 
Abb. 1. Torsionsstal M = C(x. (2) 
ID, 1. orsionsstab. Ba 


Gleichung (1) ist der Drehimpulssatz, Gleichung (2) eine rein statische Beziehung, das Elastizi- 

tatsgesetz des Stabes. Ein solcher Torsionsstab kann insbesondere stationare, freie Schwingungen 
w(x, t) = y(x) sinot, (3) 
M(«, t) = N(x) sino t 

ausfihren. 


Auf das zu dem gegebenen Torsionsstab reziproke Schwingungsgebilde kommen wir, wenn 
wir als neue Veranderliche 


pra Mo = Meee? (4) 


wahlen. Damit schreiben sich die Gleichungen (2) und (1) um in 


Py 
(Cle ae (5) 
ow* 
M* — (Ge aa 9 (6) 
wenn man noch 
fits eee 


setzt. Die so erhaltenen Gleichungen sind formal gleich gebaut wie die urspriinglichen Glei- 
chungen, kénnen also ebenfalls als Schwingungsgleichungen eines Torsionsschwingers angesehen 
werden, wenn auch die physikalischen Dimensionen andere geworden sind. So ist beim reziproken 
Schwinger das Maf fir seine ,,Massentragheit’* durch die reziproke Torsionssteifigkeit des ur- 
spriinglichen Schwingers und ebenso seine ,,Torsionssteifigkeit“‘ durch die reziproke spezifische 
Drehmasse gegeben. Diese Unstimmigkeit in den physikalischen Dimensionen ist aber unwesent- 
lich, da man ja nur die Ausgangsgleichungen in dimensionsloser Form hatte zu schreiben brau- 
chen, um dieser Schwierigkeit zu entgehen. So wie also (2) das statische Grundgesetz des wirk- 
lichen Schwingers und (1) seine aus dem fiir ihn giiltigen Hamiltonschen Prinzip folgende Eulersche 
Gleichung ist, kénnen wir diese Betrachtungsweise auch auf die Gleichungen (5) und (6) des rezi- 
proken Gebildes anwenden, und auferdem ist klar, daB insbesondere das reziproke Gebilde 
dieselben Eigenfrequenzen o wie der gegebene Schwinger hat, da es ja tatsachlich denselben 
Schwinger darstellt. Wir kénnen somit, statt von den Differentialgleichungen auszugehen, diese 
Kigenfrequenzen auch mit Hilfe des Hamiltonschen Prinzips, angewendet auf das reziproke 
Gebilde, ermitteln, sofern wir uns nur zuvor noch davon iiberzeugen, was bei der Transformation 
aus den Randbedingungen geworden ist. 


b) Randbedingungen einschlieBlich Einzelmassen. Als einfachste Randbedingung 
betrachten wir zundchst die kinematische der starren Einspannung an einer Stelle x = a 
waht OS WS Ihe 


p(a,t)=0. (8) 
Die dazu reziproke Bedingung lautet nach der ersten Gleichung (4) 


Ea, i)=0, (9) 
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ae das ist die dynamische Bedingung eines freien Stabendes. Umgekehrt sieht man sofort 
af} die dynamische Bedingung eines freien Stabendes 


9 


. cae M(a, t) = 0 (10) 
in die kinematische Bedingung der starren Einspannung 
p*(a,t)=0 (11) 


transformiert wird. 
Sei i ree , ; 
Re En aoa Stab an einer Stelle elastisch eingespannt, wobei es geniigen mag, 
: etrachtung fir das linke Stabende x = 0 durchzufiihren. Dann wird bei Auslenkungen 
(| a) ein riickdrehendes Lagerreaktionsmoment geweckt, und die entsprechende Randbedingung 
wird mit einer positiven Federungszahl ¢) vorzeichenrichtig durch 


M(0, t) = cy p(0, 2) (12) 


angegeben. Zweimalige Differentiation nach der Zeit t liefert die entsprechende reziproke Rand- 
bedingung 
Oy 1 
fa) ** . F 
oO; (at the Sa ora (13) 
x=0 Co 
ee diese besagt, daB an diesem linken Stabende nunmehr eine Scheibe mit endlicher Drehmasse 
Oy = l/c, angebracht zu denken ist. Umgekehrt zeigt man, daf} die Bedingung fiir eine endliche 
Drehmasse an einem linken Stabende 


py es 
o(qr),_ = + lO = 
in die dazu reziproke Bedingung elastischer Einspannung 
DM *(O; t) == cP wy" (0; 4) mit Cos 5 (15) 
0 


ubergefihrt wird. Es werden also auch alle Randbedingungen beim reziproken Gebilde wieder 
in mechanisch sinnvolle und im Bereich der Torsion verbleibende Bedingungen umgewandelt. 
Die Reziprozitat ist vollstandig durchfiihrbar; somit sind die Eigenfrequenzen auch am rezi- 
proken Gebilde und zwar mit Hilfe derselben Satze wie beim gegebenen Schwinger zu ermitteln. 

c) Der Rayleighsche Quotient. Wendet man das Hamiltonsche Prinzip auf stationare 
EKigenschwingungen an, so kommt man zu der Variationsaufgabe+ 


6(U—/AK)=0 mit fpeMor (16) 


wobei U die zur Amplitudenlage des Schwingers gehérige maximale potentielle, A K seine maxi- 
male Bewegungsenergie beim Durchgang durch die Null- oder Gleichgewichtslage ist. Auf 
Grund ihrer mechanischen Bedeutung kann man die Ausdriicke fiir U und K allein aus den 
statischen Gleichungen und den Randbedingungen bilden und dann, entsprechend dem gewahlten 
Lésungsverfahren, die Variationsaufgabe (16) weiter behandeln. Alles fiir das folgende Wesent- 
liche ist zu ersehen, wenn wir die Eigenfrequenz naherungsweise mit Hilfe eines eingliedrigen 
Ritz-Ansatzes M(x) fiir die Schwingungsaiplitude bestimmen. Man gelangt so in bekannter 
Weise zum Rayleighschen Quotienten, der im allgemeinen einen Naherungswert fiir das Quadrat 


der ersten Eigenfrequenz liefert: 
Ue) ie 


Dabei braucht die Vergleichsfunktion @ nur die geometrischen Randbedingungen fur die Schwin- 
gungsamplitude zu erfiillen. Man kann also stets solche Vergleichsfunktionen angeben, da in 
diesen Randbedingungen der Eigenwert A nicht auftritt. 

d) Der Grammelsche Quotient. Ubertragen wir diesen Sachverhalt auf die reziproke 


Schwingungsaufgabe, so erhalten wir den Rayleighschen Quotienten 


eee | 18 

Ay = K*[@*| : ( ) 
Die Vergleichsfunktion @* mu wieder nur etwaige kinematische Randbedingungen jetzt des 
reziproken Gebildes erfillen; bei der Torsion ware das lediglich die (dynamische) Bedingung fiir 


1 Vgl. C. B. Biezeno und R. Grammel a. a. OMBupn ec. soskid tf 
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ein freies Stabende beim urspriinglichen Schwinger. Explizit lautet dieser Quotient, wenn wir 
im Beispiel auf endliche Drehmassen und elastische Lagerungen fiirs erste verzichten, 


| I 
*  /d@*? 1 
4) ee <= *2 * 
[o(a) ‘ae [O*] dx 
Hh oe 0 * 6 aes 0 : oe ; (19) 


[ I)* G* dx [e 1)* ®* dx 


0 0 


und es besteht zwischen N* und ®* nach (6) und den transformierten Gleichungen (3) die Be- 
ziehung 
abDp* 


peer 
IN eG qe 


(20) 


Wir wollen nun zur Verbesserung des Ergebnisses eine Iteration vorschalten, indem wir nicht 
von einer zuldssigen Vergleichsfunktion ®* ausgehen, sondern umgekehrt die Momenten 
amplitude N* wahlen und dazu nach (20) unter Berticksichtigung der Randbedingungen fir @ 
diese Funktion erst bestimmen. Die Ausgangsfunktion N* kann an sich beliebig gewahlt werden, 
da das Hamiltonsche Prinzip nur der kinematischen, nicht der dynamischen Koordinate Bedin- 
gungen auferlegt. Wegen der nach der zweiten Gleichung (4) bestehenden Reziprozitat, bei der 
es auf einen konstanten Faktor im Ritz-Ansatz nicht ankommt, wollen wir aber 


NG aD (21) 
wahlen, wo ® wenigstens die kinematischen Randbedingungen des urspriinglichen Schwingers 
erfiillen soll. Damit erhalten wir aus (19) mit dem zweiten Ausdruck der rechten Seite und mit 
(7) und der ersten Gleichung (4) 
vg I 

[ 0 1G dx 

0 
Nae hme arcs (22) 
ie [MGs el adv 


Abb. 2. Biegestab. 0 


Sh] 


wo also N(x; 0 [ ®) = @* das aus (20) fiir N* = @ durch eine Quadratur ermittelte maximale 
Torsionsmoment ist. Bei dieser Integration sind natiirlich die dynamischen Randbedingungen 
des Schwingers beriicksichtigt. 


Der Quotient (22) ist nun gerade der vom Grammelschen Verfahren bei eingliedrigem Ansatz 
fiir dieses Beispiel gelieferte Quotient einschlieBlich des Iterationsschrittes, der sich beim Gram- 
melschen Verfahren so vorteilhaft auf die Giite der Naherung auswirkt. Der Grammelsche 
Quotient erweist sich daher einfach als Rayleighscher Quotient fiir das reziproke Gebilde, wenn 
man dort die Vergleichsfunktion nicht unmittelbar vorgibt, sondern nach der obigen Vorschrift 
aus der zu wahlenden dynamischen Gréfe erst konstruiert. Entsprechend lassen sich die 
Grammelschen Gleichungen mit Hilfe mehrgliedriger Ansatze herleiten, und weiter laBt sich das 
Verfahren wegen der mechanischen Bedeutung der vorkommenden Ausdriicke U* und K* fiir 
jeden beliebig gestalteten Torsionsschwinger im Prinzip ansetzen. Welche Erschwerungen in- 
dessen bei seiner Durchfiihrung auftreten kénnen, werden wir in Ziff. 4 naher untersuchen. 
Zunachst mbge der Gedankengang noch am Beispiel der Biegeschwingungen von Staben er- 
lautert werden; diese zahlen auch zu den eindimensionalen einlaufigen Schwingern. 


3. Biegeschwingungen von schlanken Stiben. a) Stetig verteilte Stabmasse (Abb. 2). 
Sei F(x) der Stabquerschnitt, a(x) die Biegesteifigkeit, M(x, t) das zur Biegelinie 7(x, t) gehérige 
Biegemoment, so haben wir jetzt anstelle von (1) und (2) die Gleichungen 


td A Ree RC (23) 


BIC) ere (24) 
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Fir die stationadren Kigenschwingungen gilt dabei 


n(x, t) = y(x)sinot, 
M(x, t) = N(x) sinot. (25) 


Auf die ents he i : . 
Bier prechenden Gleichungen des reziproken Gebildes kommen wir durch die Trans- 


ee 077 
Pi VE. M* = — = : (26) 
wodurch die Gleichungen (24) und (23) in 
+ On eM* 
(Cla sae t (27) 
0*n* 
Mo (28) 
mit 
ed, g aed 
(0 F)* = ae i oF (29) 


ubergefiihrt werden. Das sind wieder Gleichungen von der Form (23), (24). 
b) RandbedingungeneinschlieBlich Einzelmassen. Aus einer freien Auflagerung 


f= 05 M=0 (30) 
wird beim reziproken Gebilde 
M*t =0, jie S208 (31) 


das ist wieder eine freie Auflagerung. Eine starre Einspannung 


a ae Oy 
=0 « = (32) 


geht in ein freies Stabende (und umgekehrt) 


aM* 
M* =0, pone (33) 


iiber. Beides sind Sonderfalle einer elastischen Einspannung, die mit positiven Feder- 
steifigkeiten p und q bei einem linken Stabende x = 0 


Bg Li 
ee Ox’ 
(34) 
Ox a q / | 
verlangen. Die dazu reziproken Bedingungen lauten 
m* eur = oM™ mit i ss a é | 
ot? Ox Pp (35) 
a2 On f eg 1 
Ne ee * + — — 
D AP ax M mit D ra 


d. h. beim reziproken Gebilde haben wir an diesem linken Stabende eine starre Scheibe von der 
Masse m* und der dquatorialen Drehmasse D* anzusetzen. Umgekehrt verwandelt sich eine 
solche Scheibe (m, D) reziprok in ein elastisches Lager mit den Federsteifigkeiten p* = 1/m, 
q* = 1/D, und es werden linke Stabenden in linke und rechte Stabenden in rechte tibergefihrt. 
Auch hier ist also die Reziprozitat vollstandig. 


c) Der Grammelsche Quotient. Nunmehr lat sich ohne weiteres der Rayleighsche 
Quotient fiir das reziproke Gebilde ansetzen. Sehen wir vorlaufig von Einzelmassen und elasti- 
schen Lagern ab, so lautet er mit irgendeiner zulassigen Vergleichsfunktion Y* fiir die Biege- 
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amplitude ¥ 
| 


| 1 I [Y*] dx 
“/ oe 
0 ; (36) 


ee ee 
(0 F)* Y* dx 
0 


Schalten wir auch hier wieder einen Iterationsschritt vor, indem wir nicht Y* wahlen, sondern, 


ausgehend von 


mit Hilfe der statischen Biegegleichung (24), also jetzt 
APN Ge 1 7 
Ge =a =eFY, (374) 


die Vergleichsfunktion Y* = N(x;o0 F Y) durch Quadratur ermitteln, so wird aus (36), auf die 


urspriinglichen Veradnderlichen umgeschrieben, 
| 
Ol Yds 
0) 4 (38) 


A= _» 


l 
S IN2(x3 oF Y))dx 
) 


und das ist der Grammelsche Quotient. Die Funktion N(«; 0 F Y) laBt sich als ,,Biegemoment 
bei stetiger ,,Biegebelastung’’ 0 F Y deuten. 

d) Masseloser Stab mit Einzelmassen. Es ist 
lehrreich, das Verfahren auf den Grenzfall einer frei 
gelagerten, biegeelastischen, masselosen Welle anzu- 
Ly wenden, die an n Stellen x, starre Scheiben (m,, D,) 
tragt (Abb.3). Das ihr zugeordnete reziproke Gebilde 
besteht aus n-+-1 starren Teilstaben von der spe- 
zifischen Masse (9 F')* = 1/a, die an den Stellen a; 


GG 


8 


0 | Z ee : 
| -; derart federnd miteinander verbunden sind, daf} dort 
| Querkrafte Qj und Biegemomente Nj gemafs 

x I 
ae — p, (P* (a. 0) "x, 0), 
Vi : (39) 
Abb. 3. Masselose Welle mit Scheiben und reziprokes J ees —(Y* (% 0) — Y* ; aap 
: Gebilde Fel ee ae : cp my ( (% 3 ) = (a 0)) 


geweckt werden. Gehen wir also nach (36) von einer Funktion N* = Y(x) aus, die die geometri- 
schen Randbedingungen Y(0) = Y(I) = 0 befriedigt, so haben wir dazu einen an den Stellen x, 
gebrochenen Polygonzug Y*(x) aus (39) mit Ni = Y, = Y(«,), Qf = Yi = Y’(x,) so zu er- 
mitteln, dal fir ihn Y*(0) = Y*(/) = 0 wird; oder wenn wir wieder die nichttransformierten 
Veranderlichen bevorzugen, so haben wir zu einer auBeren »,Belastung’ aus freien Momenten 
mp Y;, und Querkraften D, Y;, an den Stellen x;, das Biegemoment N(x; mz Y;, Dy Yj), etwa iiber 
ein Seileck, zu konstruieren, und damit erhalten wir den Grammelschen Quotienten 
2 (m Yi + DY) 
Ay = eS a Sere (40) 
i 
a NN? (%3 my, Y), Dip). dx 
0 
Man iberzeugt sich leicht davon, daB auch der Weg ii ai 
; y ‘g uber die quellenmaBige Darstell {' 
denselben Quotienten fiihrt. Nur beilaufig sei bemerkt, daB der Rayleighsche panda are 
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wertaufgaben bei Matrizen!, wenigstens soweit es sich um ein Problem mechanischer Schwin- 
gungen handelt, nach der mechanischen Bedeutung seines Zahlers und Nenners auch als Grammel- 
scher Quotient angesprochen werden kénnte. Dies erklart sich daraus, daB in den Ausgangs- 
gleichungen solcher Probleme eine quellenmaBige Darstellung vorliegt. 


4. Besonderheiten. Es ware nicht schwer, die an den Beispielen der Torsions- und Biege- 
schwingungen von Staben durchgefiihrten Uberlegungen auf den einlaufigen, eindimensionalen 
Schwinger zu verallgemeinern. Eindimensional heift hier, daB die Bewegung des schwingungs- 
fahigen Kontinuums durch eine Gleichung der Form 


z= 2(x, t) (a <7 b) (41) 


beschrieben wird. Der Begriff einlaufig beschrankt die Schwinger auf solche Kontinua, bei 
denen die der Verformung 2(x, t) entsprechende innere Beanspruchung derart ist, da sich die 
spezifische potentielle Energie durch das halbe Produkt aus der Verschiebungskoordinate z und 
einer entsprechenden ,,Schnittkraft‘‘ S(x, t) darstellen laBt, wobei zwischen S und z noch eine 
rein statische Beziehung nach Art von (2), (24) oder auch allgemeiner besteht. Eine solche 
Verallgemeinerung wiirde aber auf eine bloBe Wiederholung unserer Uberlegungen hinauslaufen. 
Wir wollen daher darauf verzichten und dafiir die einzelnen Schritte unseres Vorgehens auf ihre 
Durchfithrbarkeit nachpriifen. 

Die Betrachtung hat sich auf eindimensionale, einlaufige Schwinger beschrankt, weil sich 
hier das reziproke Schwingungsgebilde ohne weiteres angeben lat. Wie die Reziprozitat bei 
mehrlaufigen Systemen durchzufiihren ist, ob eine solche uberhaupt aufgestellt werden kann, 
mu hier unerértert bleiben. Allererste Ansatze scheinen anzudeuten, daB die Struktur des 
reziproken Gebildes sich von der des gegebenen Schwingers starker unterscheidet, vielleicht in 
der Richtung, da Dimension und Laufigkeit des Schwingers einander reziprok gegentiberstehen. 

Verfolgt man nun das Verfahren weiter bei einlaufigen Schwingern, so bemerkt man alsbald, 
dafS§ zwei Umstande seinen Gang entscheidend beeinflussen. In das Verfahren ist ein Iterations- 
schritt eingebaut, der fiir die Genauigkeit des Verfahrens wesentlich ist; in unseren Beispielen 
handelt es sich dabei um die Quadratur der Gleichungen (20) mit N* (21) baw. (37). Diese 
Quadratur gelingt aber ohne weiteres nur dann, wenn der Schwinger statisch bestimmt gelagert 
ist. Bei unbestimmter Lagerung dagegen mite man die Integration bis zur iterierten Schwin- 
gungsamplitude weitertreiben, und dann ist es so, da der Grammelsche Quotient (allgemeiner: 
die Grammelschen Gleichungen) keine Vorteile mehr bietet, da man jetzt bei gleichem Rechen- 
aufwand den Rayleighschen Quotienten fiir die iterierte Schwingungsamplitude ansetzen kann 
und damit den Eigenwert im allgemeinen noch genauer erhalt. 

Weiter hat man zu beachten, ob beim reziproken Gebilde neben stetig verteilter Masse auch 
Einzelmassen auftreten. Dann findet sich der Eigenwert in den Randbedingungen des trans- 
formierten Problems vor. Fiir die unmittelbare Anwendung des Hamiltonschen Prinzips (16) 
bedeutet das kein Hindernis, da die Vergleichsfunktionen nur die kinematischen Randbedingungen 
zu erfiillen brauchen. Die gewiinschte Genauigkeit der Eigenwertrechnung wird dann allerdings 
mehrgliedrige Ritz-Ansiatze erforderlich machen. Beim Ubergang zum reziproken Gebilde wird 
dagegen infolge des Iterationsschrittes der Eigenwert/ unter Umstanden nicht nur in A K* 
sondern auch in U* auftreten und dadurch wird die Auflésung der Variationsgleichung 

6(U* -—/A K*) = 0 (42) 
nicht auf einen Quotienten fiihren, dessen Zahler und Nenner die seitherige Bedeutung haben. 
Dies bildet aber kein grundsatzliches Hindernis fiir die Durchfiihrung des Verfahrens. Das Ausmafs 
der Rechenarbeit bei der Iteration hangt einmal davon ab, ob der Schwinger statisch bestimmt 
gelagert ist oder nicht, und weiter von der Art der Randbedingungen. 

Als Beispiel an dem sich die verschiedenen Méglichkeiten schon ausreichend aufzeigen lassen, 
betrachten wir einen Torsionstab, der in x = 0 elastisch eingespannt ist und in x = / eine starre 
Drehmasse Q, tragt (Abb. 4). Der reziproke Torsionsstab mit 


1 : 1 ; 
(ON ace Cae; (43) 
besitzt dann an der Stelle x = 0 eine starre Drehmasse 
ot=+ (44) 
C9 


2 vel R. Zurmiihl, Matrizen, S. 155, Berlin/Gottingen/Heidelberg 1950. 
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und ist in « = / elastisch gelagert mit 


iL 

of eee 45 
Os (45) 

Damit wird 

1 
2Ut = ee [D*] dx + 4 NP [O*] (46) 
0 
und 
i 

2 K* = f{ (9 1)* O* dx + Of OF, (47) 


0 
wobei fiir N* und @* die Gleichung (20) und die Randbedingungen 
Ni =—A0%63, Nt=—c oF (48) 


gelten. 


Fur die Iteration gehen wir aus vom Ansatz 
. N* = Ox), (49) 


wobei die Vergleichsfunktion @ an keinerlei Randbedingungen gebunden ist, da bei dem gewahlten 
Schwinger (im allgemeinen) keine kinematischen Randbedingungen vorgeschrieben sind. Der 

Vergleich von (49) mit (48) zeigt aber sogleich, da8 sich diese 

Iteration bei ganz beliebiger Vergleichsfunktion ® nicht 

durchfiihren 1a8t, denn bei der Integration von (20) ist nur 
oLC £ eine Integrationskonstante frei, wahrend die Gleichungen 
(48) zwei verschiedene Bedingungen stellen, Darin drickt 
sich einfach die Tatsache aus, dafs das Problem statisch 
unbestimmt ist. Man hilft sich hier so, daB man tber @* 
hinaus die Iteration weiter bis zur kinematischen Veradnder- 
lichen des Schwingers durchfiihrt, indem man mit N = @* 
die Iterierte ® von @® aus 


Co 


Ce 


Abb. 4, Elastisch gelagerter Drehschwinger mit ma 
einer Zusatzdrehmasse und reziprokes Gebilde. d® 1 Ne 1 @* (50) 


dx C Cc 


bestimmt. Wie schon bemerkt wurde, ist dann aber der mit Hilfe dieser Iterierten D angesetzte 
Rayleighsche Quotient 


iy 2 SSS (51) 


dem Grammelschen Quotient vorzuziehen. 

Indessen bietet sich noch ein anderer Ausweg an, den wir noch verfolgen wollen, auch wenn 
iiber seine numerische Brauchbarkeit keine Erfahrungen vorliegen. Wir verzichten darauf, die 
Vergleichsfunktion ® ganz beliebig zu wahlen und untersuchen, welcher Einschrankung sie unter- 
liegt, damit das durch (20) und (48) mit N* (49) gestellte Problem lésbar wird. Mit einem noch 
unbekannten Anfangswert * = O*(0) folgt aus (20) 


x 


D*(x) = OF + | Dect (52) 


(Gee 
0 
Damit kommt aus (48) 


l 


DO) =—AOPO;, Pl) — ah Dy i & as) ; (53) 
0 


Aus der zweiten Gleichung (53) bestimmt sich ®* zum Wert 


l l 
x Dil i 
N(0) = @}= — 5 — [dx = 0,0) — [ol bade, (54) 
i) 6 
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wahrend die erste Gleichung (53) die Ausgangsfunktionen @(x) durch die Forderung 


OD eo 
oa * ; A P 
P(0) =1 0; ( a | a = ay [6 P(1) + le 1® ts) (55) 


0 


0 
einschrankt. 

Diese Forderung stellt eine Beziehung zwischen dem maximalen Drehimpuls des ganzen 
Schwingers und dem. maximalen Kinspannmoment her. Man kénnte sie durch eine einfache 
Integration der Drehimpulsgleichung (1) tber die Stablange unter Beachtung von (3) und der 
Randbedingung (14) unmittelbar herleiten. Wir genugen ihr durch den Ansatz 

C= Pi ua, Fs (56) 


mit willkiirlichen Funktionen Y,(«x), Y(«), wenn wir mit noch unbekanntem Kigenwert A 


f l 
Ma = ¢ Y(0) —A | 0, YD) + Jol, dx|=a,+/b,, 
0 
eed (51) 
[bo = — Cy ¥,(0) +A 6, v0) +] One, ix] =a,+/1b, 
0 
wahlen, so daB wir also mit den durch Y%, Y, eindeutig definierten Funktionen 
: ®,(x) =a,¥%, + a, ¥, D(a) = 6, Vo ee (58) 
au 
D=O, +10, (59) 
kommen. Dadurch ist nach (50), (52) mit 
I 
N(x) = — fe 1D; dx — 0; D,(l) (i = 1, 2) (60) 
auch die Iterierte : 
N(x; 1@) = @*(x) = Ni(x) + A Nga) (61) 


gefunden. 
Die fiir die Variationsaufgabe (42) benétigten Ausdriicke U*, K* sind jetzt von der Gestalt 


I 
2U* = [oI @ dx + 6, G1) = Ay, + AA + Appl? , | 


0 
“4 1 ed 
2K* = ao (x39 1 ®) dx + —-N® (030 1D) = By + Bul +By??. 
0 ’ 
Dabei bestimmen sich die Konstanten 4;;,, B;, eindeutig aus den Funktionen ®,, ®, und zwar 
geben die Zeiger an, welche dieser Funktionen in die Rechnung eingehen. Wir betrachten einige 
Sonderfalle. 

a) Einseitig starr eingespannter Stab mit freiem Ende (q = ~, 0, = 0). Glei- 
chung (55) vereinfacht sich auf die wohlbekannte kinematische Randbedingung ®(0) = 0. Wir 
kénnen @,(x) = 0 setzen. Die Aufgabe ist auferdem statisch bestimmt. Der eingliedrige Ansatz 
fiihrt ohne irgendwelche Schwierigkeiten zum Grammelschen Quotienten. 

b) Beidseitig starr eingespannter Stab (¢ =o, @;=o0). Nunmehr ist (55) nach (53) 
mit den beiden kinematischen Randbedingungen ®(0) = Oil) =0 gleichbedeutend. Aber 
jetzt ist ®* aus (52) nicht mehr zu berechnen, was nur ein anderer Ausdruck dafir ist, daf die 
Aufgabe statisch unbestimmt ist, und zwar in einer Weise, daB das allgemein geschilderte Ver- 


fahren undurchfiihrbar ist. Man ist gendtigt, die Iterierte PD mit D(0) = Dl) = 0 nach (50) zu 
bestimmen und wird dann den Rayleighschen Quotienten (51) auswerten. 
c) Stab mit freien Enden (c, = 0, 0; = 0). Damit die einfache Iteration gelingt, mu 


nach (55) 
l 
[oT @dx =0 (62) 
0 


sein. In dem Ansatz fir ® tritt der Kigenwert nicht auf. Aus (54) folgt O| = 0, Das Problem 
ist statisch bestimmt, daher wird man den Grammelschen Quotienten bilden. 
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d) Stab mit einem freien Ende und einer Drehmasse am andern Ende 


(co= 0, O, +0). Die Bedingung (55) lautet hier 
U 

0, O(l) + |e I@dx=0. (63) 
0 


Aus (54) folgt eindeutig @* = 0; die Aufgabe ist statisch bestimmt. Der Grammelsche Quotient 
1aBt sich sofort bilden. 

e) Starr eingespannter Stab mit endlicher Drehmasse (¢ = ©, QO, + 0). Fur P(x) 
besteht nach (55) nur die kinematische Randbedingung ®(0) = 0. Da ®} (54) nicht vom Eigen- 
wert A abhangt, gilt dies auch von der Iterierten N (61). Wiederum kann man den Grammelschen 
Quotienten anschreiben. Der Eigenwert wiirde sich hier in den Randbedingungen erst storend 
bemerkbar machen, wenn man auch die Iterierte ® von ® berechnen wollte, da fiir @ neben 


P(0) = 0 die Randbedingung 


—1A0,G(l) = N(l) (64) 
gilt. Wir haben hier also einen Fall, bei dem es zwar noch den Grammelschen, aber nicht mehr 
den Rayleighschen Quotienten (51) gibt. 

f) Stab ameinen Ende elastisch eingespannt, am andern Ende frei, (¢ 4-08 
QO, = 0). Es liegt der zu Fall e) reziproke Fall vor. Am zweckmafigsten geht man also erst zum 
reziproken Gebildé iiber, nimmt es als gegebenen Schwinger und verfahrt nach Fall e). Man 
kénnte aber auch die Vergleichsfunktion ® den Randbedingungen 


¢ ®(0) = C(0) (0), ®(l) = 0 (65) 


fiir den urspriinglichen Schwinger anpassen und dann den Rayleighschen Quotienten U[®]/K[®] 
bilden. Dagegen empfiehlt es sich weniger, das Grammelsche Verfahren unmittelbar auf den ur- 
springlichen Schwinger anzuwenden. Dazu miBte die Ausgangsfunktion der Nebenbedingung 
l 
(0) =+ le 1 ® dx (66) 
“e 
geniigen. Dann hangt aber die Iterierte N (61) vom Eigenwert 4 ab, so daB man mit einem ein- 
gliedrigen Ritz-Ansatz zu (42) eine Gleichung héheren Grades in 2 bekommt, die nicht durch 
einen Grammelschen Quotienten gelést wird. Auch darin erweist sich dieser Fall als ein Gegenstiick 
zu Fall e), bei dem diese Erschwerung umgekehrt am Rayleighschen Quotienten haftet. 

g) Stab am einen Ende elastisch gelagert, am andern Ende starr eingespannt 
(cy + 0,0, = oo). Es liegen ahnliche Verhaltnisse wie bei Fall f) vor. Um das Grammelsche 
Verfahren anzuwenden, beniitzt man am zweckmaBigsten, da Fall g) reziprok zu Fall d) ist. 

h) Allgemeiner Fall nach Abb. 4 (c« +0, 0; +0, endlich). Jetzt tritt der Eigenwert 
bei jedem Vorgehen in den Randbedingungen auf. Bei eingliedrigem Ritz-Ansatz laBt sich weder 
ein numerisch brauchbarer Rayleighscher noch der Grammelsche Quotient bilden. Wie man hier 
mit einem solchen Ansatz noch zum Ziele kommen kénnte, ist den Gleichungen (52) bis (62) zu 
entnehmen. In der Regel wird man aber die Variationsaufgabe (16) durch hinreichend vielgliedrige 
Ritz-Ansatze bequemer lésen kénnen. Die Koordinatenfunktionen dieser Ansatze brauchen 
jetzt keinen Randbedingungen mehr zu geniigen, da keine allein zu beachtenden kinematischen 
Bedingungen mehr vorhanden sind. 

Fir ein Zahlenbeispiel, bei dem man die Ergebnisse mit der strengen Lésung vergleichen 
kann, wahlen wir einen Torsionsstab mit konstantem Querschnitt (o I = konst.. C —ikonaee 
Die Zusatzdrehmasse sei 0, =o I1, d.h. gleich der Drehmasse des Stabes, die Federungszahl 
bei eJastischer Lagerung sei cy = C/I, also gleich der spezifischen Torsionssteifigkeit des Stabes. 
Die Ergebnisse sind in der Tabelle 1 zusammengestellt. Diese enthalt fiir verschiedene Schwinger- 
anordnungen zunachst den exakten Bezugswert o, Jo I I2/C der niedersten Eigenfrequenz, sodann 
die Vergleichsfunktion ®, mit der die Naherungswerte samt ihrem Fehler bestimmt sind. Dabei 
bedeutet R den Naherungswert, der sich ergibt, wenn man mit @ in die Variationsaufgabe (16) 
eingeht, also in der Regel die Wurzel aus dem Rayleighschen Quotienten; der Naherungswert G 
ist mit dem Grammelschen Verfahren berechnet, bei R;, hat man erst zu © die Iterierte @ 
bestimmt und dann den Rayleighschen Quotienten gebildet, Das Beispiel zu Fall h), in dem 
der Eigenwert in den Randbedingungen auftritt, ist zum Vergleich auch mit mehrgliedrigen Ritz- 
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Ansatzen durchgerechnet. Bei 


als wenn man das in den Gleichungen (52) bis (62) 
empfiehlt sich alsoim Falle h) nicht. 
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hat. 


weit besserem Zahlenergebnis ist der Rechenaufwand geringer, 
geschilderte Verfahren anwendet. Dieses 
Tabelle 1 zeigt tibrigens erneut, daB das Grammelsche 
Verfahren da, wo es anwendbar ist, gegentiber dem Rayleighschen (Ritzschen) Verfahren an 
Genauigkeit stark tiberlegen ist, worauf schon Grammel! hingewiesen 


Tabelle 1. Niederste Eigenfrequenzen fiir Torsionsschwinger samt Naherungswerten 


a Saye | 
Poll a, | _ | P(x) | R CG Rie 
a Ea A | 1,7321 15811 . 
ie ) : _10,3% 0,06% ae antl 
3,1623 3,143 621 3,14 181 
p sah ed) 0,66% 0,065% 0,007% 
| ede 3,464] /10 
¢ es aga 10,3% 0,66% ny 
=A 2,19 089 2.04124 
d 2,028 4x—3l , ? om 
: | a i 8,0% 0,62% 
ieee 2% 
| 0,86603 0,86 066 
e 0,86033 we 0,66% 0,038% — 
ce 
jes (2 .— | 2,8319 0,58 183 
410% 4,65% 
=== a8 
| 0,86 603 
Peas 0:55 597 #3 55,8% = 
ce | -0,55:873 =) 
| ee a 0,497% 
i nets 0.55 598 
ee eae 0,0018% 


5. Zusammenfassung. Zu den Grammelschen Gleichungen fir die naherungsweise Bestimmung 
von Eigenwerten gelangt man auf zwei Wegen, entweder itiber die quellenmaBige Darstellung 
des Problems, d.h. iiber eine Integralgleichung oder iiber ein System von endlich vielen linearen 
Gleichungen, die man in bestimmter Weise umformt, — oder iiber die Differentialgleichung und 
ein besonderes Variationsprinzip, das an eine Kigenschaft der Schwarzschen Quotienten ankniipft. 
In der vorliegenden Arbeit ist ein dritter Zugang zu den Grammelschen Gleichungen angegeben, 
wenigstens fiir einlaufige eindimensionale Schwinger: Man fihrt das zum Schwinger widerstands- 
reziproke Gebilde ein und wendet darauf das Hamiltonsche Prinzip an. Damit ist gegeniber 
den ersten beiden Verfahrensweisen einerseits der Anwendungsbereich der Grammelschen Glei- 
chungen erweitert, da das Hamiltonsche Prinzip als Grundgesetz der Mechanik auf wesentlich 
weniger Voraussetzungen aufbaut als die mathematische Eigenwerttheorie (die dafiir allerdings 
dann auch genauere Aussagen iiber den Naherungscharakter der Lésung machen kann); anderer- 
seits wird die Betrachtung auf einlaufige, eindimensionale Schwinger eingeschrankt. Bei der 
Durchfithrung der Uberlegungen zeigt es sich, daB durch die Art der Randbedingungen ent- 
schieden wird, welches Verfahren bei der Lésung des Eigenwertproblems im Rahmen der tech- 
nischen Genauigkeitsanspriiche den Vorzug verdient. Es ist dies in den meisten Fallen das 
Grammelsche Verfahren. Soweit aber Erschwerungen auftreten, kann man stets mit mehr- 
gliedrigen Ritz-Ansatzen zum Ziele kommen, da Koordinatenfunktionen, wie sie vom Hamilton- 
schen Prinzip gefordert werden, stets angegeben werden kénnen. 


(Eingegangen am 21. Marz 1956.) 


Anschrift des Verfassers: Dozent Dr. K. Zoller, Stuttgart-Sillenbuch, Kirchheimerstr. 21. 
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The Motion of Automobiles in Unbanked Curves 


By P.R. Paslay and A. Slibar 


1. Introduction. The motion of an automobile in a curve has been treated both experi- 
mentally and theoretically by various authors; their analyses have shown the complexity of 
the problem and its many close relations to other problems in automotive engineering. For 
many problems connected with the 


Voss motion of an automobile in a curve the 
a a As) complete solutions are still unknown. 
wh Two such problems are the behavior 
ot of tires and the response of an auto- 
lace ! 6x9) mobile going over a bump. In spite of 
Rye the lack of information in so many 
ye a related fields, the treatment given 
0 


pole here shows that for driving through 
a | | 948 an unbanked curve theoretical results 
which check existing experimental re- 
sults can be obtained by an approxi- 
mate theory by making use of previous 
experimental work. 


2. Review of Previous Work. The 
mat best way to review work related to this 
problem is to first outline the results of 
experimental work which involves the 
behaviour of tires and then to outline 
work which involves the response of 
the vehicle itself. 


a) Tire Behaviour. In 1940 L. Hu- 
700 IL ber’ published results of a very exten- 
ded experimental project dealing with 
the characteristics of tires. The most 
important results of his treatment, as 


ee : : 
r regards the present work, is that in. 

order for a tire to be able to transmit 

| a force not acting in the plane of the 

| o wheel the plane of the wheel has to 

a 5 1 % degr be inclined at an angle to its forward 
Fig.1. The sideforce S as a function of the slip angle x of the wheel velocit ¥ This fact was pr 

a for different normal loads, from L. Huber [1]. et . a Ly ee Dy. se 

Die Seitenkraft S als eine Funktion des Schriiglaufwinkels ~ des Rades on models of tires of different sidewise 

fiir verschiedene Normallasten nach L. Huber [1]. stiffness even for cases of high coeffi- 


cient of friction. The angle at which 
the plane of the wheel is inclined to the forward velocity is called the Slip angle“ and is 
denoted by the letter a here (see Fig. 1). 


For fast driving in a curve it turns out from experiments as well as from present theory, that 
the sideforce from the ground to the tires necessary to overcome the centrifugal force requires 
slip angles up to 10°. Sidewise wind forces alone during straight ahead driving can create slip 
angles? up to 5°. Huber found that in general stiffer tires and higher tire pressures reduce the 
slip angle necessary to transmit a given sideforce. 


1 L, Huber, Die Fahrtrichtungsstabilitat des schnellfahrenden Kraft D 
deen Ee a enden aftwagens, Dtsch. Kraftfahrtforschung, 


2 O. Dietz und R. Harling, 


Heft 44, Berlin 1940, Die Fahrlage des Kraftwagens in der Kurve, Dtsch. Kraftfahrtforschung, 
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The sideforce transmitted by a tire at a given slip angle is not proportional to the load the 
tire carries. The sideforce increases less than proportional to the load on the tire in this case. 


Fig. 1, taken from L. Huber, shows the relation between the sideforce S and the slip angle « 
for various normal st 


loads on the tire. The 
reader should note 
that in Fig. 1 the 
acting sideforces are 
given with no depen- 
dence on the relative 
velocity between the 
tire and the ground. 
Pig. 2, showing the 
sideforce S as a func- 
tion of the slip angle 
a for a constant load 
and a wide range of 
relative velocities, 
points out that the 


60 


eqesoroor poco no -- = f--------- 


a 
sideforce is practi- 7% a a 6 8 10 12 14 degr 
call independent of Fig. 2. The sideforce S as a function of the slip Die Seitenkraftt S als Funktion des Schraglauf- 
yi P 
inc ° angle o for different values of relative velocity winkels x fiir verschiedene Werte der Relatiy- 
ve ocity. Huber’s ex- with constant normal load, from geschwindigkeit bei konstanter Normal!last nach 
periments further- L. Huber [1]. L. Huber [1]. 

h that th The parameters of relative velocity are: Die Parameter der Relativgeschwindigkeit sind: 
more show : a € 5 mph — dash dot curve, = 8,0 km/h — strich-punktierte Kurve, 
curves of Fig. 1 are ze a = — pe 30,0 mae — fest ausgezogene Kurve, 

= mph — dashed curve, = 60,0 km/h — strichlierte Kurve, 
not appreciably chan- 55 oe — dotted curve. =~ 90,0 km/h — punktierte Kurve. 
he : By Dhabas the Ubstt 4 Aligning Torque 
side stiffness of the tire. Asasecondary 
: j 700 at lia — 
result experiments of different authors, 5° 
for example *, have shown that a self Ue 
aligning torque is created by the non- 7 Ga 
uniform distribution of the pressure Vy 
acting between tire and ground. This 80 cae We g 
self-aligning torque may be represented We 
by a vector perpendicular to the plane : 
of the ground. The sense of the torque 
is to turn the plane of the tire toward 
the direction of the impressed relative 
velocity. Fig. 3 shows the dependence 
of the aligning torque on the normal 
load for different values of slip angles 
taken from M. Olley. 
Other authors like J. Bradley and 
R.F. Allen4, R. D. Evans®, A. W. Bull® 
and W. F. Billingsley, R. D. Evans, W. 
A. Hulsmit, and E. A. Roberts’? have 
published additional experimental re- 
sults dealing with tires. The most com- 
plete reviw in English literature of all 
these results has been given by M.Olley. | 
oe ST a aes a el 5 esa a lee 
3 M. Olley, Proc. Inst. Automobile Eng., Y 500 1000 1500 lbs 
41 (1946), P- 147. ‘A Fig. 3: The Aligning Torque T as a function of the normal load for different 
dle Bradley and R. F. Allen, Automobile values of slip angles ~ as parameters, from M, Olley [3]. 
ung. 5 BY Das riickdrehende Moment T als Funktion der Normallast fiir verschiedene Werte 
Bes A (Cua Properties of Tires Ef- sae des Schraglaufwinkels « nach M. Olley [3]. 

51D. . 
fecting Ridin Steering and Handling, ; 

SAE-Prans., Feb. 1935. 6 A. W. Bull, J. Soc. Automotive Eng., 44 (1939), p. 344. 


7 W.F. Billingsley, R. D. Evans, W. A. Hulsmit and E. A. Roberts, Rolling Resistance of Pneumatic 
Tires as a Factor in Car Economy, SAE-Trans., Feb. 1942, p. 37. 
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b) Automobile behaviour. Experimental research on the behaviour of an automobile has 
been carried out by 0. Dietz and R. Harling. This problem was given a theoretical treatment 
by P. Rieckert and T. E. Schunck § and additional research has been carried out lately by Eberan 
v. Eberhorst 9 and E. Fiala 2. 

Dietz and Harling reported results of extended experimental work on the subject of auto- 
mobiles in curves. They present various diagrams showing that the amount of skidding is of 
the same order of magnitude as the angles involved in steering. The results show that the amount 
of slipping is very sensitive to the forward velocity of the automobile for a prescribed track. 
This effect is well known by automobile engineers. Dietz and Harling found that tire pressure 
may be treated as a secondary effect. L. Huber pointed out that it is advantageous to have 
the front wheels carry more load than the rear wheels from the standpoint of more easily resisting 
gusts of wind. 

Rieckert and Schunck study the motion of an automobile in a curve by considering a model 
which is simplified by having only one rear wheel and one front wheel mounted in the center 
of each axle. Linearizing the equations of motion of this model and making use of Huber’s 
results for tire characteristics they treat the following two problems: first, the steering angle 
required to follow a given circular track and, second, to find the track when the steering angle 
is prescribed as a function of time. For these two problems the stability criteria are given, 
no comparison with experimental results is included. 

Recently Eberhorst compared the influence on the slip angles of the different ways of distri- 
buting the loads on front and rear wheels and of different ways of connecting the axles to the 
body of the car. 

The scope of this work is to extend the existing theory by treating a more realistic system 
which, contrary to the model treated by Rieckert and Schunck, includes the effects of four wheels. 
As the influence of rolling of an automobile in a curve is not considered in this work we may 
neglect the relative motion between the axles and the body of the automobile. 


3. Nomenclature. The following nomenclature is used in this treatment: 
M = mass of wheels and body, 
i, = radius of gyration of wheels and body with respect to the vertical axis through 
center of gravity of wheels and body, 
I,, = moment of inertia of each wheel about axis through center of gravity of wheel 
perpendicular to plane of wheel, 
I, = moment of inertia of each wheel about axis in the plane of the wheel through 
the center of gravity of the wheel, 
1, 2, 3, 4 = subscripts referring to various wheels (see Figure 4), 
r, O, z = cylindrical ccordinates, 
a; = slip angle of i-th wheel measured between the velocity of the wheel and the 
plance of the wheel in a horizontal plane, 
ry) = distance from center of curve to center of gravity of automobile, 
{2 = angular velocity of radius r,, 
a = distance along axis of car from front axle to center of gravity of body, 
b = distance along axis of car from rear axle to center of gravity of body, 
c = one-half the tread (equal for front and rear), 
Aaja te, B=/jeP+e, 
e = vertical distance of center of gravity above the plane of the axles of the auto- 
mobile, 
jf = radius of each wheel, 
fp = steering angle of front wheels, 
( = swim angle, angle between the long axis of the automobile and the velocity 
of the center of gravity of the body (this angle is sometimes called tke side 
slip angle or the attitude angle) 


8 P. Rieckert and T. E. Schunck, Ing.-Arch. 11 (1940), p. 210. 


9 Eberan v, Eberhorst, Automobiltechn. Zeitschr. 55 (1953 > p. 246; Aut bile Rev il 
10, Fiala, Automobiltechn. Z. 56 (1954), p. 5. : Be iain aH: 
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Vp oe pe ine 
@; = angular velocity of i-th wheel around its axle, 
ec 


~ dé : gers F 
é= Ape first derivative with respect to time, 


V = velocity of center of gravity of automobile, 
V; = velocity of center of gravity of i-th wheel, 


D; = drag force on the i-th wheel, 


F = driving force acting in plane of road on third and fourth wheel, 


P; = vertical load supported by i-th wheel, 


S; = sideforce transmitted from the ground to the i-th wheel orthogonal to the 
plane of the wheel, 


T; = torque transmitted from the ground to the i-th wheel, its vector perpendicula 
to the ground, 


Po> Poo Pi, Fy = constant mean values of the corresponding variables, 
€1> 61» Pi» fy = deviation from the constant mean value of the corresponding variables, 


Che kek = experimentally given constants used to represent the behaviour of the force and 


moment transmitted from the road to the wheel, 
k 


= = ratio of the spring stiffnesses between the body and front axle and the body 
"and the rear axle. 


4. Derivation of Equations of Motion. Since this problem is mechanically very complex we have 
to simplify the given system quite drastically in order to get useful results with a reasonable 
amount of analytical work. 

We assume the body of the automobile to be rigid, neglect the relative motion between the 
axles and the body and consider steering to be effected by the two front wheels only turning 
through the same angle. We furthermore assume that there is no tow-in or camber. This model 
moves on a horizontal flat, circular track without air resistance. Fig. 4eshows schematically the 
system treated. ; 

By neglecting the relative ,,rolling“‘ motion between the axles and the body we imply two 
things: 1) the change in radius of curvature of the center of gravity of the car due to rolling is 
small compared to the radius of the path, and 2) the change in load distribution on the tires due to 
the sidewise motion of the center of gravity of the car is small compared to its mean value. For 
modern automobiles the two implications given above are certainly suitable assumptions for the 
problem treated here. ; 

The equations of motion for our model (Fig. 4) with its center of gravity traveling on a cir- 
cular path with constant magnitude of velocity are 


— M.r, 2 = —(S, + S,) cos (8 + p) — (S; + S,) cos p + (D, + Dz) sin (8 + 9) 
SCD, Deer) sine. (1) 


0 =—(S,+$,) sin (6 +)—(S; + S,) sin p—(D,|+ D,) cos (6 +¢) —(Ds + Di —2 F) cosy, (2) 


Mg=P,+ P,+ Ps + Py, (3) 


Mi?@ + 2 Tew B = S, A cos (f Vad + S, A cos (y,— B) — (Ss + Su) B cos y, — (4) 
A [D,sin (8 +y,)—Dasin (¥p— A) +B (D,—D,) sin yr—(T, + T, + Tr + Th), 


+e) (D, + Dy) cos B—I,, (@ +2) cos B +(5, + S,) (f+) pip [Da O72 P| (ane) 
_ Bee ree bo (p+B+ 2) sin B —(P, + P,) a t (Ps; +4 oe 


i i Do) Si e) cos S,)(f +e) 
Se D, + D,) sin B + Ty (1 + 2) sin B + (Si, + 82) (f+ @) B + (Ss + Su é 
ee oy cen kone OP, Fe, Bend} 


29 
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Since the problem is statically indeterminate one additional equation is required. We make 


use of the following relation: ; 
(Ee ee (7) 


which is derived from the fact that all four wheels must remain in contact with the ground. 
The effect of shock absorbers is not of primary importance and has been neglected in deriving 
Equation (7). 


Fig. 4, Schematic diagram of the model treated. Schema des behandelten Systems. 


From the geometry of the system we may write for the slip angles 


A (p + 2) cos (vp — yp) | 
V—AQ@+2)sn—y,| 
A (p + 2) cos (p + yp) 
Vie a mene 
B(~ + Q) cos p +- y,) (8) 
V+ Bo +2) ee 
B (p + 2) cos (p — y,) 
V+ B@ + Q) sin es) 


a =9 + p— ten | 


% =p + P— tan | 


Ky =P tan-| 


Oy =o+ tan 
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and for the angular velocities of the wheels 

o, = S4/[V 4 AG + Dany — OP HAW +O) cos o,— HF. | 

0, = V[V—A p+ Qin (yy, + PP FLAG FD) c08 (y, + OP: 

= VV + Bp + Q) sin (y, + (E+ [B@ + Q) c08 (y, +9)? 3 


; COS K, 


fF VV — Bp + 2) sin, — @)P + [Bp + Q) 008 (y, — QP - 


(4 


Now making use of the experimental results given by L. Huber we may write the sideforce S; 
transmitted from the ground in the following form: 


i= (ce —c P;) Pra;. (10) 
For the torque JT; transmitted from the road we have the relation 
T, = (k’ +k P) Pa. (11) 


Dealing with the automobile on a horizontal road we assume the rolling resistance and friction of 
the wheel bearing may be neglected so that the drag forces D;, acting on the circumference of the 
wheel, may be written as 
Bye “ Qj - (12) 
Equations (1) through (12) govern the problem we are dealing with. These equations are highly 
non-linear. In order to get a solution of physical meaning with a reasonable amount of analytical 
work we have to simplify the given equations in such a way that the principal characteristics of 
the behavior of the system are preserved. 

The non-linearity of this governing set of equations is not only due to trigonometric terms but 
is also introduced by applying experimental results for the sideforces S; and torques T; trans- 
mitted from the road to the tires. 


5. Approximate Solution. A perturbation method may be applied to get the first approxi- 
mative solution for the non-linear simultaneous set governing the given system. In this pertur- 
bation scheme we split the unknowns into constant mean values and small deviations from these 
amounts that are functions of time. For this purpose we write 


P=Po +> B=Potfi> P, = P; + pi; F=Fy+fi- (13) 


In the Equations (13) the first terms on the right side represent the constant mean values and 
the second terms stand for the variable parts of the unknowns. We further restrict the following 
quantities 11: 


eG) a el Rs OLB. 
Gs max Q ee S;A Po max a Bo max 2 max S;A Bo max Mi P; max By max 


to have values small compared to unity. 

Substituting the expressions (13) into the Equations of motion (1) —(6) and Equation (7) and 
carrying out a standard perturbation procedure we get the set of zero order equations which have 
to be satisfied by the constant mean values of the system. These equations are 


hae ED ie [[e’ == Pe Bs O19 + (c 3 P,) P, oq cos (By + Po) + [(c’ =e Ps) P; X30 | (La) 
=i te ea Pp, P, 4g] COS Po + 2 F'sin gy » 


11Due to the generality of this problem the actual mathematical restrictions to justify our procedure 
are quite complex. The conditions (14) are suitable for this treatment on the average modern passenger vehicle. 
The restrictions are obtained from the requirement that the contributions in the various equations from terms 


including 9, Py» Pv Py By B,, p; and f, are small compared to terms not dependent on the time. 
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OR [(c’ == EP ) Py O19 (c’ —c P,) P, Og» | sin (By + Po) + [(e’ —ec ip. P, X39 | (2a) 
+ (c’ —e OP, P, Ogg] Sin Py — 2 F cos q > 


= # [(e’ == 6 le) P, Oy) £08 (By + yr) + Cie Cc P,} P, Oleg COB (Bo —ve)] 
— B[(c' —€ Ps) Ps og) + (c’ —¢ Py) Py 4] 008 Yr — (ke oe hb abeor (3a) 
— (k’ + k Py) Pe O69) — (k’ +k 2) 1A Oh a (k’ +tkP ye ae 


O=(f+e) [(c’ —@ P,) Py aq + (c' —@ P,) Pha Ota9| sin By — (P, + P,) a+ (P, + P,)b 


Q2 ae = : = 4 A cos (yr — %>) (4a) 
+ Fy f sin By 72+ A? -+ 27, Asin (VF — Go) * 608 |%¥ +6 — tan™ |; + A sin (yp — %) 


0 eee Ft _{ 400s (ve + $0) 
Ly “sin Bo {yr + A®— 2A rysin (yr + Po): 608 | Po + Bo — tan \7, — Asin (yp + H)/II’ 


0 =(f + e) [(e’ =) ig) Py o49 + (¢’ -o) Pp Pay Comey = (ee) [(c’ = Hig 12) Peon 
+ (e' ¢ Pa) Pa oa] 
Q? aT = : (408 Vr = 9) 
+ Tw cos Bo yr? + A? +2 Argsin(yr—p)° 608 |Po + Bo — tan |7, 1 A sin (vp — 0) 


(4.208 (ve + 90) 
cos VFB AP—2 Ar, sim (ye 1) £05 |» + Bo—tan-* (“asin Gre wo) (5a) 


“F 

Q2 ; ; 4 B cos (% + Y,) ) 
Lue /r+ B + 2 Brosin (y, +o) * C08 |Yo + tan ro + B sin (Y + 7) 
2 

a, 


ee _, (_ B08 Go =) 
yr B + 2 Bry sin (Pp —Yr)* COS | Po + tan To + B sin (% —?,) 


+ ¢(P, 4 7 "ig P;), 


ak 


WARE I SR Ses eps Be fe lie (6a) 
Preset (7a) 
P,—P,— 


In the Equations (la)—(7a) the a;) are the mean values of the a; as follows: 


_| 4008 Vr — ¥o) 
O19 = Po + By — tan ro + Asin (yr — %)|’ 
_A.cos (yr + %) 
Arg = Po + Po — tan ro —A sin (yp + %)| ’ 
| [_ Boos (rr + 90) 
Aso = Yo + tan “|r, + Bain (y, +:%)|° 
B cos (Po ie 9) 
Xgq = Po + tan ro +B sin (Po— Y)| 


The eleven Equations (la)—(8a) give the solutions for the mean constant parts 049, M9, A392 X49» 


Pa P,, P;, P,, F,,~) and By. The first order equations derived from the perturbation process 
become quite complicated. As solving of the first order equations of the system is routine work 
once the equations are given and involves a large amount of algebraic and numerical work, the 
equations only are given in the Appendix and a typical problem has been solved below. 


6. Numerical Solutions of Zero Order Equations. In order to allow a comparison of existing 
experimental work with the results predicted by the set of zero order equations the authors have 
carried out numerical calculations for a particular automobile representing a typical combination 
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of the different parameters characterizing the system. The automobile characteristics taken aro 
as follows (see Fig. 4): 


Mg = weight = 4000 Ib = 1810 kg, 


600 — 16 tires (properties given in Figures 5 and 3), 


Ubs 4 
800 sas 
No 
9 Qo 
£00 3° 
7 o 
6 ° 
400' 
5 °° 
Shpargle 
Je 4 
ce] 
200 = 2 
ie 0 
4 = 
load 
Q 500 1000 7500 Lbs 


Fig. 5. The sideforce S as a function of the normal load on the tire for different 
values of slip angels as parameters, from J]. Olley [3]. 


Die Seitenkraft S als Funktion der Normallast des Fahrzeugreifens 
fiir verschiedene Werte des Schraglaufwinkels « nach M. Olley [3]. 


Table 1. Calculated values of zero order solution for chosen numerical example, 
Gerechnete Werte der Ndherungslésung fiir gewdahltes Zahlenbeispiel. 


60 | 60 30 30 30 
18,2 18,2 9,1 9,1 | * 9,1 
250 500 150 250 500 
16 152 4560) 116.0 nls2 
Py 12 |e) 0258 be 00 0s 
A | V7 4,75 | 2,6 1,32 
86 | 2,25 2,55 it 0,49 
8,7 2,25 2,48 | 1,07 | 0,49 
487) 0,72) = 2,05 | —0,97.| 0,435 
Tae 20.79 2,09 1 0,435 
1491 1291 | 1257 | 1191 1140 
| 678 ry am (sy | ee 517 
6Re 89 | 923 | 989 | 1040 
eae ire ries | ern ee 
Tsien, (SioT7) stor 1.900 
595 | 504 | 488 | 459 | 436_ 
/ 509 | 709 | 743 | 809 | 860 
je ae ee 367 | 390 
ee pe is eS ee 
| 50,8 8 6,62 1,75 | 0,41 
tel Oe ata 15 12,5 0 49 Sis ee 8.9911 1.58 0,42 | 0,098 
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I,, = mass moment of inertia for each wheel = 2 lb- ft - sec? = 0,276 kgmsec?, 
oo itt o2m. 

b = 6 ft. = 1,83 m, 

c = 2.5 ft. = 0,764 m, 

e = 1 ft. = 0,305 m, 

f=1.25 ft. = 0382)m, 

Seal 

k 


We made use of the experimental results obtained by M. Olley shown in Fig. 5 and 3. 


degr} 5, Go 700(305) 
60(18,2) 
8 
60(18,2) 
6 
A091) 
4X 
pie) 


4006 1ft 
Ye 


see 


Fig. 6. The mean values of the steering angle f, and the swim angle , as a 

function of the reciprocal of the radius of the curve for various forward velocities. 

B, is given on dotted curves and 9, is given on solid curves. Numbers on curves 
refer to forward velocities in feet per second. 


Die Mittelwerte des Lenkwinkels 8, und des Schwimmwinkels gy, als Funktion 
der Kriimmung der Kurve fiir verschiedene Werte der Fahrgeschwindigkeit. 
B, ist gegeben durch die strichlierten und g, durch die voll ausgezogenen Linien. 
Die eingeklammerten Zahlen an den Kurven geben die Fahrgeschwindigkcit 
in m/sec, 


Calculations were made for different driving conditions, i. e., variable amounts of speed and 
variable radii of curvature. It should be emphasized that we treated the motion of the automobile 
in an unbanked curve. The equations were solved for the unknown mean values of the slip angles 
19> X29 X3q> X4g3 the unknown mean value of the swim angle q,; the unknown mean amplitude of the 


steering angle 6); the mean values P,, P,, P,, P, of the loads of the tires and furthermore for the 
force F’, transmitted from the road to the tire acting in the plane of the road. From the latter 


result for F the horsepower requirement to drive the automobile at the given velocity and on 
the given path was determined. 

The results of this calculation are shown in Table 1 and in addition the most interesting results 
are represented by the curves in Fig. 6 and 7. 
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ihe Comparison of Zero Order Results with Experiments and Conclusions. The results of the 
above calculation allow us to compare our predictions with experimental results available from 
Dietz and Harling. Fig. 8 shows some of the experimental results referred to. The results of Dietz 
and Harling are seen to be in good agreement with the theoretical results stated in Table 1. The 
deviations of the theoretical from the experimental results may be explained by the fact that 
Dietz and Harling used a different model automobile with front wheel drive in carrying out their 
tests than was used in our calculations. 


HP4 Power 


_tott/sec(3n5m/sec) 60 ft/sec (18,2m/sec) 


dott /sec (G7m/sec) 
po ae to 
0 0.002 0.004 006 7/*t 
Fig.7. Horsepower requirement for driving through a horizontal unbanked 


curve as a function of the reciprocal of the radius of the curve @ for various 
forward velocities. 


Leistungsbedarf zum Durchfahren einer horizontalen, nicht iiberhéhten Kurve 
als Funktion der Kriimmung der Kurve fiir verschiedene Werte 
der Fahrgeschwindigkeit. 


The calculation proves that for high velocities and small radii of curvature the demand in 
power to keep the automobile going with a given speed gees up quite strikingly, a fact well known 
to the experienced driver. The steep rise in power with increasing driving velocity and decreasing 
radius clearly points out the fact of large tire wear under these driving conditions. 


Appendix. This section gives the first order equations obtained from the perturbation proce- 
dure used in Equations (1) through (12). As the authors point out in the text, solving of these 
equations generally is routine but quite involved algebraically. These equations are included as 
they are certainly of interest and to the authors’ knowledge have not been presented yet. The slip 
angles a; and the velocities w; of the wheels can be written in the following form: 


Oy = +P, +%% +0, ¥%, 

Oy =O +61 +% 2, +7 Pe, 

Os = Os +919; + 1 Vs: (1A) 
Oy = Oy +92, + (1 Yi» 

re TE OK, acy Le. (EA 


The Equations (1A) show that by linearization we are able to split the amount of the slip angles 
a; as well as of the angular velocities @; into parts due to the zero order solution, the steering 
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degrh p 


0 
50 60 70 80 90 ft/sec 


Fig. 8. Steering angle f and swim angle @ as a function of the forward velocity V 
for constant radius of curvature 9 ~550ft., from Dietz and Harling [2]. 


dotted curves — ordinary road tires, tire pressure = 28 psi, solid curves — snow 
tires, tire pressure = 28 psi, dash-dotted curves — snow tires, tire pressure = 35 psi. 


Lenkwinkel 8 und Schwimmwinkel y als Funktion der Fahrgeschwindigkeit V’ 
fiir konstanten Kriimmungsradius g ~ 170m nach Dietz und Harling [2]. 


strichlierte Linie — gewéhnliche Strafienreifen, Reifendruck ~ 2 atii, voll 
ausgezogene Linie — Schneereifen, Reifendruck ~ 2 atii, strich-punktierte 
Linie — Schneereifen, Reifendruck =~ 2,5 atii. 


angle B,, the swim angle ¢,, and its first derivate y,. In the Relations (1A) the factors ®;, W;, H;, 
K; and L; are functions of the given parameters of the system and the zero order solution. Using 
the above abbreviations we get the following set as the first order equations: 


O0=+[(e’—e 1) Prag + (c’ —€ Ps) Pyoing] sim (By + Po) (Bx +9) + [(e’ —€ Py) Py so 
aie (c’ —¢ P,) P, O49] (sin Go) + Y, — 2 F (cos Hp) - p1 — 2 (sing) ft 


red Bi 2 : 
2 Ff Sin (Bo + Yo) (Ki + Ke) + Gy + (Ly + L,) 91] 


Sia - sin @y [(K, + Ky) g, + (LZ, + Ly) 9] 

ec OS (Bo + Po) {(c' —26 Pr) or Ps +(c’—2¢ P.) dg Po ae [(c’ —e P,) ie 

+ (c’ — © P,) P,] B,} — cos (By + 9) {[(e’ —@ Br) Py , + (c’ —e P,) P, Oy] gy 
+ [(c'—@ PE. bai (c’—¢P,) P,P.] G1} — cos po {(c'—2¢ P,) X30 Pa 

+(c’—2e Ps) O49 Pa + [(e’ Gi) P; ®; + (c’ Sate) P, P,) 91) 

— cos Y[(c’ —€ P;) P, P+ (c’ —e P,) P, Y,| V1» 
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O— |{c'—e P,) Pion, + (c’ Se P,) P Py A459] cos ( (By + Go) * (Bx +o) + [le — Py) Ps ds9 
qe =¢ P,) P, O49] (cos Yo) - PRONE) #1 —2 (008 gr) fi + Fc08 (By + Po) 
x [(K 1 + K,)g, + (L 1 + Ly) gi] + 608 gy [(Ke + Ky) p, + (Lz + Ly) @1] 


+ sin ( (Bo + Po), {(c’ — 26 Py) ayy py +(e’ — 2€ P») bo pp + [(e" See (3a) 
+(e’ —¢ P,) Ps] By} + sin (By + 0) {[(’ —@ P,) P, B, + (/ —€ P,) P, ,] 9, 

+ [(e —e P,) P, WY, + (c' —¢ P,) P, Po] Gi} + sin yp {{c' —2¢€ ae X39 Ps 

+ (c’ —2€ P,) oy ps, + [(c’ —é P,) P, ®, + (c’ —¢ P,) P, D,] 
+ sin po [(c’ —¢ P,) P, VY, + (c’ —¢ P,) P, Yo; 


Se 


Miz @, + 2 Ino B, = — A (c’ —€ P,) Py o49 sin (By + yr) Pr + A[(e’ —2 ¢ Py) O19 Px 
+(e! —eP,) P, (Ai +919, +9, ¥. 1)] ¢08 (By +7) — A(c’ —€ P,) P,o%9q sin (Bo —yr)B 
+A [(e’ oe P,) a Xs Po 1 (c’ — P,) P, (Bi + 1 @. + G1 P2)] cos ( (By — yr) 
— Boos y, {(c — Je TP.) aay Pee (e Lips taie —= ie) P, DB; 
+ (¢ —€P,) PD], + [e' —FP) PY +(e 7B) AY) A)—F A 
x {LK, sin (By + YF) ++ K, sin (By) —yr)] G1 + [Li sin (By + yr) + Lz sin (Bo — vr) | ¢1} (4a) 
+3 BI[Ks— K;) sin y, Gy Ae (Ti Ty einy = pl k’ LO. Pi 
kb +2 hk P,)o Bop Pah 2b P,) ogo pe — (FE 2k P,) oy Ps 


—[(k’ +k P,) Py + (k’ +k P,) P,] By —[(k’ +k P,) P.O, + (W +k P,) P,® 
+ (k’ +k P,) P, ®, + (k’ +k P,) P.O.) o,—[(k +h Py) PY 4( +h P,) PLY 
ee Pee, (hk! bP) P, P| 9%, 
0 = wp (By + Ky) cos By: g, + (Ly nad as aie 
+ (ce —¢ P,) P 2 Xoo] (f + e) cos By- Py + e) sin By {(c’ —2€ Py) O19 Pi 
+(e 6" Ps) Gog Ps + tie — € P,) P, ee ¢ P,) P,| Bi} + (f + 2) sin Bo 
xAe ane bP D, +(e —¢ P,) P, PD.| gp, + [(c’ —e P,) P, le ores Polo} (5a) 
Ss lw [(Ks + Ks) ¢ @y + (Ls + Ls) Qi] + Lo Q sin By [(K, + K,) 9, + (ZL, + L,) o] 


+ I, (Hy + Hz) 2 (cos Bo) By + 1, sin Bo (Aa + £ ig) (G1 + Bx) —2(f+¢h 
— (p; + Ps) @ + (Ps + Pa) 5 


O = — Iw -¢ sin Bo [( [(K, + Ke) o, + (Lr + Ly) Ga] —(f + 2) sin By [(c’ — ¢ Py) Py or 
1 @—7F) Fimdl+04 ¢) 008 By {(c’ — 26 Py) ayo Py +(e’ —2¢ Ps) ony P2 
+ [(c’ — eP,) Py + (c’ —€ Ps) Po] By} +(f+e) ) cos By {[(c’ —¢ P,) P, D, 
4 (c’ —6 P) Py ®,] p, + [(e’ —¢ Py) Pi + (c --6 Pe) P,P a} + (f+ 4) (6a) 
x {(c’ — 2¢ Ps) digo Ps + (¢' — 26 Py) O49 Pa + [(c’ —¢ Ps) Ps Ds +(e 0 P,) P, D,| 9; 


+ [(c’ —eé Pp) P; Y% +(e’ —e P,) P,P] G1} + Lo {— (Hy + Ay) @ sin By: By 
+ (H, + H,) cos Poe (4 oa bx) + 2 cos By [(Ky + K,) g, + (Li + £2) Pil} 
+ Io {(H; + Hy) gy + 2 [(K; + K,) gi + (43 + I,) fil} +- € (— py + P2— Ps Pas 
Pit Pe t+ Ps + Ps = 9; 7a) 


k 
Bees Py aaa (8a) 
Ps — P3 k, 
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One important use of this set of equations should be mentioned: the proof of the stability of 
the system. The zero order solution given in the previous section is only valid if the system is 
stable, this fact being up to now purely a matter of experience. 
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